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1.0  Introduction 


The  turbulent  flow  around  a  wing-body  junction  is  a  very  complex,  highly  three 
dimensional  flow,  yet  it  is  commonly  encountered.  Figure  1.1  shows  a  simpified  wing- 
body  junction  flow.  Note  that  the  tunnel  coordinate  system  is  shown,  which  will  be  used 
to  present  the  majority  of  the  data.  The  dominating  flow  characteristic  is  the  mean 
secondary  flow  strucuare  known  as  the  horseshoe  (or  necklace)  vortex.  Near  the 
appendage  nose,  the  flow  is  highly  unsteady,  and  velocity  histograms  show  bimodal 
behavior,  termed  the  zero-flow  and  backflow  modes  (Devenpon  and  Simpson,  1990b). 
Flow  visualizations  performed  in  this  region  using  hj  drogen  bubbles  in  water  (Kim,  1991) 
show  a  very  complex,  highly  turbulent  flow,  where  the  instantaneous  ve’ocity  distribution 
is  often  very  unlike  the  mean  flow  patterns  associated  with  either  flow  mode.  Tne  flow 
field  resulting  as  this  unsteady  vortical  structure  is  stretched  and  skewed  around  the  body 
is  complex  and  not  well  understood.  The  author’s  study  attempts  to  gain  insight  into  the 
physical  nature  of  the  junction  vortex  flow  by  examining  and  comparing  his 


Introduction 


1 


measurements  to  those  of  other  researchers  with  the  same  body  shape,  but  varying 
approach  boundary  layer  flow  conditions. 

The  horseshoe  vortex  is  aeated  by  a  combination  of  two  effects,  the  skewing  and 
stretching  of  the  transverse  vorticity  present  in  the  incoming  turbulent  boundaiy  layer  as 
it  passes  the  wing,  and  the  rolling  up  of  the  incoming  fluid  along  the  centerline  due  to 
the  large  adverse  pressure  gradient  created  by  the  wing  geometry.  The  term  horseshoe 
vortcx  is  a  misleading  description,  since  it  is  not  a  typical  "tornado "  type  vortex.  The 
mean  vordcity  distribution  tends  to  be  flattened  into  a  more  elliptical  shape,  and  the 
majority  nf  the  streamwise  vorticity  (ii,),  except  near  the  appendage  surface,  is 
contributed  by  the  transverse  skewing  term,  9W/9Y.  Dickinson’s  (1986a,b)  work  in 
^^L^culitr  shows  those  results  to  be  true. 


McMahon  et  al.  (1987)  discusses  the  two  general  types  of  secondary  flov  ,  first 
described  by  Prandtl,  flows  of  the  "first  kind"  and  "second  kind".  Secondary  flows  of  the 
first  kind  arc  caused  by  a  shear  layer  which  becomes  skewed  about  a  streamwise  axis. 
The  subsequent  transport  of  vorticity  into  tite  streamwise  direction  results  in  the 
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cur/ature  of  the  streamlines  and  resulting  pressure  gradients  skew  the  shear  layer  and 
stretch  tire  resulting  streamwise  vorticity.  Secondary  flows  of  the  second  kind  arc  caused 
by  Reynolds  stress  gradients  in  planes  normal  to  the  streamwise  direction.  This  type  of 
flow  occurs  in  comers  created  by  intersecting  semi-infinite  planes,  and  is  described  by 
Gessner  (1973). 
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Tiiere  are  many  initiatives  for  researching  junction  flows.  Junction  flows  are 
commonly  present  at  aircraft  wing  roots,  at  ship  and  submarine  appendage-hull  junctions, 
at  bases  of  bridge  supports,  in  turbo-machinery  flows,  and  also  in  flows  with  injection, 
or  a  component  of  the  injection,  normal  to  the  main  flow.  The  vortical  flow  structure  can 
cause  problems  such  as  scouring  and  deterioration  of  support  foundations,  excessive  heat 
transfer,  and  excessive  noise  and  vibration  due  to  the  creation  of  additional  turbulence  in 
the  junction.  Further  research  leads  to  2  greater  understanding  of  the  physical  nature  of 
juncture  flows,  and  also,  vortical  and  3-D  turbulent  shear  flows  in  general,  enabling 
control  of  their  adverse  effects.  Also,  additional  experimental  work  gives  CFD 
researchers  a  larger  data  base  for  comparison  of  numerical  results. 

The  author’s  investigation  of  a  wing-body  junction  flow  was  performed  by  taking 
U,  W,  u',  w\  and  -iiw  hot-wire  measurements  in  6  YZ  planes  adjacent  to  and  downstream 
of  an  appendage-flat  plate  junction  (Figure  1.1),  These  planes  are  designated  5  through 
11.  The  locations  of  planes  12-15  are  given  in  Table  1.1.  The  appendage  geometry  was 
a  3:2  elliptical  nose  joined  at  its  maximum  thickness  to  a  NACA  (X)20  section  which 
formed  the  tail.  A  total  of  64  profiles  of  approximately  30  points  each  were  taken  using 
a  boundary  layer  type  single  hot-wire  probe.  Much  farther  downstream,  U  and  u  data 
were  obtained  in  3  YZ  planes  (planes  13-15)  located  in  the  far- wake  region,  again  using 
a  single  hot-wire  probe.  A  total  of  45  profiles  were  taken  in  these  planes.  Two  of  these 
far  wake  planes  were  located  in  a  diverging  section  of  the  wind  tunnel,  in  order  to  study 
the  effects  of  an  advene  pressure  gradient  on  the  trailing  legs  of  the  horseshoe  vortex. 
This  configuration,  described  fully  in  section  2.2,  simulates  the  conditions  over  the  aft 
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portion  of  a  body  with  a  tapering  surface,  such  as  an  aircraft  fuselage  or  the  hull  of  a 
submersible. 

The  main  emphasis  of  this  study,  however,  is  the  comparison  of  data  from  junction 
flow  experiments  involving  the  same  body  shape,  the  aforementioned  3:2  elliptic  nose 
with  a  NACA  0020  tail.  The  purpose  of  the  comparisons  is  to  determine  the  approach 
boundary  layer  effects  on  the  junction  flow.  Section  1.1.2  describes  these  experiments, 
and  section  1.2  outlines  the  author’s  research  goals.  The  literature  review  focuses  on 
junction  flow  experiments  which  examined  the  effects  of  varying  the  approach  boundary 
layer,  since  these  studies  arc  most  relevant  to  the  author’s  work. 


1.1  Review  of  Previous  Work 


A  great  amount  of  previous  work  has  been  done  researching  the  flow  in  a 
turbulent  wing -body  junction.  Past  research  has  involved  various  measurement  methods 


and  also  vanous  ajj^enuage  geomcuics.  Juncture  ilow  studies  can  be  divided  into  two 


broad  categories  based  on  geometry,  1)  flov/s  around  a  cylinder  mounted  normal  to  a  flat 
plate,  and  2)  flows  around  a  more  streamlined  or  airfoil  shape  mounted  normal  to  a  flat 
plate. 
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I.IA  Cylindrical  Appendage -Flat  Plate  Juncture  Flow  Studies 

Measurements  for  these  flows  are  usually  taken  only  upstream  of  the  maximum 
thickness  due  to  the  large  scale  separation  and  unsteadiness  present  in  the  large  wake 
region.  These  experiments  are  useful  in  studying  the  formation  of  the  horseshoe  vortex, 
and  the  flow  structure  in  the  steady  flow  region  is  similiar  to  that  in  the  nose  regions  of 
more  streamlined  shapes. 

Baker  (1980)  studied  the  changes  in  the  horseshoe  vortex  structure  at  the  base  of 
a  cylinder  due  to  variations  in  the  incoming  boundary  layer.  Using  oil-flow  visualization, 
he  examined  the  change  in  distances  Xs  and  Xy  as  Rcg,  and  6*/D  changed.  Xj  is  the 
location  of  the  saddle  point,  or  primary  separation,  and  Xy  is  tlic  location  of  the 
secondar>'  separation,  which  indicates  the  main  vonex  position.  Re^*  is  the  Reynolds 
number  based  on  the  approach  boundary  layer  displacement  thickness  (5*),  and  D  is  the 
cylinder  diameter.  Figure  1.2  shows  sketches  of  the  oil- flow  pattern  observed,  and  the 
significance  of  Xg  and  Xy.  Baker  conjectured  that  a  4-vortex  mean  flow  model  was 
applicable  for  this  pattern,  also  shown  in  Figure  1.2.  By  measuring  mean  pressure  along 
the  flow  centerline,  he  found  a  local  minimum  in  pressure  which  corresponded  to  the 
primary  vonex  location,  Xy  (see  Figure  5.3).  From  these  mean  pressure  measurements. 
Baker  postulated  that  changes  in  Res*  or  S'/TD  do  not  change  the  venical  dimension  of 
the  vonex  system  in  the  plane  of  symmetry.  From  the  oil-flow  results,  he  concluded  that 
the  4-vonex  mean  flow  model  existed  along  the  centerline  for  4,000  <  Reg,  <  90,000  and 
.0333  <  5‘/D  <  .25. 
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Later  studies  by  Baker  (1985)  attempt  to  further  quantify  the  effects  of  Res.  and 
5*/D  on  the  locations  of  minimum  and  maximum  corresponding  to  Xg  and  Xy 
respectively.  Using  Ids  data  as  well  as  data  from  Ram  (1963),  Langston  and  Boyle 
(1982),  and  Peake  and  Galway  (1965),  he  attempted  to  correlate  Xy/D,  X^/D  =  f(Res., 
6*/D)  for  laminar  and  turbulent  boundary  layers.  For  turbulent  boundary  layers,  his 
results  were 

Xy/D  =  0.100(Re5.)  ‘  For  1,000  <  Reg.  <11,000 

Xg/D  =  0.373(Re5.)°^^  For  1,500  <  Rcs,  <  10,000 

Xg/D  =  15.6(Re6.)-^'  For  10,000  <  Reg.  <  18,000 

where  Xy  and  Xg  are  measured  from  the  cylinder  leading  edge.  Note  the  lack  of 
dependence  on  5’/D.  Baker  (1985)  also  presents  a  theoretical  investigation  supporting 
the  lack  of  dependence  on  5*/D.  Note  that  the  author  has  included  the  references  for 
Ram  and  Peake  and  Galway  for  informational  purposes  only,  as  he  was  unable  to  obtain 
copies  of  these  references. 

Belik  (1973)  performed  research  similiar  to  Baker’s.  The  goal  of  his  experimental 
work  was  to  correlate  Xg  /D,  found  using  oil-flow  visualization,  for  a  cylinder-end  wall 
junction  flow,  as  functions  in  the  form  ffRe^,  Shy)  or  f(Shy,  Foy).  Shy  is  the  vortex 
Strouhal  number,  defined  by  Belik  as  U,ef/(DQ),  where  Q  in  this  case  is  the  vorticity  of 
the  horseshoe  vonex  system.  He  assumed  H  to  be  constant  in  a  plane  perpindicular  to  the 
cylinder’s  surface.  FOy  is  the  vortex  Fourier  number,  defined  as  vD/(U,^f)^.  Belik 
assumes  that  -  cOx  -  (3(U)/^Y),^J.,  For  a  typical  turbulent  boundary  layer,  he  assumes 
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(3(U)/dY)„^  =  .5(U„f/8).  After  substitution,  vShy  =  25/D,  and  Foy  =  .5(v/U„f5).  It  can 
then  l-e  seen  that  Foy  ~  Rcg"’. 

Belik  also  found  Xj  /D  for  a  variety  of  artificially  produced  turbulent  boundary 
laycre  with  differing  (5(U„f)/5Y),ve  values,  and  theiefore  different  values  of  Foy.  He 
reached  the  final  functional  relationship  of  Xg  /D  =  f(Fov,  Rcq)  for  the  range  36,tX)0  < 
Reo  <  220,000  and  .4  x  10"^  <  Foy  <  3.5  x  10■^  The  function  correlating  the  data  is 
Xs/D  =  (10-*/FOy)'%516  4-  .lll(lO-^ReD))  (a  =  .016) 

This  is  a  very  weak  function  of  RCp.  The  approximation  to  Xj/D  =  f(Fov)  is 
Xs  /D  =  .57(10-^/Fov)  *’  (CT  -  .0202) 

One  can  see  that  Xj/D  «  (Rcg)  '^.  This  is  the  same  type  of  functional  relationship  which 
Baker  later  observed. 

Belik  also  measured  mean  static  pressures  at  the  wall  along  the  centerline.  For 
the  cases  presented,  the  local  minimum  in  Cp  corresponding  to  the  vortex  core  location 
was  relatively  constant,  remaining  at  Xy/D  =  .21  (again  comparable  to  Baker’s  results). 
The  oil-flow  visualizations  shown  seem  to  indicate  a  vortex  system  similiar  to  Baker’s, 
linforiunaieiy,  Beiik  does  not  give  details  of  the  incoming  boundary  layer  parameters  such 
as  >3  or  5*  for  a  more  direct  comparison  to  Baker’s  work.  The  estimated  range  of  8/D 
is  .('/97  to  1.48,  and  from  this,  an  estimate  of  the  5*/D  range  would  be  approximately  .015 
<  67D  <  .22.  TInis  co’TCsponds  to  a  range  of  540  to  48,400. 

Eckeile  and  La.ngston  (1987)  also  performed  studies  on  a  junction  flow  around  a 
circular  cylinder.  They  performed  surface  static  pressure  measurements,  surface  flow 
visualizations,  and  mean  velocity  measurements  with  a  five-hole  prolie.  Their  test 
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conditions  wc^:.  =  30.5  m/s,  Rcp  =  550,000,  Rcg*  =  8,470,  RCq  -  6,270,  6/D  =  .099, 

and  6*/D  =  .0154,  where  D  =  29.8  cm.  Tlie  saddle  point  location,  Xs/D,  for  these  flow 
conditions  was  approximately  .4,  while  Xy/D  was  approximately  .22.  Baker’s  correlation 
gives  values  of  .503  and  .25  respectively.  Eckerle  and  Langston  concluded  that  for  their 
case,  the  single  primary  vortex  model  applied,  and  that  the  horseshoe  vortex  structure 
became  fully  formed  between  the  5°  and  25°  planes  (centerline  plane  at  0°).  Note  that 
dieir  vorte.'  mci^el  conclusion  does  nut  violate  Baker’s  earlier  findings,  as  5*/D  =  .0154, 
which  is  outside  the  range  given  by  Baker. 

LaPleur  find  Langston  (1990)  present  an  intriguing  study  of  an  iceformation  design 
of  a  hull  near  a  cylinder/hull  juncture.  The  geometrical  features  of  the  resulting 
icefonnation  are  directly  related  to  flow  and  themul  boundary  conditions.  The  test 
parameters  varied  were  0^.  =  (Tp  -  To)/(Tg  -  Tr)  and  Rcu,  where  Tp  =  temperature  of  fluid. 
To  =  freezing  tcisiperature,  and  Tj^  =  temperature  of  boundary  surface.  The  Rcp,  range  was 
737  to  1,84'3. 

LaReur  and  Langston  present  sketches  of  their  results,  showing  separations  as 
ridges  in  die  ice,  aid  valleys  being  formed  by  reattaching  flow.  For  Bj-  fixed  at  .37,  the 
Rcd  =  737  case  resulted  in  only  1  visible  vonex  in  the  plane  of  symmetry,  while  die  other 
cases  produced  2  predominate  vorticies  (Baker’s  4  vortex  model).  Finally,  they  found 
when  using  an  iceform  contour  in  the  junction,  a  5  to  6%  drag  reduction  over  the  flat- 
plate  case  could  be  achieved. 

Agui  and  Anureopoulos’  (1990)  study  of  the  flow  around  a  cylinder  was  concerned 
with  the  mean  and  fluctuating  wall  static  prcs.sures.  Their  main  goal  was  an  increased 
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physical  understanding  of  junction  flows.  They  were  hoping  to  gain  information  about 
the  vorticity  flux  out  of  the  wall  by  measuring  the  wall  static  pressures.  They  describe 
the  relationship  between  the  wall  vorticity  flux  and  wall  pressure  in  a  later  paper  (1991). 
From  the  momentum  equation  evaluated  at  the  wall,  the  equations  are 

d(p)  _  _ 

ax  ■  ar  d.i, 

dZ  ^  dr 

These  relationships  show  the  importance  of  die  wall  pressure  gradient  caused  by  the  body. 

Two  cylinder  junction  flow  cases  aie  presented,  Re^  =  1.0  x  10^  and  2.2  x  10^ 
(Agui  and  Andreopoulos,  1990).  The  given  Xj/D  values  for  these  cases  are  ,26  and  .32, 
respectively.  The  time  averaged  flow  structure  was  taken  to  be  represented  by  the  2 
vortex  model  (1  main  vortex).  An  estimate  of  the  boundary  layer  thickness  is  given  as 
6/D  =  .1.  It  is  not  clear  if  this  is  constant  for  both  Reynolds  numbers,  and  one  would 
think  that  6  would  decrease  as  Reo  increased.  However  6/D  =  .1  gives  Rcs.  estimates 
of  1,500  and  3,3(X)  for  675  ratios  of  .15,  and  5*/D  =  .015,  which  shows  increasing 
as  Reg*  increases,  agreeing  with  Baker’s  and  Belik’s  results. 

I.IJ  Streamlined  Appendage-Flat  Plate  Juncture  Flow  Studies 

The  past  decade  has  seen  much  research  concerning  the  juncture  flow  around  a 
Streamlined  appendage.  The  pressure  gradients  around  the  aft  portion  of  most  streamlined 
shapes,  unlike  a  cylinder,  are  usually  mild  enough  to  delay  separation  until  the  trailing 
edge  region.  As  oil-flow  patterns  show  (Figure  3.10a),  most  streamlined  appendage 
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juncture  flows  have  a  distinctive  flsh-tai!  wall-limiting  streamline  pattern,  resulting  from 
separation  at  the  trailing  edge.  Because  of  the  lack  of  imge  scale  separation,  and 
corresponding  unsteadiness,  measurements  can  be  made  adjacent  to  and  in  the  wake  of 
the  appendage. 

Dickinson  (1986a, b)  presents  a  juncture  flow  study  on  a  body  with  a  3:2  elliptic 
nose  joined  to  a  0020  tail  section  at  their  respective  maximum  thicknesses.  He  first  tested 
a  (X)20  appendage  shape,  but  he  reports  that  this  more  streamlined  shape  did  not  produce 
crossflow  velocities  as  great  as  the  modified  0020  appendage.  The  larger  crossflow 
velocities  caused  by  the  more  rounded,  blunter  nose  enables  the  trailing  horseshoe  vortex 
structure  to  be  more  easily  identified. 

Devenport  et  ai.  (1990)  examined  the  horseshoe  vortex  legs  produced  by  the 
modified  NACA  0020  geometry  in  a  YZ  plane  at  XJC  =  3.00.  Two  junction  geometries 
were  used:  the  unmodified  junction  (appendage  surface  normal  to  floor  at  wing-body 
intersection),  and  a  junction  wrapped  by  a  fillet  of  radius  3.81  cm  (.53  T).  Two  different 
approach  boundary  layers  were  produced.  The  thick  boundary  layer  case  is  identical  to 
the  author’s  approach  flow.  The  thin  boundary  layer  case,  created  using  a  false  floor,  was 
approximately  50%  as  thick  as  the  fomier  case  (see  Table  3.3).  Devenpon  et  ai.  mainly 
studied  the  effects  of  adding  the  fillet  to  the  appendage  base.  They  found  that  the  vonex 
legs  were  further  apart  and  that  the  region  of  boundary  layer  distortion  was  larger  when 
using  the  fillet.  They  concluded  that  the  fillet’s  primary  consequence  was  to  increase  the 
effective  appendage  nose  radius.  It  was  also  seen  that  the  horseshoe  vortex  was  only 
slightly  effected  by  the  changing  approach  boundary  layer  thicknesses. 
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In  addition  to  the  above  studies,  many  other  experiments  have  been  performed 
using  the  modified  NACA  0020  appendage  shape,  which  has  become  somewhat  of  a 
standard  test  case  within  the  past  decade.  Aside  from  the  studies  done  using  it  witiiin  the 
Aerospace  and  Ocean  Engineering  Department  at  Virginia  Tech,  measurements  were  done 
on  this  shape  in  the  Mechanical  Engineering  Department  at  Virginia  Tech  (Shin,  1989). 
McMahon  et  al.  (1987)  also  performed  junction  flow  studies  using  the  modified  NACA 
0020  geometry.  Table  1.1  lists  the  locations  and  types  of  measurements  taken  in  juncture 
flows  with  this  same  geometry.  Note  that  Table  1.1  includes  measurements  made  in  YZ 
planes  only,  which  are  comparable  to  the  author’s  measurements.  Complete  data  sets 
from  Dickinson  (1986b),  McMahon  et  al.  (1987),  Devenport  and  Simpson  (1990a),  and 
Devenport  ei  al.  (1990)  were  made  available  to  the  author  on  disk,  while  tabular  data  and 
figures  from  Shin  (1989)  were  available.  Selected  profiles  of  Shin’s  data,  with  Z/T 
locations  matching  the  author’s  locations  almost  exactly,  were  entered  into  computer  data 
files  for  plotting  and  numerical  calculations.  The  availablity  of  data  from  the  other 
researchers  made  comparison  of  their  experimental  results  possible. 

Oilier  research  done  on  the  junction  flow  produced  by  the  modified  NACA  0020 
appendage  shape  includes  upstream  boundary  layer  profiles  along  a  streamline,  showing 
die  pressure  induced  skewing  of  the  incoming  flow  (Olcmen,  1990),  and  skin-friction 
measurements  made  with  an  oil-flow  laser  interferometer  (Ailinger,  1990).  Shin  (1989) 
reports  on  the  effects  of  various  short,  strake-like  leading  edge  fairing  geometries  added 
to  the  modified  NACA  0020  appendage.  Devenport  and  Simpson  (1990b)  made  LDV 
measurements  in  the  plane  of  symmetry  at  the  appendage  nose.  These  results,  as 
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mentioned  at  the  beginning  of  the  chapter,  clearly  show  the  bimodal  characteristics  of  the 
flow  in  the  nose  region.  The  wind  tunnel  and  inflow  boundary  layer  parameters  for 
Devenport  and  Simpson’s  LDV  measurements  were  the  same  as  for  the  author’s  work. 
Rife  (1991)  studied  the  relationship  between  velocity  and  surface  pressure  fluctuations  in 
the  bimodal  flow  region  using  an  LDV  and  microphones.  Kim  (1991)  also  made 
measurements  on  the  same  appendage  shape,  but  in  a  water  tunnel.  Using  an  LDV,  he 
made  measurements  in  the  plane  of  symmetry  at  two  different  Reynolds  numbers.  His 
data  also  show  the  bimodal  characteristics  of  the  flow  quite  clearly. 

Many  other  experimenters  have  researched  wing-body  flows  around  various 
streamlined  appendages.  Mehta  (1984)  studied  the  effects  that  the  appendage  nose  shape 
had  oil  the  transport  of  vorticity  in  the  strcaniwibe  direction.  Shin  (1989)  also  examined 
the  effects  of  nose  geometry  on  the  horseshoe,  vortex  mean  flow  structure.  The  results 
of  these  researchers  are  used  in  chapter  6. 

Kubendran  et  al.  (1986)  examined  the  turbulent  juncture  flow  around  a  constant 
thickness  body  with  a  3;2  elliptic  nose.  They  documented  the  mean  vortex  path,  and 
noted  similarities  between  the  turbulent  shear  stresses  and  the  mean  flow  strain  rates. 
Similarities  between  profiles  of  u  and  w'  were  also  found.  Shabaka  and  Bradshaw’s 
(1981)  measurements  in  a  wing-body  junction  revealed  that  the  eddy  viscosities  for  the  - 
uv  and  -uw  Reynolds  stresses  are  anisotropic  and  are  negative  over  large  regions.  They 
concluded  that  calculation  methods  based  on  eddy  viscosities  are  likely  to  be  unsuitable 
for  juncture  flows. 
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WcMxl  and  Westphal  (1991)  present  measurements  taken  in  a  lifting  wing-body 
junction.  They  report  that  the  trailing  suction-surface  leg  is  the  stronger  of  the  two 
trailing  vorticies.  Since  streamlines  passing  over  the  suction  side  of  the  appendage  are 
more  distorted,  this  is  an  intuitive  result.  Moore  and  Forlini  examined  the  junction  flow 
around  a  Rankine  half-body,  a  geometry  with  a  well  known  potential  flow  solution.  Their 
measurements,  which  include  mean  velocity  results  in  4  cross-stream  planes,  were 
primarily  intended  for  validation  of  numerical  calculations. 

Rood  (1984)  investigated  the  temporal  flow  features  inherent  in  a  wing-body 
junction  flow.  He  found  distinct  flow  structures  in  the  approach  boundary  layer,  in  the 
wing-body  horseshoe  vonex,  and  in  the  junction  vortex-wake  flow.  These  flow  structures 
were  all  due  to  the  piesence  of  the  wing,  one  resulting  from  the  distortion  of  pre-existing 
structures  in  the  approach  boundary  layer,  and  the  others  newly  created  by  the  wing. 
Rood  reports  that  the  structures’  characteristic  frequencies  were  at  least  an  order  of 
magnitude  lower  than  the  Strouhal-type  shedding  from  the  wing  trailing  edge. 

Several  researchers  have  performed  experiments  involving  vonices  embedded  in 
turbulent  boundary  layers.  Cutler  and  Bradshaw  (1986)  studied  the  pair  of  trailing 
vortices  generated  by  a  delta  wing  along  a  flat  plate,  simulating  a  close-coupled  canard 
and  aircraft  wing.  Westphal  et  al.  (1985)  examined  the  behavior  of  a  vortex  embedded 
in  a  turbulent  boundary  layer  subjected  to  an  adverse  pressure  gradient,  using  zero 
pressure  gradient  results  for  comparison.  Shizawa  and  Eaton  (1990)  investigated  the 
interaction  of  a  vortex  embedded  in  a  skewed  turbulent  boundary  layer,  a  configuration 
siiniliar  to  the  trailing  edge  region  of  a  streamlined  appendage-flat  plate  junction. 
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Westphal  et  al.  and  Shizawa  and  Eaton  both  used  half-delta  wings  protruding  from  the 
flat  plate  to  generate  the  embedded  vortices.  Even  though  the  vorticity  magnitudes 
encountered  and  distribution  patterns  produced  by  these  flows  are  somewhat  different  than 
encountered  in  a  junction  vortex,  the  results  of  these  studies  are  relevant  to  the  author’s 
work  and  are  discussed  in  chapter  4. 

1.1.3  Numerical  Work  in  Appendage-Junction  Flows 

Some  numerical  work  has  been  perfonned  on  the  juncture  flow  problem  for  both 
laminar  and  turbulent  flows.  While  accurate  soludons  for  the  laminar  cases  are  relatively 
common,  the  turbulent  flow  solutions  suffer  from  the  lack  of  an  adequate  3-D  turbulence 
model. 

Visbail  (1991)  studied  the  laminar  juncture  flow  around  a  circular  cylinder.  He 
obtained  solutions  for  6/D  =  .1  and  Rep  values  in  the  range  of  500  to  5,400.  He  observed 
a  non-unique  relationship  between  wall  limiting  streamline  patterns  and  the  flow  above 
the  plate.  He  also  observed  that  the  flow  becarne  unsteady  above  Rcd  =  4,000,  with  the 
resulting  Strouhal  number  (defined  as  /D/U„f)  equal  to  .21.  This  result  agrees  quite  well 
with  Sh  =  .20,  found  by  experiment  (referenced  in  Visbail,  1991).  Visbail  reports  that 
this  unsteady  behavior  is  not  due  to  Karman  vonex  shedding  in  the  cylinder’s  wake. 

Deng  et  al.  (1990),  performed  calculations  on  the  turbulent  flow  in  the  junction 
around  an  appendage  with  the  same  geometry  studied  by  the  author  and  several  others. 
The  inflow  boundary  layer  conditions  matched  Dickinson’s  thick  boundtu^  layer  test  case 
(1986a,b,  see  Table  3.3).  The  streamwise  velocity  contours  compare  well  qualitatively 
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with  experimental  results,  but  quantitatively,  the  calculated  isovels  are  consistently  less 
than  the  experimental  results.  This  may  be  due  to  appendage  blockage  effects  present  in 
Dickinson’s  test.  These  effects  are  discussed  in  section  4.4.  No  turbulent  fluctuation 
results  are  presented. 

Sung  and  Lin  (1988)  obtained  solutions  to  the  Navier-Stokes  equations  using  the 
Baldwin-Lomax  turbulence  model  for  the  flow  through  a  wing-body  junction.  Their 
primary  interest  was  to  investigate  the  effects  of  adding  triangular  shaped  fairings  to  the 
leading  and  trailing  edges  of  the  appendage,  similiar  to  fairings  used  by  Shin  (1989). 
From  the  three  different  leading  edge  fairing  configurations,  they  found  that  a  longer 
fairing,  which  extended  further  upstream,  w'as  optimum  in  terms  of  reducing  the  mean 
velocity  distortion  in  the  wake.  They  also  determined  that  trailing  edge  fairings  tu'e 
marginally  beneficial  when  combined  with  a  leading  edge  fairing.  Maximum  total  drag 
reduction  was  5%  less  than  the  no  fairing  case. 


1.2  Author’s  Research  Program  and  Objective.' 

1.2,1  Experimental  Research  Goals 

The  experimental  research  program  of  the  author  was  designed  to  meet  several 
goals.  1)  To  furtlicr  expand  the  3:2  elliptic  nose,  0020  tail  junction  flow  database  with 
high  quality  fiow  uncertainty)  measurements.  2)  To  examine  the  data  closely  to  gain  a 
more  complete  understanding  of  the  physical  nature  of  the  flow  3)  Mo.st  iinpon:.ntly,  to 
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compare  these  data  to  other  measurements  taken  in  juncture  flows  with  the  same 
geometry.  The  other  data  sets  to  be  used  for  comparisons  are  shown  in  Table  1.1.  Table 
3.3  shows  the  inflow  conditions  for  each  data  set,  and  section  3.2  discusses  the  varying 
approach  boundary  layer  parameters.  Through  comparisons,  the  effects  of  Reg*,  Re^, 
5‘/T,  U„/r,  Ur^6,  and  other  parameters  on  the  junction  flow  will  be  more  fully 
understood. 

A  final  goal  is  to  compare  the  author’s  ’.tot-wire  data  to  that  obtained  by  LDV 
where  possible.  Olcmen  (1990)  noticed  discrepancies  in  comparing  the  results  of  hot-wire 
and  LDV  measurements.  He  found  that  near  wall  turbulence  measurements  and  mean 
velocity  measurements  very  near  the  wall  (Y’’  <  50)  did  not  agree  between  the  two 
experimental  techniques.  Additional  information  is  sought  concerning  this  problem.  Cf 
estimates  were  also  made  and  compared  to  Ailinger’s  (1990)  results. 

1,2  J  Description  of  Numerical  Work 

In  addition  to  wind  tunnel  measurements,  numerical  computations  using  vortex 
panels  were  used  to  gain  insight  into  die  effects  of  geometiy  on  the  horseshoe  vortex 
structure.  This  was  accomplished  by  using  a  2-D  inviscid  vortex  panel  calculation  given 
by  Kuethe  and  Chow  (1986)  to  predict  the  velocity  distributions  around  various 
appendage  shapes. 

The  geometry  of  an  appendage  effects  the  distribution  and  total  amount  of 
vorticity  (circulation)  present  in  a  horseshoe  vortex  in  two  ways.  High  velocity  gradients 
around  the  appendage  nose  lead  to  high  rates  of  vortex  stretching,  increasing  the  peak 
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vorticity  and  affecting  the  shape  of  the  vonex  structure.  In  addition,  the  geometrical 
bluntness  of  the  appendage  causes  a  stronger  adverse  pressure  gradient,  increasing  the 
initial  vortex  roll-up,  which  increases  the  circulation  in  tlte  horseshoe  vortex. 

Mehta  (1984)  experimentally  studied  the  horseshoe  vorticies  produced  by  different 
appendage  nose  shapes.  He  verified  that  increased  appendage  bluntness  causes  a  stronger, 
more  structured  vortex.  Also,  by  comparing  their  work  to  other  researchers,  Kubendran 
et  al.  (1986)  found  that  the  horseshoe  vortex  path  and  its  strength  aie  related  to  the  nose 
slenderness  ratio.  They  report  that  the  secondary  flow  in  a  junction  created  by  an 
appendage  with  a  large  nose  slenderness  ratio  (12:1)  is  mainly  due  to  the  cross-stream 
gradients  of  the  Reynolds  stresses  (Prandtl’s  "second  kind"  secondary  flow). 

The  numerical  work  was  carried  out  by  computing  the  velocity  distributions 
around  a  wide  variety  of  appendage  shapes.  The  final  results  include  correlations  between 
appendage  geometries  and  average  vortex  stretching  rates,  and  between  a  vorticity 
transport  parameter  (determined  numerically)  and  a  non-dimensional  circulation  estimate 
from  Mehta’s  experimental  work  (1984).  The  numerical  research  goals  and  results  are 
fully  discussed  in  chapter  6. 
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2.0  Experimental  Apparatus  and  Techniques 


Unless  otherwise  noted,  all  data  presented  were  taken  in  the  Virginia  Tech  Low 
Speed  Boundaj."y  Layer  Tunnel.  This  tunnel  has  been  used  in  much  previous  work  at  VPI 
and  SU  as  well  as  Southern  Methodist  University  (Ahn,  1986,  Devenport  and  Simpson, 
1990a,  1990b,  Devenport  et  ai,  1989,  Simpson  et  al.,  1980). 

This  chapter  describes  the  appendage,  wind  tunnel,  instrumentation  and 
experimental  techniques,  data  collection  and  reduction,  and  uncertainty  estimates. 


2.1  Appendage  Description 

The  appendage  (or  wing)  used  in  this  study  was  made  up  of  a  3:2  elliptic  nose 
joined  to  a  NACA  0020  tail  section  at  its  maximum  thickness.  The  section  coordinates 
are  given  in  Table  2.1,  and  Figure  2.1  shows  its  cross  section.  As  mentioned  in  chapter 
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1,  this  body  shape  has  been  studied  in  many  other  experiments.  It  has  a  maximum 
thickness  (T)  of  7.17  cm,  a  chord  (C)  of  .W.5  cm,  and  is  22.9  cm  high. 

To  prevent  unsteadiness  due  to  natural  flow  transition,  boundary  layer  trips  made 
from  6.35  mm  wide  strips  of  120  grade  sandpaper  were  attached  to  the  appendage  with 
the  sandpaper  leading  edge  positioned  28.2  mm  downstream  of  the  wing  leading  edge. 
Appendage  boundary  layer  data  were  presented  by  Devenport  and  Simpson  (1990a)  using 
these  trips.  The  appendage  boundary  layer  thickness  at  its  trailing  edge  is  =  .25  T  in  the 
two  dimensional  flow  region.  The  symmetry  of  the  appendage  boundary  layer  was  taken 
to  be  a  good  indication  of  overall  flow  symmetry  (Figure  3.7).  Devenport  and  Simpson 
also  report  no  separation  is  visible  in  the  two  dimensional  flow  region  of  the  appendage 
(1.06  <  Y/T  <  2.12)  using  oil-flow  visualization.  Section  3.1,2  presents  the  data 
describing  the  appendage  boundary  layer. 


2.2  Wind  Tunnel 

ITie  Virginia  Tech  AOE  Low  Speed  Boundary  Layer  Tunnel  test  section  is  shown 
in  Figure  2.2.  The  test  section,  which  measures  8  m  long  by  .91  m  wide,  can  be  divided 
into  3  regions  as  shown.  These  different  regions  will  be  described  in  following  sections. 
The  test  flow  for  this  open  circuit  tunnel  is  driven  by  a  25  horsepower  centrifugal  blower. 
The  flow  passes  through  an  air  filter,  blower,  fixed- .setting  flow  damper  (to  control  flow 
speed),  a  section  of  honeycomb  to  sr.-aighten  and  remove  flow  swirl,  7  turbulence  screens, 
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and  a  two  dimensional  4:1  contraction  nozzle  before  entering  the  test  section.  The  test 
section  floor  is  constructed  of  fin-forrn  plywood  19  mm  thick,  except  for  two  access 
panels  around  the  wing,  which  are  6.4  mm  plexiglas.  Tne  side  walls  are  float  plate  glass, 
and  the  ceiling  is  made  of  6.4  mm  plexiglas  reinforced  with  aluminum  channel. 

2.2.1  Test  Section  Description,  Regions  1  and  2 

Within  the  first  1.63  m  of  region  1,  the  flow  undergoes  a  further  1.5:1  contraction, 
ending  at  a  throat  25.4  cm  in  height.  From  the  throat  to  the  end  of  region  2,  the  upper 
wall  diverges  slowly  to  a  height  of  26.67  cm  over  a  length  of  3.54  m.  This  divergence 
reduces  the  variation  in  Cp  of  the  test  flow  due  to  boundary  layer  growth  ttlong  the  floor 
and  ceiling.  A  6.3  nun  forward  facing  step  trips  the  boundary  layer  at  the  beginning  of 
region  1. 

Plexiglas  wall  liners  6.4  mm  thick  support  the  ceiling  in  regions  1  and  2. 
Positioning  the  wall  liners  330  mm  (X/C  =  1.08)  upstream  of  the  appendage  leading  edge, 
and  203  mm  (X/C  =  .605)  downstream  of  the  trailing  edge  allowed  the  appendage 
blockage  effects  to  be  'duced  by  crudely  approximating  its  two  dimensional  potential 
flow  streamlines  (Figure  2.3).  Note  that  the  discontinuities  caused  by  the  wall  liners  were 
faired  over  vdth  tape.  By  comparing  potential  flow  results  to  measurements  of  the  static 
pressure  at  the  wall,  Devenport  and  Simpson  (1990a)  show  the  absence  of  any  substantial 
blockage  effects  in  the  flow. 
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2.2  J,  Test  Section  Description^  Region  3 

Re,^on  3  of  the  test  section  was  modified  by  the  author  for  gathering  appendage 
far-wake  data.  It  was  designed  to  study  the  decay  of  the  trailing  vortex/wake  structure 
in  the  presence  of  an  adverse  pressure  gradient.  This  region  utilizes  the  tunnel’s  active 
boundary  layer  control  to  bleed  off  a  portion  of  the  incoming  side  wall  boundary  layers. 
Figure  2.4  shows  oil  flow  results  of  the  areas  around  the  boundary  layer  control  suction 
slots.  The  oil  flows  and  tufts  were  used  to  adjust  the  mass  flow  through  tlie  slots  so  that 
the  inner  side  wall  boundary  layer  would  be  bled  off  while  maintaining  a  two  dimensional 
flow  in  the  freestream.  The  active  boundary  layer  control  system  is  further  described  by 
Simpson  et  al.  (1980).  A  boundary  layer  "scoop"  was  built  into  the  ceiling  of  region  3 
(Figure  2.5).  Its  piuposc  was  to  skim  off  the  inner  2  cm  of  the  ceiling  boundary  layer, 
and  thus  prevent  separation  in  the  diverging  section. 

Figure  2.5  also  presents  the  dimensions  of  the  diverging  section.  Research  by 
Reneau,  Johnston,  and  Kline  (1967),  shows  that  for  a  diffuser  with  these  dimensions,  no 
separation  should  occur  downstream  of  the  throat.  The  inner  side  wall  liners  were  made 
of  9.5  mni  ihick  piexigias,  while  the  ceiling  sections  were  made  of  6.4  mm  thick  plexiglas 
reinforced  with  alummum  channel.  The  scoop  section  was  constructed  of  plexiglas  3.2 
ram  thick  to  facilitate  the  needed  ceiling  curvature.  The  leading  edge  was  tapered  to  1.6 
mm,  rounded,  and  angled  upwards  approximately  2°  to  prevent  separation  off  the  lower 
side  of  the  scoop.  Figure  2.6  shows  a  portion  of  the  results  from  an  oil-flow  visualization 
performed  at  the  scoop  leading  edge.  No  separation  was  in  evidence  at  the  scoop  leading 
edge. 

Experimental  .Apparatus  and  Techniques  21 


A  back-pressure  screen  was  constructed  to  match  the  static  pressure  in  region  2 
to  the  atmospheric  pressure  and  minimize  air  leakage  into  the  test  section.  After  adjusting 
the  open  area  ratio  of  the  screen  with  tape,  the  final  pressure  difference  was 
approximately  -.1  inches  of  water  (atmospheric  pressure  slightly  greater  than  region  2 
pressure).  Tliis  allowt-d  for  a  small  factor  of  safety,  as  a  positive  pressure  difference 
would  tend  to  raise  tlie  v;st  section  ceiling  in  regions  1  and  2.  The  ceiling  section  in 
region  3  had  to  be  helo  .r  place  due  to  the  increasing  internal  pressure  as  the  flow 
diverged.  Oil  flow  vis’ioliztttion  was  used  extensively  in  region  3  to  assure  that  the  flow 
was  not  separating;;,  p;u-iicularly  in  the  side  wall/floor  junction  areas.  The  oil  flows 
showed  no  substantial  separation  present  for  X/C  <  12,  where  C  is  the  appendage  chord 
length. 

2.2.3  Appendage  Mounting 

The  leading  edge  of  the  appendage  was  located  1.39  m  downstream  of  the  throat 
in  region  2..  The  ceiling  curvature  defining  the  start  of  region  3,  begins  2.24  m 
downstream  of  the  appendage  leading  edge.  The  body  was  mounted  at  zero  incidence  and 
sweep  along  the  tunnel  centerline.  A  37  mm  gap  was  left  between  the  end  of  the 
appendage  and  the  tunnel  ceiling  to  prevent  the  formation  of  a  second  junction  vortex 
along  the  ceiling. 
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2.2.4  Wind  Tunnel  Performance 

The  flow  characieristics  of  regions  1  and  2  have  been  documented  in  previous 
work.  Devcnport  and  Simpson  (1990a)  report  that  the  freestream  flow  of  region  2  is 
uniform  to  within  .5%  in  the  spanwisc  (Z)  direction  and  1%  in  the  vertical  (Y)  direction, 
with  a  rms  turbulence  intensity  of  .2%  at  =  27  m/s.  Detailed  documentation  of 
regions  1  and  2  was  performed  by  Ahn  (1986).  He  reports  that  the  boundary  layer  in 
regions  1  and  2  is  an  cquilibiium  boundary  layer  and  that  it  closely  satisfies  the  two 
dimensional  momentum  integral  equation.  Also,  the  velocity  and  pressure  spectra  show 
no  preferred  frequencies.  Figure  2.7  shows  mean  velocity  measurements  in  the  potential 
core  of  regions  1  and  2,  and  displays  the  lack  of  any  substantial  pressure  gradients  past 

V  : - 1.-_  /I  >~rrt  \ 
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The  flow  characteristics  in  the  diverging  section  (region  3)  were  documented  by 
the  author  and  will  be  discussed  in  section  3.1. 


2.3  Experimeiitai  Techniques  and  Instrumentation 

2.3.1  Coordinate  Systems 

Four  different  coordinate  systems  were  used  to  reduce  and  present  the  data;  wire, 
local  stream,  local  freestream,  and  tunnel  coordinate  systems.  The  Y-axis  for  all  of  these 
systems  is  always  perpendicular  to  the  floor.  Figure  2.8  shows  the  relationship  between 
tunnel  and  local  freestream  coordinates.  The  local  freeso'eam  coordinates  are  defined  by 
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setting  VV  =  0  in  the  freestream.  This  coordinate  system  is  fixed  with  respect  to  Y,  but 
iti  orientaL'on  changes  as  X  and  Z  change.  The  tunnel  coordinate  system  is  fixed  for  all 

X,  Y,  and  Z,  and  its  orientation  is  shown  in  Figure  1,1.  The  angle  a  is  defined  positive 
by  right-hand  rule  convention,  and  when  subscripted,  the  subscript  refers  to  the  coordinate 
system  from  which  the  angle  is  measured  (the  angle  in  Figure  2.8  would  be  otj-e  )•  'i’he 
majority  of  the  data  is  presented  in  tunnel  coordinates. 

Figiii'e  2.9  shows  two  coordinate  systems  fundamental  to  the  initial  data  reduction, 
the  wire  and  local  stream  coordinate  systems.  The  wire  system  is  fixed  with  respect  to 
the  hot-wire  probe.  The  local  stream  coordinate  system  is  defined  by  W  =  0  at  every  X, 

Y,  and  Z  location.  Outside  the  boundary  layer,  the  local  freestream  and  local  stream 
coordinates  should  be  identical,  figure  2.9  also  shows  the  local  flow  angle  being 
measured  from  some  arbitrary  coordinate  system.  The  importance  of  the  wire  and  local 
stream  coordinate  .systems  can  be  seen  m  section  2.5.2. 

2.12  Instrumentation  anC  Data  Collection 

All  data  collection,  with  the  exception  of  the  streamwise  traverses  made  in  the 
potential  core  of  region  3,  was  performe<i  with  a  TSI  model  1218-T1.5  boundary  layer 
hot-wire  probe.  The  strci'mwise  traverse  data  was  collected  with  a  TSI  model  1210'T1.5 
straight  hot-wire  probe.  Fach  pmbe  sensor  is  approximately  ,004  mm  in  diameter  and  1.3 
mrn  in  lengdi. 

For  both  probes,  a  Miller-t  .'pe  integrated  circuit  constant  temperature  hot-wire 
anemometer  was  used.  Details  of  this  anemometer  iire  given  by  Miller  (1976).  A  method 
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of  frequency  response  testing  using  a  sine  wave  generator  is  given  by  Simpson  et  al. 
(1979).  Employing  this  method,  the  author  found  the  frequency  response  of  the 
anemometer  to  be  in  excess  of  20  kHz  for  an  overheat  ratio  of  1.7  and  U^ff  =  27  m/'s. 

Calibrations  were  made  witii  a  TSI  model  1 125  calibrator  from  flow  speeds  of 
approximately  9.5  m/s  to  29.5  m/s.  Using  King’s  Law  for  a  constant  temperature 
anemometer  in  incompressible  flow, 

=  A*  +  (2.1) 

the  velocity  versus  output  voltage  relationship  could  be  linearized  within  the  data 
reduction  computer  program.  The  piobe  was  calibrated  by  fixing  the  exponent  n  at  .45 
and  then  finding  the  coefficients  A*  and  B‘  using  a  linear  least-squates  curve  fit.  Tlie 
computer  then  varied  n  to  obtain  the  maximum  correlation  coefficient.  The  final  value 
of  n  was  usually  in  the  range  .39  <  n  <  .41.  Correlation  coefficients  in  excess  of  .99999 
were  attainable  using  this  method. 

Figure  2.10  shows  a  block  diagram  of  the  equipment  used  for  calibration  and  data 
acquisition.  An  EBM  AT  computer  with  a  12  bit  Data  Translation  DT2801  A/D  converter 
was  used  to  sample,  reduce,  and  record  the  anemometer  signals.  The  offset  amplifier, 
with  a  cutoff  frequency  of  approximately  40  kHz,  served  as  a  signal  conditioner  for  the 
A/D  board  and  also  as  a  signal  buffer.  The  voltmeters  were  used  to  calibrate  the 
amplifier  (which  was  done  every  60  to  90  minutes)  and  also,  along  with  the  oscilloscope, 
to  monitor  the  anemometer  signal. 

Us,  m's,  w's,  -tnvs,  and  cqc  measurements  were  taken  at  a  single  location  by 
yawing  the  hot-wire  (held  in  the  XZ  plane)  about  the  Y-axis  (perpendicular  to  the  floor) 
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to  three  different  angles.  These  angles  were  approximately  ±30°  and  0°  relative  to  the 
local  streamwise  direction.  This  procedure  is  explained  in  greater  detail  in  section  2,5. 
At  each  orientation,  7000  samples  were  taken  at  5  kHz,  for  a  record  length  of  1.4 
seconds.  The  effective  record  length  was  longer  than  this,  since  the  computer  processed 
blocks  of  1000  samples,  each  block  taking  approximately  5  to  8  seconds  to  process. 
Therefore,  the  7000  samples  were  distributed  over  35  to  50  seconds.  No  time  series  data 
were  recorded,  as  all  the  initial  data  processing  was  done  by  the  data  acquisition  program. 

To  take  measurements,  the  probe  was  positioned  in  the  flow  on  a  rotating  support 
(Figure  2.11).  Because  the  probe  sensor  is  offset  from  its  support  axis,  a  "dog-leg  ”  probe 
holder  was  used  to  pmsition  the  sensor  at  the  center  of  rotation.  The  probe  support  was 
held  by  a  lead  screw  traverse  system  positioned  above  the  test  section  ceiling,  enabling 
movement  in  the  YZ  plane.  The  initial  probe  position  from  the  wall  was  determined 
using  a  cathetometer  to  measure  the  distance  from  the  probe  to  its  reflected  image  from 
the  floor.  Table  2.2  presents  the  estimated  X,  Y,  and  Z  location  uncertainties.  Note  that 
the  uncertainties  in  X  and  Z  increase  with  height  above  the  floor  due  to  slight  probe 
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2.4  Uncertainty  Estimates 

One  of  the  goals  for  these  measurements  to  obtain  low  uncertainty  data.  The 
"jitter"  program  method,  presented  by  Moffat  (1982),  was  used  to  estimate  the  95% 
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confidence  levels  for  the  reduced  data  using  the  data  reduction  technique  discussed  in  the 
next  section.  The  experimental  uncertainties  are  shewn  in  Table  2.3.  Nominal 
uncertainty  estimates  for  the  reduced  data  (tunnel  or  freestream  coordinates)  are  shown 
in  Table  2.4,  Although  die  jitter  program  estimates  a  5%  uncertainty  for  the  w'  and  -im 
data,  the  scatter  in  the  results  seem  to  indicate  that  7%  to  10%  may  be  a  better,  or  even 
optimistic,  uncertainty  estimate. 

To  aid  m  future  measurements,  the  contributions  of  each  of  the  experimental 
uncertainties  to  the  total  uncertainty  of  the  reduced  results  were  broken  down  by 
percentage.  The  major  contributors  for  each  quantity  can  be  seen  in  Table  2.5.  Of  note 
is  the  major  contribution  (greater  than  95%)  of  probe  yaw  angle  uncertainty  to  the  total 
uncertainty  in  otre  and  W^c. 

It  should  be  noted  that  these  uncertainty  estimates  do  not  account  for  systematic 
or  bias  errors  in  the  experiment.  These  figures  merely  represent  the  likely  scatter  due  to 
random  effects.  Systematic  errors  were  hard  to  avoid  when  aligning  the  probe  with  the 
tunnel  coordinate  system.  Using  a  2-D  inviscid  flow  calculation,  attempts  were  made  to 
identify  and  correct  for  these  errors.  Section  4.3.3  discusses  the  Gorrections  made  for 
these  alignment  errors.  Another  known  source  of  systematic  eiTors  can  be  seen  in  some 
of  the  profile  plots  of  m  jc'^ref  ^rtd  (plane  10  data  in  Figures  4.99  and  4.101,  for 

example),  A  small  amount  of  noise  (»  3  mV  rms)  could  be  seen  using  an  oscilloscope 
when  the  hot-wire  probe  was  in  the  freestream.  The  origin  of  the  noise  was  finally  traced 
to  a  computer  monitor,  and  moving  the  monitor  solved  the  problem,  The  other 
fluctuating  quantity  {-uw)  was  not  affected  due  to  the  nature  of  the  data  reduction 
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equations  and  only  the  freestream  data  seemed  to  be  affected  for  u'  and  w'.  To  check 
repeatibility,  a  test  was  performed  by  taking  multiple  and  ou-c  measurements  at  X/C 
=  1.05,  Z/T  =  -.2.  The  test  showed  the  mean  streamwise  velocity  results  to  be  repeatible 
within  1%  of  U„f,  with  nominal  differences  in  o^q  of  1.5°. 

Finally,  the  errors  of  the  other  data  sets  examined  in  this  study  need  to  be 
estimated  to  compare  with  the  author’s  data.  Table  2.6  lists  the  errors  given  in  the 
references  for  each  of  the  data  sets.  McMahon  et  al.  and  Dickinson  both  give 
uncertainties  in  the  same  form  as  the  author,  as  a  percentage  of  measured  value. 
Devenport  and  Simpson  report  uncertainties  as  absolute  quantities,  regardless  of  the  value 
being  measured.  The  errors  listed  for  McMahon’s  case  are  the  maximum  possible  errors, 
and  90%  of  the  uncertainties  are  only  half  of  the  given  maximum  values.  Dickinson 
reports  uncertainties  only  for  mean  quantities.  Not  listed  in  Table  2.6  are  Ailinger’s  skin 
friction  uncertainties.  For  the  locations  in  planes  5,  8,  and  10,  the  uncertainties  of  his 
data  ranged  from  2%  to  10%.  A  good  nominal  uncertainty  estimate  would  be  5%  or  6% 
of  the  measured  Q  value. 

The  experiment dl  techniques  Vcirieci  for  each  dma  set,  3.nd  this  in  itself  ma.y  have 
brought  about  bias  errors,  or  "bias  differences"  between  each  data  set.  A  discussion  on 
this  subject  is  beyond  the  scope  of  this  thesis,  but  die  reader  should  be  aware  of  this  fact, 
it  is  the  author’s  opinion  that  none  of  the  data  sets,  as  a  whole,  suffer  from  uncenainties 
or  bias  errors  that  would  invalidate  general  comparisons. 
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2,5  Data  Reduction 


A  method  previously  used  by  Devenport  and  Simpson  (19?0a)  was  used  to  obtain 
the  desired  quantities.  This  method  avoids  the  problem  of  having  to  linearize  every 
voltage  sample  and  thus  speeds  up  tiie  data  reduction. 


2.5.1  Calculation  of  U,^  and  u',^ 
By  letting 


f(E'^  can  be  expanded  in  a  Taylor  series  about  E.  In  doing  so.  U..,  and  u  „  n.m  be.  found 
in  terms  of  E,  E^,  E^,  and  E'*  (note  that  these  are  mean  quantities,  e.g.  the  average  of  the 
square  of  the  voltages).  These  quantities  can  be  obtained  by  looping  through  the  voltage 
arrays  only  one  time.  A  more  complete  listing  and  derivation  of  the  equations  involved 
is  given  in  appendix  A. 

To  test  the  accuracy  of  this  method,  a  simulation  was  done  using  a  computer 
routine  that  generated  a  set  of  random  samples  with  a  Gaussian  distribution.  Since  the 
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mean  and  variance  of  this  simulated  velocity  distribution  were  known  a  priori,  the  errors 
using  this  method  could  be  calculated,  given  nominal  values  of  the  calibration  coefficients 
A*,  B*,  and  n.  A  record  length  of  5,000  data  samples  were  used,  and  the  coefficients 
were  A*  =  .1225,  B*  =  .09,  and  n  =  .41.  No  attempt  was  made  to  simulate  the  skewness 
present  in  a  real  boundary  layer  velocity  histogram. 

Figures  2. 1 2a  and  2. 12b  show  the  results  of  the  simulations  using  various  Uj^cajAiref 
values.  The  calculated  errors  diverged  rapidly  when  m'AJioc»i  became  greater  than 
approximately  .25.  The  absolute  values  of  the  errors  are  plotted,  so  it  should  be  noted 
that  a'AJref  always  underestimated  by  tlte  method  in  Appendix  A.  As  seen  from 
Figure  2,12b,  both  of  the  error  functions  seem  to  have  a  correlation  in  the  form  log(eiTor) 
=  log('f/Ujoj^)  +  constant.  Below  the  divergence  limits,  the  absolute  value  for  U^rropAJiocji 


was  less  than  .0005  and  was  less  than  .003.  The  maximum  local  turbulence 

value  ever  encountered  is  usually  around  .30  for  an  equilibrium  boundary  layer,  so  this 
method  is  quite  accurate  at  common  m'AJiooii  values.  The  maximum  local  turbulence 
levels  encountered  by  the  author  in  the  vicinity  of  the  wing/body  junction  was  .27,  in 

1 0  of  fK/»  »i/ir*rr  M~oilinf»  fr^ovirniino  T»oln#>p  ol-d/nx/c  n<»oT*  fVi<»  xi/oll 

^  AVU  AW  /Vf  AAAW  Uk  i.  A  WK^^W  »  A  A  AW  i  A  AkA./ X  A  A  A  ACA  A  i  A  T  VAW  >T  WA  W  LAA  T  T  LA  ^  1  A  WAAA  AilW  TT  AAAi.^ 

usually  at  Y/T  <  .01.  Nominal  near  wall  values  of  in  plane  10  were  .l4  to  .25. 

In  other  planes  near  the  appendage,  the  nominal  near  wall  values  were  rarely  greater  than 
.20.  In  region  3,  the  near  wall  values  of  were  much  higher,  with  nominal  values 

of  .2  to  .31,  and  a  maximum  value  of  .39  in  plane  15.  For  this  reason,  attempts  were 
made  to  conect  the  data  in  region  3. 


:‘0 
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The  corrections  were  made  by  curve  fitting  the  error  functions  (Figure  2. 12b).  The 
final  equations  were 


UenAJlocU  =  (e-"''*'")(n'/U,„^ 

PxY  =  .9973 

(2.5a) 

PxY  =  -9984 

(2.5b) 

The  correction  criteria  used  was  a  correction  greater  than  1%  of  the  value  being  conected, 
and  a  local  tiurbulencc  intensity  of  more  than  .05.  Using  these  criteria,  no  U/U„f  data  was 
corrected.  In  plane  13,  31  corrections  were  made,  with  an  average  correction  of  -.001. 
In  plane  14,  96  corrections  were  made,  with  an  average  of  --.0012,  and  in  plane  15,  180 
corrections  were  made,  with  an  average  of  -.001 1.  The  maximum  correction  ever  made 
was  only  2.5%  of  the  original  value,  with  the  nominal  collections  in  the  range  1%  to 
1.5%.  The  effects  of  tliese  corrections  can  be  seen  in  die  profile  plots  of  the  u'  data  in 
planes  13-15  (see  section  4.3.2). 

2.5.2  Calculation  of  Ug,  u's,  h*  5,  and  ctj 

As  mentioned  earlier,  the  procedure  for  obtaining  Us,  u's,  w'g,  -wvvs,  and  ttj  in 
planes  5  through  11  is  the  same  as  descn’oed  by  Devenport  and  Simpson  (1990a).  This 
procedure  is  similiar  to  the  "conventional"  method  described  by  Rodi  (1975).  He  reports 
that  this  method  is  accurate  for  turbulence  intensities  below  25%  of  U,(,f.  The  two 
coordinate  systems  used  for  the  calculations  are  shown  in  Figure  2.8.  The  angle  a  is 
positive  as  shown  followmg  right  hand  rule  convention.  Using  the  equation 
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’’r. 

■'i 

c/^  =  ul  ^ 

(2.6) 

• 

as  shown  in  Jorgensen  (1971),  the  expression  for  the  effective  cooling  velocity  can  be 

written  as 

• 

C/Jj  =  (J7^a  +  Pr^cE)^  + 

(2.7) 

+  (-f/jsina  +  PK^osa)^ 

where  a  has  been  measured  with  respect  to  the  wire  coordinate  system. 

One  can  expand 

• 

equation  (2.7)  into  a  more  convenient  form: 

■■ : 

c4  =  AU]  +  iJFj  +  CPP^  +  FU^W^ 

(2.8) 

• 

where 

.  -’r 

A  =  cos^a  +  /t^sin^a 

B  =  h} 

• 

C  =  sin^a  •*  it^cos^a 

F  =  2(l-it^)cosa  sina 

Since  Ug,  Vg,  and  Wg  are  components  in  the  time-mean  flow  coordinate  system. 

# 

Wg  =  0  by  definition,  and 

'v'*? 

(2.9) 

'  •< 

II 

•V 

If  it  is  further  assumed  that 

'  ■  '"j" 

Mg  ,  ,  Wg  ,  Fj  <  .25f/g 

then 

• 

•t' 
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A 

Ug  ,  Vj  ,  Wg  ,  UWg  ,  KjVj  » 

After  dropping  these  second  order  terms,  the  relevant  equations  become 

(x'a,)^  =  ^  “""j 

4  vl 

where  A  and  F  arc  functions  of  the  local  stream  mean  flow  angle  measured  with  respect 
to  the  probe,  which  changes  as  the  probe  is  yawed.  Using  equation  2.10  and  two 
measurements  taken  at  yaw  angles  of  approximately  ±30°,  the  flow  yaw  angle  with 
respect  to  an  arbiaary  datum  and  the  magnitude  of  the  mean  velocity  vector  can  be  found. 
The  angle  damm  was  selected  such  that  a  =  0°  when  the  velocity  vector  v/as  aligned  with 
Xjc-  After  taking  another  measurement  at  a  third  yaw  angle,  preferably  aligned  with 
a  simultaneous  system  of  equations  can  be  solved  to  yield  values  for  u's,  w'j,  and  -uw^, 
witli  redundant  values  of  Uj  and  ots-  The  coefficient  of  axial  sensivity,  k,  was  taken  to 
be  .13  for  these  measurements.  There  is  a  great  amount  of  uncertainty  in  this  figure,  but 
as  Table  2.5  shows,  this  uncenainty  does  not  greatly  effect  the  reduced  quantities. 

The  output  from  the  data  acquistion  program  (mean  velocity  vector  magnitude  and 
direction,  and  3  components  of  the  Reynolds  stress  tensor  in  the  XZ  plane)  was  given  in 
the  local  stream  coordinate  system.  This  output  was  then  transformed  into  the  tunnel 
coordinate  system.  Since  the  initial  probe  rlignment  was  done  by  eye,  a  relatively  large 
error  could  propagate  through  the  results  of  data  presented  in  tunnel  coordinates.  Profile 
data  presented  in  the  local  frecstream  coordinate  system  would  eliminate  the  alignment 
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error,  since  the  relative  angular  measurements  for  different  Y  locations  in  each  profile  are 
more  accurate.  The  tunnel  coordinate  system  was  used,  however,  since  it  was  considered 
better  for  comparing  data  sets.  Note  that  for  planes  13,  14,  and  15,  the  probe  was  not 
rotated  and  only  Ujc  and  u\q  were  measured.  This  was  done  mainly  in  the  interest  of 
time,  as  this  meant  only  a  third  as  much  time  was  spent  taking  data.  Also,  it  is 
reasonable  to  assume  that  much  of  the  secondary  flow  structure  present  from  the 
horseshoe  vortex  becomes  diffused  and  washed  out,  particularly  as  the  flow  diverges  in 
the  adverse  pressure  gradient. 
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3.0  Description  of  Test  Conditions 

An  effort  was  made  to  keep  the  daily  test  conditions  as  uniform  as  possible.  The 
wind  tunnel  lab  area  was  air-conditioned,  and  temperature  equilibrium  was  established  at 
approximately  25°  C  after  running  the  wind  tunnel  and  laboratory  air-conditioner  for  30 
minutes.  Constant  temperature  was  important  for  accurate  hot-wire  measurements  and 
also  to  keep  flow  characteristics  (flow  properties,  separation  locations,  etc.)  from  varying. 
The  daily  variation  in  temperature  was  estimated  at  ±  .5“  C.  Note  that  a  temperature 
uncertainty  of  ±  1°  C  was  used  in  table  2.3.  This  figure  was  used  as  a  worst  case 
estimate. 

The  daily  atmospheric  pressure  was  approximately  945  ±10  mb,  and  the  daily 
reference  velocity  was  approximately  26.75  ±  .75  m/s.  Small  changes  in  were 
constantly  updated  as  measurements  were  made. 
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3.1  Flow  Characteristics 


3.1,1  Empty  Tunnel  Flow  Characteristics 

To  document  the  flow  characteristics  in  the  adverse  pressure  gradient  region, 
potential  core  and  boundary  layer  measurements  were  taken  in  regions  2  and  3  without 
the  wing.  When  taking  these  measurements,  side  wall  inserts  were  used  in  region  2  to 
fill  die  gaps  used  to  approximate  the  potential  flow  streamlines  when  the  wing  was 
installed. 

Data  from  the  potential  core  streamwise  traverses  are  shown  in  Figure  3.1.  A  1- 
dimensional  inviscid  continuity  calculation  was  done  for  comparison  to  the  UAJref  and  Cp 
results,  and  is  shown  in  Figure  3.1  by  the  solid  line.  The  variation  is  mainly  due  to  the 
additional  blockage  caused  by  the  growing  boundary  layer.  An  increase  in  the  freestream 
value  of  u'  from  .2%  to  .7%  occurs  as  the  flow  passes  through  region  3.  Note  the  slight 
"hump"  in  the  UAJfef  and  Cp  data  where  region  3  begins.  This  is  due  to  the  blockage 
effects  of  the  suction  slots  and  the  boundary  layer  bleed  scoop  at  the  ceiling.  The  change 
in  Cp  with  respect  to  X/C  was  estimated  using  central  differencing.  The  results  for 
d(Cp)/dQC/C)  are  also  shown  in  Figure  3.1. 

Figure  3.2  shows  U/U„f  and  u'/U^f  profile  data  for  3  locations  in  region  3.  The 
mean  velocity  seems  to  decrease  almost  uniformly  as  X/C  increases.  The  turbulence 
intensity  in  the  inner  boundary  layer  region  decreases  as  Ups  decreases,  contrary  to  the 
turbulence  behavior  in  the  freestream.  As  seen  from  Figure  3.2,  the  boundary  layer 
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profiles  seem  to  have  similiar  shapes.  Calculations  were  done  to  determine  if  these 
profiles  resulted  from  an  equilibrium  boundary  layer  which  could  be  described  by 

« (X-XX 

as  given  by  Rotta  (1962).  These  profiles,  however,  failed  to  meet  the  criteria  required  for 
an  equilibrium  turbulent  boundary  layer. 

Using  Spalding’s  wall  law  (Spalding,  1961)  for  tne  inner  region,  Q  estimates  were 
calculated  for  the  three  boundary  layer  profiles  measured  in  region  3.  Spalding’s  equation 
is  given  as 

Y*  =  )  (3.2) 

G  21  31  41 

where  k  =  0.407  and  G  =  10.  Figure  3.3  gives  an  example  of  how  well  equation  (3.2) 
fit  these  profiles.  Table  3.1  lists  the  profile  parameters.  The  displacement  thickness  (5*) 
results  were  used  to  perform  a  more  accurate  1-D  '  inviscid"  calculation  for  the  velocity 
distribution  in  region  3.  The  dashed  line  in  Figure  3.1  shows  the  results  of  the 
calculation.  Tbe  calculated  U/U,^f  values  with  the  5*  correction  are  closer  to  the 
measured  values,  but  as  X/C  increases,  the  results  diverge.  To  perform  the  calculation 
coiTccdy,  the  boundary  layer  growtli  along  the  tunnel  roof  and  sidewalls  must  also  be 
accounted  for.  A  simple  calculation  shows  that  although  the  physical  angle  of  divergence 
of  region  3  is  4.996°,  the  vs.  X/C  potential  core  data  show  a  mean  "inviscid” 

divergence  angle  of  3.276°  with  a  standaid  deviation  of  1.16%  for  X/C  >  8. 
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3.12  Appendage  Flow  Characteristics 

As  previously  mentioned,  a  large  amount  of  data  are  presented  regarding  this  flow 
in  Devenport  and  Simpson  (1990a).  Some  of  the  data  att  pre.sented  here  to  give 
additional  information  about  the  wing-body  junction  flow  chaj  .v'Uvristics. 

Figure  3.4  shows  a  surface  plot  of  the  Cp  distribution  cn  the  .cst  wail  surrounding 
the  appendage,  and  is  good  for  showing  the  qualitative  features  of  the  pressure  field. 
Note  the  large  adverse  pressure  gradient  imposed  on  the  incoming  boundary  layer  at  the 
nose  of  the  appendage.  As  the  flow  is  accelerated  around  the  appendage  sides,  the 
pressure  reaches  a  minimum  value  just  forward  of  the  maximum  thickness.  At  the 
minimum  pressure  locations,  an  extremely  large  adverse  pressure  gradient  exists  in  the 
diicciion  of  incieasing  jZ/Tj.  As  the  flow  deceieraies,  in  me  presence  of  a  more  mild 
adverse  pressure  gradient  in  the  streamwise  direction,  positive  values  of  Cp  occur  at  the 
trailing  edge.  As  oil-flow  visualizations  show,  separation  occurs  in  this  region  due  to  the 
effects  of  the  pressure  gradient  and  the  horseshoe  vortex.  Separation  does  not  occur  in 
the  2-dimensional  flow  region  toward  the  mid-span  of  the  appendage.  A  comparison  of 
the  wail  pressure  distribution  and  that  produced  by  an  unbounded  potential  flow  is  shown 
in  Figure  3.5.  While  the  qualitative  distribution  is  the  same,  the  magnitudes  of  the 
measured  Cp  disttibution  are  less  than  the  calculated  values  by  approximately  25%  to 
30%.  Figure  3.6  shows  a  comparison  of  the  actual  Cp  distribution  at  2/T  =  ±3.17  with 
that  predicted  from  potential  flow.  This  data  demonstrates  the  lack  of  significant 
blockage  effects  in  the  test  flow. 
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As  mentioned  in  section  2.1,  symmetry  of  the  appendage  boundary  layer  was  taken 
as  a  good  indication  of  the  flow  symmetry  (Figure  3.7).  Figures  3.8  and  3.9  show  the 
results  of  mean  velocity  and  turbulent  fluctuation  measurements  1  mm  downstream  of  the 
trailing  edge  at  various  Y/T  locations.  The  2-dimensionality  of  the  boundary  layer 
indicates  that  the  junction  flow  is  independent  of  the  height  of  the  wing  (Devenport  and 
Simpson,  1990a).  These  measurements  can  be  used  to  determine  what  features  of  the 
junction  flow  can  be  attributed  to  the  wing  boundary  layer  and  also  for  determining  the 
dissipation  characteristics  of  the  wake. 


3.1,3  Oil-Flow  Visualization  Results 


Oil-flow  visualizations  were  performed  by  previous  researchers  (Devenport  et  al., 
1990,  Devenport  and  Simpson,  1990b,  and  Olcmen,  1990),  and  Figures  3.10a,b  display 
the  results.  Note  that  the  individual  measurement  locations  are  indicated  by  the  white 
dots,  and  the  individual  XY  planes  are  also  indicated  in  Figure  3.10a.  To  produce  the  oil- 
flow  visualizations,  a  mixture  of  kerosine,  titanium  dioxide  powder,  and  oleic  acid  w'as 


spread  on  sheets  of  adhesive  glossy  black  plastic  sheeting 


VV-/1W111V11,  1 \j5ijig  lui:^ 


technique,  the  original  oil  flows  could  be  preserved,  since  the  sheet  could  be  removed 


from  the  wind  tunnel  when  the  oil  had  dried. 


Examing  Figiur  3.10b,  one  can  identify  the  primary  separation  point  (saddle  point) 
and,  closer  to  the  appendage,  the  collection  of  pigment  along  what  can  be  called  a  line 
of  low  shear  (LOLS)  (Devenport  and  Simpson,  1990b).  This  line  is  not  a  separation  line 
since  oil  streaks  pass  through  it  near  the  plane  of  symmetry.  This  is  in  contrast  to 


Description  of  Test  Conditions 


39 


previous  researchers’  conclusions.  The  4-vortex  structure  deduced  by  Baker  (1980)  and 
others  may  exist  instantaneously,  but  time  mean  velocity  measurements  along  the 
centerline  upstream  of  the  vtdng  show  only  one  vortex  in  the  plane  of  symmetry 
(Devenport  and  Simpson,  1990b).  It  can  be  postulated,  however,  that  the  LOLS  location 
and  Baker’s  Xy  location  are  directly  related,  since  away  from  the  plane  of  symmetry,  the 
LOLS  more  closely  resembles  an  ordinary  separation  line.  This  theory  is  supported  by 
the  fact  that  for  the  thin  and  thick  boundaiy  layer  cases  documented  by  Devenport  et  al. 
(1990),  the  variation  in  Xj  and  X^ols  fodow  the  same  trends  as  Belik  (1973)  and  Baker 
found  for  Xs  and  Xy.  Note  that  the  oil-flow  visualization  does  show  evidence  of  a  small, 
counter-rotating  vortex  very  close  to  the  base  of  the  wing,  suggesting  a  mean  two  vortex 
flow  structure. 

At  the  appendage  nose,  an  ordinary  separation  line  emanates  from  the  singular 
separation  point.  Figui’c  3.11  presents  a  sketch  showing  how  the  initial  vortex  roll-up 
causes  this  line.  Also  note  the  skewing  of  the  outer  streamlines  shown  in  this  figure.  As 
this  initial  voritical  structure  wraps  itself  around  the  appendage,  the  separation  line  and 
initial  LOLS  merge  to  form  what  appears  to  be  an  ordinary  separation  line,  at  least  in  the 
region  adjacent  to  the  appendage.  This  blending  suggests  that  the  '  vortex  structure 
present  at  the  nose  is  of  no  consequence  downstream  of  the  wing  maximum  thickness, 
since  any  individual  vortices  (aside  from  the  small  comer  vortex)  combine  to  form  one 
primary  vortex.  At  the  appendage  ti-ailing  edge,  a  "fish  tail"  wake  pattern  occurs.  This 
separation  pattern  is  caused  by  the  pressure  recovery  at  the  trailing  edge,  and  the 
secondary  flow  structure  helps  to  shape  it.  Figure  3.12  gives  the  dimensions  of  the 
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characteristic  oil-flow  pattern,  and  Table  3.2  lists  the  locations  of  the  LOLS  for  planes 
5-10. 


3.2  Comparison  of  Incoming  Boundary  Layer  Test  Parameters 


One  of  the  main  objectives  of  this  research  is  the  detennination  of  the  effects  of 
varying  inflow  conditions  on  the  properties  of  a  turbulent  wing-body  junction.  The  proper 
comparison  of  inflow  conditions  for  the  various  data  sets  becomes  a  key  issue  because 
of  this  objective. 

Approach  boundary'  layer  data  from  Dickinson  (1986a, b),  McMahon  et  al.  (1987), 
and  Dcvcnport  and  Simpson  (1990a)  was  used  to  evaluate  the  incoming  boundary  layer 
test  parameters.  Note  that  the  boundary  layer  parameters  for  the  author’s  data  is  exactly 
the  same  as  for  Devenport  and  Simpson’s  data.  Incoming  boundary  layer  data  for  Shin’s 
(1989)  research  was  unavailable  to  the  author,  and  the  given  values  will  be  used  for  his 

Figure  3.13  shows  the  approach  boundary  layer  profiles  for  the  three  main  test 
conditions.  This  data  was  taken  with  the  wing  in  place.  Devenport ’s  and  McMahon’s 
data  is  nearest  the  wing,  and  therefore  is  affected  by  the  adverse  pressure  gradient  caused 
by  tlie  wing  geometry  more  than  Dickinson’s  data.  Since  Dickinson’s  and  McMahon’s 
profiles  do  not  include  data  below  U/U^^f  =  ,5,  Spalding’s  wall  law  (equation  (3.2))  was 
used  to  estimate  additional  data  points  closer  to  the  wall  to  obtain  more  accurate 
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momentum  and  displacement  thicknesses  (as  discussed  in  section  4.1,3).  By  fitting  this 
equation  to  the  data*  an  estimate  of  Q  is  also  found.  Figures  3.14  (a-c)  show  the  results 
of  using  this  equation  to  extrapolate  additional  data  points  down  into  the  laminar  sublayer. 
For  each  profile,  4  additional  points  were  calculated,  and  are  indicated  on  the  figures  by 
the  solid  squares.  Devenport’s  data,  taken  v/ith  a  hot-wire  held  horizontal  to  the  floor, 
shows  quite  good  agreement  with  equation  (3.2).  Dickinson’s  data  seems  to  agree  well 
also,  but  suffers  slighdy  from  having  few  points  in  the  inner  region.  A  good  fit  to 
McMahon’s  data  using  equation  (3.2)  could  not  be  found,  so  the  Q  value  which  matched 
the  slope  of  tlie  linear  portion  of  the  inner  region  was  used  to  extrapolate  data  closer  to 
the  wall.  Several  possible  explanations  exist  as  to  why  McMahon’s  data  and  Spalding’s 
equation  do  not  agree  well.  The  data  was  obtained  using  a  digitizer,  which  could  lead 
to  some  inaccuracies.  The  size  and  design  of  the  probes  may  have  been  a  problem  near 
the  wall,  since  one  of  the  probes  was  a  slant  wire  with  its  supports  protruding  from  the 
floor,  a  design  which  could  have  possibly  biased  the  measurements.  Ine  combination  of 
a  thick  boundary  layer  and  the  close  proximity  of  the  wing  creating  an  adverse  pressure 
gradient  could  invalidate  equation  (3.2)  near  the  wall.  The  latter  reason  seems  the  most 
likely.  Ludwicg-Tillmann’s  skin-friction  correlation,  given  by  White  (1974),  was  used 
to  check  for  internal  consistency  of  the  Q  results.  The  Ludwieg-Tillmann  Q  results  were 
very  near  those  calculated  from  Spalding’s  equation,  with  errors  of  4%  for  Devenport’s 
and  Simpson’s  data,  and  2%  for  both  Dickinson’s  and  McMahon’s  data. 

After  fitting  equation  (3.2)  to  the  three  profiles,  5*  and  9  were  calculated.  The 
extra  points  increased  tliese  thicknesses  appoxirnately  5%  for  McMahon's  and  Dickinson’s 
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data,  but  the  changes  in  Devenport’s  results  were  negligible.  McMahon’s  and  Dickinson’s 
profile  results  were  corrected  to  X/T  =  -2.146  to  match  Devenport’s  location.  Power-laws 
were  used  to  correct  the  data.  The  relationships  (from  White,  1974)  used  were 

Cf  =  .0592  Rex'  ^ 

8/X  «  .37  Rex  ^  (3.3) 

G/X  -  .036  Rex'  ^ 

where  X  is  measured  from  a  virtual  origin.  The  shape  factor,  H,  was  assumed  constant 
and  used  to  evaluate  5*.  Table  3.3  lists  the  approach  boundary  layer  parameters  for  all 
data  sets  considered  in  this  study.  Two  important  facts  to  note  are  (1)  tliat  in  Shin’s 
investigation,  the  approach  boundary  layer  parameters  were  measured  at  X  =  0,  with  no 
wing  in  place,  and  (2)  that  Devenport’s  data  at  XJC  =  3.00  (plane  12)  is  from  what  was 
called  the  thin  boundary  layer  case  (Dcvenport  etal.,  1990),  those  parameters  being  listed 
in  the  last  column.  The  first  column  lists  approach  boundary  layer  comparisons  for  the 
author’s  (JUF)  hot-wire  data,  Devenport  and  Simpson’s  (WJD)  LDV  data,  and  Devenport 
et  al.’s  (1990)  hot-wire  data  for  the  thick  boundary  layer  case. 
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4.0  Experimental  Results 


This  chapter  presents  tlie  experimental  results  of  the  author’s  research.  The  like 
quantities  are  grouped  together,  with  the  contour  plots  and  secondary  velocity  vector  plots 
shown  first,  followed  by  selected  individual  profile  plots  which  enable  a  more  detailed 
examination  of  the  results.  Because  c.''  the  large  volume  of  results,  no  tabular  listings  are 
given.  The  data  files  are  available  on  magnetic  disk,  an,d  Appendix  B  describes  their 
organization. 

As  mcntionfxl  previously,  tlie  author’s  experimental  results  include  Ujc,  w'ye, 
^'tc  information.  In  addition  to  these  fundamental  quantities, 
(3(W)/5Y)6/ll,jf,  heiicity  density  (h),  the  distortion  function  (Id  ),  5*,  9,  and  Cf,  the  skin- 
friction  coefficient  based  on  U  are  also  presented.  The  calculation  of  these  quantities 
is  described  in  the  next  section.  Hodograph  plots  axe  also  presented,  displaying  the 
relationship  between  W,;g  and  Upg. 

Note  that  all  data  were  taken  on  the  starboard  side  of  the  appendage.  The  contour 
plots  of  data  in  a  YZ  plane  are  or  .nted  as  if  the  observer  were  looking  upstream,  with 
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the  appendage  on  the  left  side,  and  -Z/T  increasing  to  the  right.  All  contour  and 
secondary  velocity  plots  arc  scaled  to  1  inch  =  1.0  T. 


4.1  Description  of  Calculated  Quantities 

The  following  quantities  discussed  were  all  estimated  from  experimental  data. 

4.1.1  Streamwise  Vorticity  (d(W)faY)  and  Helicity  Density  (h) 

The  streamwise  vorticity  was  calculated  since  it  is  a  direct  indication  of  the 
strength  of  the  secondary  flow.  The  streamwise  vorticity  is  given  by 

Q  (4.1) 

*  sr  az 

where  ail  quantities  are  in  tunnel  coordinates.  The  differentials  were  estimated  by  fitting 
a  parabola  through  5  consecutive  data  points  using  a  least  squares  method,  and  then 
e'Mluating  the  derivative  of  this  curs'e  fit  at  the  central  (third)  point.  For  the  author’s 
data,  only  the  contribution  due  to  W  could  be  calculated,  rs  no  V  information  was 
available.  This  problem  was  not  serious,  however,  because  the  mean  secondary  flow 
structure  in  the  junction  tends  to  be  flattened  into  an  elliptical  shape  near  the  wall 
(Dickinson,  1986a, b).  Within  the  secondary  flow  region,  the  contribution  of  the  W  term 
to  is  usually  at  least  an  order  of  magnitude  iarger  than  the  V  contribution.  Figure  4.1 
illustrates  this  imbalance  by  comparing  d(W)/dY  and  -3(V)/^Z  in  planes  5  and  8  for  data 
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sets  where  V  data  is  available  (Devenport  and  Sirnpson,  1990a,  and  Dickinson,  1986a,b). 
The  W  term  has  the  larger  contribution  by  far  in  both  planes,  to  the  extent  that  it 
characterizes  the  vorticity  inherent  in  die  secondary  flow  structure  by  itself.  For  this 
reason,  3(W)/dY  was  the  quantity  plotted  in  contour  and  profile  plots,  enabling  direct 
comparisons  between  all  the  data  sets.  Since  the  term  is  proponional  to  the 

voritici  y  present  in  die  approach  boundary  layer  data,  all  vorticity  data  was  non- 
dimensionalized  by  dividing  through  by  this  term,  unless  otherwise  specified.  Figure  4.1 
V  as  non-dimensionalized  in  this  way. 

Hclicity  is  dircedy  related  to  the  vorticity,  and  is  defined  by 

H  ^  J  V  •  d(Volume) 

The  integrand  of  equation  (4.2)  is  known  as  the  helicity  density,  and  is  designated  h. 
Hussain  (1986)  discusses  the  importance  of  helicity  in  the  study  of  coherent  structures. 
He  explains  that  when  structuies  with  large  values  of  helicity  are  present,  they  will  have 
low  dissipation  rates  and  the  breakup  of  large  scale  into  smaller  scale  motions  is  slowed. 
This  occurs  as  a  consequence  seen  in  a  form  of  the  Navier-Stokes  equation,  and  the 
reader  is  referred  to  Hussain  (1986)  for  further  reading. 

With  regards  to  the  author’s  work,  the  non-dimensional  helicity  density  was 
csdmated  for  all  data  sets  by  letting 
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h 


,  ar '  dY 


6 


(4.3) 


This  approximation  results  from  ignoring  the  mean  velocity  partial  derivatives  in  the  X 
and  Y  directions.  In  a  fashion  similar  to  i2x>  ^  non-dimensionalized  by  5/(Urgf)^. 
Larger  values  of  helicity  density,  from  the  preceeding  discussion,  will  indicate  longer 
lasting  vortical  structures,  with  less  dissipation  due  to  turbulent  momentum  transport. 
Helicity  density  contours  are  presented  in  chapter  5. 


4.12  Distortion  Function  (f^) 

The  distortion  function  was  used  to  separate  the  effects  of  the  secondary  flow 
structure  from  the  effects  of  the  incorning  boundary  layer  and  of  the  appendage  geometry. 
The  distortion  function  is  defined  by  the  equation 


for  a  constant  X/C  value.  Ug,  the  local  boundjuy  layer  edge  velocity,  was  predicted  by 
a  potential  flow  calculation.  These  poieriiiai  flow  caicuiations  compare  favorably  to  the 
experimental  results,  as  can  be  seen  in  Figure  4.4.  Tlie  function  g  is  defined  by  the  outer 
’'2-D"  mean  flow  distribution  in  the  boundary  layer  in  the  given  YZ  plane  location.  One 
can  see  that  f^  describes  the  effects  of  the  horseshoe  vortex  on  the  mean  streamwise  flow 
distribution,  and  its  value  should  go  to  1  at  large  Z/T  or  Y/T. 
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4.1.3  5*,  0,  and  Estimates 

Results  from  estimating  Q,  8*,  and  0  were  used  to  compare  and  contrast  the 
effects  of  the  horseshoe  vortex  structure  in  the  different  data  sets.  These  results  are 
discussed  in  section  5.3.4.  All  results  were  calculated  in  the  local  freestream  coordinate 
system  ((We)ps  =  0). 

The  local  boundary  layer  displacement  and  momentum  thicknesses  in  the  local 
freestream  direction  are  given  by 


A1 

T 


r'Jo  I,*' 


(4.5) 


fA  A\ 


These  expressions  were  integrated  using  the  trapezoidal  rule.  The  integrand  of  equation 
(4.5)  is  a  maximum  at  Up^E  “  so  using  the  trapezoidal  rule  to  evaluate  5*  for  profiles 
lacking  near  wall  data  was  considered  a  reasonable  approximation.  The  integrand  of 
equation  (4.6)  reaches  a  maximum  at  U/Ue  =  .5,  and  goes  to  zero  as  Ups,TJe  approaches 
zero.  Therefore,  using  the  trapezoidal  rule  to  evaluate  0  for  profiles  where  (Ups/T1i:)min 
>  .5  can  result  in  significant  en'ors.  When  profile  data  did  not  meet  die  (UitsAJe)min 
.5  criteria,  the  Ups/Ue  =  .5  location  was  roughly  estimated  using  linear  interpolation,  and 
then  equation  (4.6)  was  used  to  estimate  0  for  all  data  sets. 

Several  methoc's  were  used  to  estimate  Q  in  the  local  freestream  direction  ((Cf)i.s) 
for  each  profile.  Johnston  (1960)  discusses  a  technique  involving  what  he  refers  to  as 
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hodograph  or  polar  plots  of  WpsAJ^f  vs.  Ups/U„t  to  obtain  a  wall  law  relationship  for  3- 
dimensional  turbulent  boundary  layers.  He  suggested  the  form 


U* 


=  -i- In(y")  +  5.0 

.41 


(4.7) 


where  Kwall  is  measured  with  respect  to  the  local  freestream  coordinate  system,  and 
is  defined  using  the  magnitude  of  the  wall  shear  stress  vector.  This  angle  can  be  found 
by  evaluating 


_ ^ 

V  ^FSJ 


(4.8) 


as  Ups  approaches  zero,  and  is  presumably  also  the  direction  of  the  local  wall  shear  stress 
wciwi,  j.  igaiu  i-iwiiji  vYiiiic  iiiuiiuiiicb  uic  uixjvc  rciauonsnip,  Jonnscon  reports 

some  success  using  equation  (4.7)  for  boundary  layers  with  small  secondary  flows. 

Unfortunately,  the  hodograph  plot  idea  and  Johnston’s  wall  law  breaks  down  when 
bidirectional  or  "cross-over"  profiles  are  encountered,  as  shown  in  Figure  4.2,  or  for 
extreme  flow  angles.  The  horseshoe  vortex  creates  such  bidirectional  polar  plots.  The 
skewing  created  by  liic  strung  secondary  now  made  estimates  of  Q  obtained  using 
equations  (4.7)  and  (4.8)  highly  erratic.  The  results  at  some  locations,  usually  well  away 
from  the  wing,  were  encouraging,  but  overall  results  were  unacceptable, 

The  other  methods  used  to  estimate  Cj  in  the  local  freestream  direction  were  two 
skin-friction  correlation  formulas,  given  by 
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(4.9) 


Cj  =  0.246 


0.3c 


-;.33ff 


Oogioi^e) 


(1.74t0.31«> 


(4.10) 


where  H  is  the  shape  factor.  The  former  equation  is  the  classical  Ludwieg-Tillmann 
relation,  while  the  latter  is  a  skin-friction  correlation  given  by  White  (1974)  using  Coles’ 
wall  law  coefficients.  These  two  correlations  are  considered  accurate  only  for  mean  2- 
dimensional  flows,  but  were  used  to  provide  "last  resort"  estimates  to  facilitate  at  least 
a  qualitative  comparison  of  skin-friction  results  between  data  sets.  Note  that  (Cf)ps  can 
be  used  to  find  the  total  skin-friction  coefficient  if  is  known. 

Figure  4.3  shows  an  example  of  the  typical  Q  estimates  using  the  different 
equations.  Note  that  the  results  from  equation  (4.9)  and  (4.10)  are  almost  identical. 
Equation  (4.10)  and  the  local  velocity  profile  properties  were  used  to  generate  Cj  data  for 
the  local  streamwise  component  of  the  wall  shear  stress.  These  results  are  compared  to 
Ailinger’s  (1990)  Q  results  obtained  using  laser  interferometry,  and  are  discussed  in 
section  5.3.6.  Since  the  Q  estimates  from  equation  (4.10)  are  based  on  the  local 
freestream  dynamic  pressure,  the  skin-friction  coefficients  were  multiplied  by  (Uj/U,,f)^ 
to  non-dimcnsionalize  Cf  on  the  reference  dynamic  pressure. 
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4.2  Near  Wing  Data  Corrections 


A  2-D  potential  flow  solution  was  obtained  for  the  appendage  shape  studied  using 
the  vortex  panel  solution  described  in  chapter  6.  The  comparison  of  the  author’s 
experimental  results  to  the  2-D  numerical  solution  (Figure  4.4)  led  to  corrections  of  the 
flow  angle  datum  in  planes  9,  10,  and  11. 

Plane  5  comparisons  (Figure  4.4a)  show  good  agreement  for  freestream  values  of 
UpsAlref-  Th®  author’s  freestream  values  of  Ojc  do  not  compare  very  well,  and  are 
approximately  1.5°  to  2.0°  off  at  most  locations.  Initially,  a  correction  was  thought  to  be 
in  order,  but  after  examination  of  the  profile  plots  of  WTc/U^jf  in  plane  5  (Figure  4.32), 
One  can  see  that  ctpc  never  reaches  a  constant  or  what  could  be  considered  a  freestream 
value.  Considering  this  fact,  no  corrections  could  be  justified. 

In  planes  6  and  7  (Figures  4.4b, c),  good  overall  agreement  of  both  the  freestream 
velocity  and  flow  angle  was  observed,  and  no  corrections  were  needed.  In  plane  9 
(Figure  4.4e),  the  flow  angles  were  corrected  by  -2°  for  the  outer  7  profiles,  these 
pfOuics  weie  measured  almost  2  weeks  after  the  iniicr  3  profiles  (which  were  not 
coi  ected),  and  show  a  systematic  bias  error  in  the  probe  alignment  of  -2°.  Figure  4.4e 
shows  the  corrected  results. 

Comparisons  in  planes  10  and  11  (Figures  4.4f,g)  show  good  agreement  for  the 
freestream  velocities,  but  the  author’s  experimental  freestream  Otc  results  show  some 
scatter.  These  results  also  showed  a  possible  bias  error  for  all  profiles  in  these  two 
planes.  Tests  for  symmetry  in  plane  11  showed  that  a  correction  of  -1.5°  would  lead  lo 
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symmetric  flow  angles,  and  repeatibility  tests  showed  the  (Xj-^  measurements  to  be  precise. 
Therefore,  a  -1.5°  correction  in  Opc  was  made  for  planes  10  and  11. 

The  effects  of  the  corrections  in  flow  angle  for  planes  9,  10,  and  11  aie  visible 
only  in  and  otpc;.  The  other  quantities  measured  were  effected  very  slightly,  which 
is  expected,  since  the  cosine  of  1.5°  is  .99966. 


4.3  Examination  of  the  Present  Flow  Characteristics 

The  results  of  the  author’s  measurements  are  discussed  in  this  section.  Devenport 
and  Simpson’s  (1990a)  LDV  data  and  Devenport  et  al.’s  (1990)  data  from  the  thick 
boundary  layer  case  is  used  for  continuity  from  planes  5-15.  The  LDV  data  also  reveals 
additional  detail  since  it  contains  V  velocity  component  information.  Appendix  B 
describes  the  organization  of  the  data  on  magnetic  disk.  The  plots  are  organized 
according  to  the  quantity  being  presented,  with  the  contour  and  velocity  vector  plots 
shown  flrst,  followed  by  boundary  layer  profile  plots  from  selected  X/'C,  Z/T  locations. 

The  inclusion  of  individual  boundary  layer  profile  plots  presents  the  data  in  such 
detail  that  to  fully  discuss  the  results  is  beyond  the  primary  objectives  of  the  autlior’s 
research.  These  plots  are  included  primarily  for  the  interested  reader,  and  therefore  will 
not  always  be  nienrioned.  Ai  umes,  however,  the  profile  plots  show  needed  detail 
unavailable  from  the  contour  plots,  particularly  near  the  wall,  and  they  also  display  results 
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not  presented  in  the  contour  plots.  In  chapter  5,  the  profile  plots  aid  in  the  comparison 
of  data  sets,  as  well  as  the  determination  of  the  flow  scaling  parameters. 

4.3.1  Near  Appendage  Data,  Planes  5-11 

Figures  4.5a  and  4.6a  display  the  contours  of  mean  velocity  and  distortion 
function,  while  Figure  4.7  shows  the  secondary  velocity  vectors  where  available.  The 
term  "secondary  velocity  vector"  is  used  loosely  here  in  describing  the  vectors  formed  by 
the  Wtc  and  Vtc  data.  Ideally,  the  velocity  vectors  should  be  presented  in  a  way  that 
distinguishes  the  actual  secondary  flow  due  to  the  horseshoe  vortex  from  the  spanwise 
flow  due  the  appendage.  Presentation  of  the  vectors  in  a  plane  orthogonal  to  the  vortex 
axis  (or  line  of  low  shear)  or  in  freestream  coordinates  may  be  more  appropriate  than 
tunnel  coordinates,  which  is  the  coordinate  system  used.  This  system  was  chosen 
primarily  to  enable  clear  comparisons  with  other  data  sets  (chapter  5). 

Contours  of  U-j-cAJfef  for  planes  5-11  show  the  significance  of  the  line  foraied  by 
the  merging  of  the  line  of  low  shear  location  (LOLS)  and  the  separation  line.  It  is 
designated  by  the  dashed  line  in  Figure  4.5a,  while  the  solid  line  indicates  the  appendage 
surface  location.  A  "tongue"  or  plume  of  low  momentum  fluid  grows  as  the  flow  moves 
past  the  junction.  This  upward  distortion  of  the  boundary  layer  is  always  near  the  LOLS. 
The  secondary  velocity  vector  plots  which  contain  V  infonnation  show  the  maximum 
vertical  velocities  near  the  LOLS,  which  also  correspond  to  the  distortion  function  contour 
patterns.  From  examining  the  relationship  between  Uj^  and  the  secondaty  velocities,  one 
can  easily  see  the  characteristic  pumping  action  of  the  horseshoe  vortex  at  work,  a 
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mechanism  which  brings  higher  momentum  fluid  into  the  junction,  and  ejects  lower 
momentum  fluid  from  the  near  wall  region.  Cutler  and  Bradshaw’s  (1986)  experimental 
results  display  the  same  general  features,  with  an  initially  large  crossflow  leading  to  the 
transport  of  fluid  away  from  the  wall. 

The  tongue  of  low  momentum  fluid  extending  from  the  near-wall  boundary  layer 
changes  in  shape  as  X/C  increases,  growing  narrower  as  the  appendage  trailing  edge  is 
approached,  and  then  broadening  aft  of  the  trailing  edge  separation  location.  The 
contours  of  fj,  change  from  a  more  vertical  orientation  and  begin  to  tilt  inwards  toward 
the  centerline  as  the  distortion  pattern  broadens. 

The  flow  created  around  an  appendage  is  subjected  to  skewing  and  hence  should 
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embedded  in  a  skewed  boundary  layer  created  by  deflecting  the  freestream  flow  30°  in 
the  transverse  di  ection,  and  examined  the  cases  resulting  from  the  two  possible  vortex 
rotational  directions.  It  was  found  that  the  vortex  decayed  much  quicker  in  a  skewed 
rather  than  2-D  boundary  layer  case,  and  that  the  vortex  persisted  longer  and  produced 


loTj^ci  flow  diStOiiiOuS  ifi  the  uuSc  whcic  ilic  skewed 


fluid 


vciOciiy  near  liie  wan  was  or 


opposite  sense  to  the  flow  induced  by  the  embedded  vortex.  This  flow  was  designated 
case  2  by  Shizawa  and  Eaton.  Tlie  distortions  produced  by  the  opposite  case,  designated 
as  case  1,  decayed  the  quickest  because  the  near-wall  secondary  flow,  moving  in  the  same 
direction  as  tlie  vortex  induced  flow,  suppressed  separation  and  inhibited  the  creation  of 
the  characteristic  plume  of  low  momentum  fluid.  Shizawa  and  Eaton  observed  that  fluid 
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containing  vorticity  of  opposite  sign  was  able  to  pass  underneath  tlie  remains  of  the  case 
1  vortex  core  as  the  flow  developed  downstream. 

ITie  region  upstream  of  the  trailing  edge,  in  planes  6-8,  and  also  in  planes  9  and 
10  to  a  lesser  extent,  exhibits  characteiistics  of  Shizawa  and  Eaton’s  case  2  flow.  The 
secondary  velocity  vectors  clearly  show  the  vortex-induced  skewing  opposing  the 
boundary  layer  skewing  caused  by  the  appendage  geometry.  These  conflicting  flow 
patterns  merge  to  form  an  energetic  separation,  as  seen  by  the  V  velociues  ii  plane  8 
(X/C  =  .75).  This  accounts  for  the  more  upright  distortion  function  contours  adjacent  to 
the  wing. 

Case  1  simulates  the  wake  region  of  a  junction  flow,  and  produces  distonion 
patterns  in  mean  velocity  similiar  to  die  present  experiment  Plane  11,  at  X/C  =  1.50 
shows  very  similiar  secondary  velocity  and  mean  velocity  distortion  patterns.  These 
observations  suggest  that  the  same  vortex  decay  mechanisms  are  at  work  in  the  two  flows, 
both  adjacent  and  downstream  of  the  body,  even  though  me  skewing  for  Shizawa  and 
Eaton’s  flow  is  much  stronger.  Profile  plots  of  vorticity  from  plane  1 1  (Figure  4.90), 
showing  near-wall  regions  of  negative  vorticity,  provide  further  evidence  of  similiar 
vortex  decay  mechanisms  at  work.  The  adverse  pressure  gradient  encountered  at  the 
appendage  trailing  edge  also  encourages  the  broadening  of  the  flow  distortions  and 
diffusion  vortical  flow  stmetures. 

Looking  again  at  the  planes  adjacent  to  the  appendage,  a  region  of  negative 
voracity  can  be  seen  near  the  wall  outboard  of  die  primary  vortex  stiaicture  in  Figure  4.8 
and  in  the  profile  plots  of  vorticity.  This  region  is  larger  than  that  observed  by  Cutler  and 
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Bradshaw  (1986)  for  a  2-D,  unskewed  boundary  layer,  but  not  as  large  as  seen  for 
Shizawa  and  Eaton’s  case  2.  These  two  studies  show  that  the  magnitude  of 
"background"  boundary  layer  skewing  plays  a  major  role  in  determining  the  behavior  of 
any  embedded  vortices. 

One  must  be  careful  about  viewing  a  wing-body  junction  flow  as  a  vortical  flow 
structure  superimposed  over  a  "background"  turbulent  boundary  layer  flowing  around  an 
appendage.  These  two  features  of  a  junction  flow  are  inseparable,  and  both  influence 
each  other.  To  illustrate  this  fact,  consider  the  following  relationships.  Without  an 
approach  boundary  layer,  the  appendage-generated  pressure  gradients  could  not  create  a 
root  junction  vortex.  The  resulting  vortical  flow  structure  affects  the  boundary  layer  mean 
flow  distribution,  which  in  turn  changes  the  wall  pressure  gradients,  which  can  then  affect 
the  vorticity  flux  out  of  the  wall  (Agui  and  Andreopoulos,  1991). 

Examining  the  vorticity  distribution,  shown  by  contours  of  (3(W)/9Y)6AJref  in 
Figure  4.8,  one  observes  that  the  vortex  is  ellipdcally  shaped  and  near  the  wall.  This  is 
similiar  to  Dickinson’s  (1886a)  earlier  findings.  The  initial  time-averaged  elliptical  shape 
of  tile  iTidxiiniini  vorticity  contours  in  pluncs  5-8  suggests  thnt  the  vortex  niiiy  be 
meandering.  Beginning  in  plane  9,  near  the  trailing  edge,  the  core  vorticity  contours 
begin  to  appear  more  circular  in  nature.  This  apparent  stabilizing  effect  may  be  due  to 
the  increase  in  distance  from  the  bimodal  flow  in  the  nose  region.  Note  tliat  the  LOLS 
is  always  located  outboard  of  the  vortex  core,  as  expected  from  the  previous  examination 
of  the  flow  separation  mechanism. 
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Comparing  Figures  4.5a  and  4,8,  one  can  see  that  the  vortex  core  locations 
coincide  with  regions  where  the  streamwise  mean  flow  gradients  in  the  spanwise  (Z) 
direction  are  large  (the  isovels  of  Ujq  are  almost  vertical).  This  observation  clearly 
demonstrates  the  relationship  between  the  redistribution  of  the  mean  velocity  in  the 
boundary  layer  and  the  vortex  posidon.  The  vonex  depresses  the  contours  of  Utc  by 
drawing  in  high  momentum  fluid  from  the  outer  boundary  layer,  displacing  low 
momentum,  near  wall  fluid  upward,  and  therefore  creating  large  spanwise  mean  flow 
gradients.  Table  4. 1  lists  the  maximum  levels  of  (3(W)/8Y)5/U„f  found  in  the  vortex 
core  at  each  measurement  plane.  The  decrease  in  vorticity  from  the  maximum  thickness 
location  to  the  trailing  edge  displays  evidence  of  the  diffusion  and  decay  occuiring  in  the 
vortex  core. 


Due  to  the  no  slip  condition  at  the  wall,  a  very  tiiin  region  of  negative  vorticity 
is  concentrated  underneath  the  outer  edge  of  the  positive  vorticity  region  near  the  wall. 
In  plane  5,  this  vorticity  is  confined  to  a  region  below  Y/T  =  .04  (Figure  4.84),  and 
extends  outwaid  from  approximately  the  LOLS  location.  As  the  flow  develops,  the  area 

OF  nponfive  cmHinllir  inr*Tv»rtc#ao  nnH 

T  Ai«w«wt^owo  LUIVJ.  lAiv  1.  A  1 1  ^  .t  *  IV/ LIV  A  IW  VY  iivlo  Vliv  liilUl  J,'' 

vonex  away  from  the  wall.  Similiar  behavior  is  observed  by  Culler  and  Bradshaw  (1986). 
Profile  plots  (Figures  4.84-90)  and  enlarged  contours  plots  of  (5(?tV)/c)Y)6/TJ,^f  (Figures 
5.27-5.39)  may  be  referred  to  for  details  of  the  negative  vonicity  region. 

The  characteristic  patterns  formed  by  contours  of  the  square  rewt  of  the  Reynolds 
normal  stress  are  shown  in  Figure  4.9a.  Note  how  the  local  peaks  in  turbulence 

intensity  are  centered  slightly  outboard  of  the  LOLS  location,  except  at  X/C  ■-  .93,  which 
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is  uiar  Lhe  trailing  edge  separation.  T  location  of  the  L.OLS  at  X/C  =  .93  is  a 

minimum,  therefoie  the  relationship  between  the  primiir}'  vortex  position  and  the  LOLS 
may  change  near  this  location.  Devenport  et  al.  (1990)  suggests  several  causes  for  the 
local  turbulence  peaks.  These  are  a)  increased  mixing  of  the  boundary  layer  fluid  by  the 
aorseshoe  vortex,  b)  spanwise  meandering  of  the  trailing  vortex  legs,  and  c)  a 
destabilizing  of  the  boandary  layer  caused  by  the  lifting  action  of  the  secondary  flow  near 
the  LOLS,  simiiiar  to  tlie  effects  of  an  adver.se  pressure  gradient.  The  first  and  third 
causes  seem  to  be  the  most  likely.  The  coincident  locations  of  tite  steep  spanwise 
gradients  in  Uxe  and  tiie  vortex  coves,  both  positioned  itibcard  of  the  LOLS,  would  seem 
to  discount  the  second  possible  caense  of  the  turbulence  peaks. 

The  individuaJ  profile  plots  show  the  effects  of  the  secondary  flow  s'xucture  in 


greater  detail.  The  kinlung  of  the  Ux,/U„.f  profiles  (Figures  4.10'4,23)  shows  the  effects 
of  adding  higher  momentum  fluid  near  the  wall,  and  displacing  the  lower  momentum  fluid 
further  out  of  the  boundary  layer.  As  the  flow  moves  downstream,  the  distortion  of  these 
profiles  becomes  much  greater  Plots  of  Yxe/U^ef,  available  from  Devenport  and 
Simpson’s  (1990a)  LDV  data,  show  local  minima  ^or  locations  inboard  of  the  rime- 
avemged  vortex  center,  and  local  maxima  outboard  of  the  vortex  .  enter  (I-igiires  4.29- 
4.31).  Also  note  that  the  V  velocities  tUways  approach  zero  near  the  wall,  even  nc-.r  the 
LOLS  location. 


Profile  plots  of  WxcA.^  4f  (Figures  4.32-4.38)  and  cr,;-  (I'igures  4. 39-4, 4,3)  show  the 
large  changes  the  spanw;.se  flow.  In  plane  5,  the  spanwise  skewing  exceeds  40%  of 
In  the  downstream  planes,  b  -directional  skewing  is  often  present,  particularly  inboard  of 
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the  LOLS  location.  I'he  complexity  of  the  3-D  boundary  layer  at  these  locations 
invalidates  triangular  flow  models  based  on  hodograph  plots  (Johnston,  1960).  The 
hodogtaph,  or  polar  plots,  are  displayed  ir.  Figures  4.77'4.83,  and  show  the  complex 
relationships  existing  between  Wjrs  and  Up^  in  the.  junction  vortex  flow  region.  At  larger 
values  of  |Z/T|,  however,  the  boundary  layer  skewing  often  resembles  Johnston’s  (1960) 
experimental  flow,  in  which  a  duct  end  wall  created  the  skewed  boundary  layer. 

Profiles  of  the  fluctuating  quantities  k^c.  v'xc>  are  given  in 

Figures  4.46-4.76.  The  profiles  of  v'^-c  and  uv-^-q  were  obtained  from  Devenport  and 
Simpson’s  (1990a)  LDV  data.  The  details  of  the  mean  flew  influence  on  the  distribution 
of  these  quantities  within  the  boundary  layer  are  displayed  in  these  figures.  Near  the 
LOLS,  the  profile  plots  usually  show  local  ab,sol!uc  maxima.  This  is  panicalariy  evident 
in  plane  5,  where  2-nd  w\q  .show  large  levels  near  7/T  =  -.88,  tfie  approximate  LOLS 
location.  The  large  values  in  occur  near  the  wall  (Y/T  <  .01),  while  smaller 
increases  in  w'  and  large  peaks  in  v'p,-  occur  at  YfY  of  approximately  .09.  This  location 
coinciaes  precisely  with  the  peaks  in  uVe  iri  plane  5.  As  the  flow  moves  downstream, 
she  distri'ouiions  of  the  normal  str6.ises  are  more  uniform,  but  small  peaks  are  still  present, 
and  seem  to  be  related  to  one  another,  'ITiis  appa.'ent  relationship  suggests  that  the  How 
mechanism  respon-sible  for  altering  the  distribution  of  u\q  affects  v'-j-c  and  w'-^c  in  u 
similiar  fashion. 

It  is  interesting  to  comp.'irc  the  norm.aI  stresses  in  plane  5  to  the  nomiruil  values 
occuring  for  a  2-D  flow  .over  a  flat  plate.  Klcbanoffs  classical  turbulent  boundary  layer 
over  a  flat  plate  (IVum  Schetz,  1984)  shows  maximum  values  of  n'/U|.;,  v'/Uj;  and  'v7l.l|. 
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to  be  around  .11,  .04,  and  .06,  respectively.  Allowing  for  the  fact  that  the  present 
experimental  data  are  non-dimensionalized  by  U„f,  not  Ug,  the  values  for  »'xc/Ug  near  the 
wall  are  approximately  tiie  same  as  Klebanoffs  results,  but  die  other  two  turbulent 
quantities  are  much  nigher.  The  maximum  values  of  present  in  plane 

5  are  well  over  twice  those  given  by  Klebanoffs  data.  Presenting  the  data  in  freestream 
coordinates  does  not  affect  this  comparison  with  2-dimensional  data.  As  X/C  increases, 
the  local  maxima  in  v'^-c  and  w'^c  decrease  to  more  comparable  levels,  and  by  plane  8, 
these  maximum  values  are  approximately  .05  and  .07,  respectively. 

The  mean  Reynolds  shearing  stresses  and  MWjc/CUrcf)^  profile  plots  are 

shown  in  Figures  4.67-4.76.  Tlie  effects  of  the  altered  mean  flow  gradients  are  revealed 
in  these  figures.  Large  negative  values  of  uw-yc  exist  underneath  the  mean  vortex 
structure  because  of  the  spanwise  gradient  in  Uxc  at  Jus  location,  and  also  because  of  the 
large  transverse  velocities  near  the  wall.  The  large  vaiues  in  tovxc  n^iay  also  be  indicative 
of  vortex  meandering.  Positive  or  negative  regions  of  uvvxc  can  usually  be  explained  by 
examining  the  spanwise  gradient  in  Uxc-  Local  maxima  or  minima  seem  to  occur  in 
regions  marked  by  high  normal  stresses  as  well. 

Interestingly,  no  significant  regions  of  positive  (TP-jc  wcie  found  in  plane.s  5,  8,  or 
10.  Local  maxim,a  of  nvxc  which  tended  towai’d  zero  can  be  easi’y  identified  at  Z/T 
locations  strongly  effected  by  the  horseshoe  vortex. 

Figures  4.91  and  4.92  show  the  similiar  variations  of  and  O/Qait*  in  the 

spanwise  direction  for  each  of  die  VZ  measurement  planes.  Those  arc  freestrearn 
coordinate  system  results.  Planes  5  and  6  display  a  shaqi  drop  in  6'  and  0  as  the  body 
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is  approached,  which  is  a  combination  of  the  thinning  effects  of  the  near  wing  pressure 
gradients  and  the  developing  horseshoe  vortex.  As  |Z/T|  increases,  the  boundary  layer 
thicknesses  are  relatively  constant  in  these  planes.  Beginning  with  plane  7,  tlie  variation 
of  5*  and  9  show  local  maxima  near  the  LOLS  locations.  This  result  is  expected  from 
earlier  observations  of  the  streamwise  mean  velocirv  contours.  The  local  spanwise 
variations  in  6*  and  0  increase  monotomically  with  downstream  distance.  Figure  4.94 
shows  the  vtiriation  in  A(5‘/T)  and  A(6/T)  vs.  X/C.  These  parameters,  which  are  a 
measure  of  the  extent  of  the  distortive  effects  due  to  the  horseshoe  vortex,  increase  as  the 
flow  moves  downstream.  Tlie  initial  distortion  rate  grows  very  quickly,  but  levels  off  to 
an  almost  constant  growth  rate.  Tlie  rapid  distortion  growth  is  an  indication  of  the 
skewing  direction  of  the  primary  flow  over  the  aft  portion  of  the  body,  which  tends  to 
increase  the  di.stortions  (case  2,  Shizawa  and  Eaton,  1990).  The  reversal  in  the  primary 


flow  skewing  direction  decreases  the  rate  of  growth  of  A(5*/T)  and  A(6/T)  past  X/C  =  1 . 
The  initximum  changes  in  5'/T  and  0/T  for  planes  5-12  are  listed  in  Table  4.1. 

Lceal  Cf  estimates  were  calculated  using  equation  (4.10),  and  the  results  are  shown 


!n  Figure  4.93.  These  results  are  presented  using  the  frccstrcam  coorduiatG  systcrri.  Fiorn 


the  form  of  equation  (4.10),  the  variation  in  Cf  is  expected  to  decrease  when  0  increases. 
This  is  in  fact  what  happens,  with  a  local  minima  occuring  near  the  oil-flow  visualization 
LOLS  location.  Of  course,  it  is  the  local  shear  minima  that  causes  this  line,  and  ideally, 
the  (Cf)vii»^  and  LOLS  locations  should  coincide  precisely,  llic  expected  general  decrease 
in  skin  friction  magnitudes  is  observed  as  the  boundary  layer  continues  to  develop  in  the 
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4.32  Far  Wake  Data,  Planes  12~15 

This  section  discusses  tlie  author’s  data  taken  at  X/'C  locations  of  5.89,  9.14,  and 
11.56,  and  the  data  from  Devenpon  et  al'^  (1990)  tliick  boundary  layer  at  X/C  =  3.00. 
The  two  most  downstream  positions  are  in  the  adverse  pressure  gradient  section  (region 
3,  discussed  in  section  3.1.1).  Only  UTcA^re>  measurements  were  taken  in 

these  planes.  The  data  planes  spanned  0  <  -2/1  <  4.5,  which  is  a  large  region  considering 
that  the  half  width  of  tlie  wind  tunnel  is  approximately  6  wing  thicknesses. 

Contours  of  UjcAIrsf  and  wTUref  planes  12^1 5  arc  given  in  Figures  4..5b  and 
4.9b.  Tlie  continuing  distortion  patterns  arc  well  defined  in  these  far  wake  planes.  It  is 
interesting  to  note  that  for  the  Los  Angeles  (688)  class  and  Ohio  class  submarines,  the 
propulsor  planes  are  located  at  approximately  X/C  -  9  and  11,  respectively. 

From  X/C  -  1.05  to  9.14,  the  vortex  induced  m  an  flow  distortions  move  away 
from  the  centerline  due  to  the  induced  flow  of  the  wall  image  vortex.  From  X/C  =  9.14 
to  X/C  =  11.56,  tlic  spanwise  position  of  the  di.stortion  patterns  changes  only  slightly,  if 
at  all.  The  vonical  flow  .structure  appears  to  be  quite  diffused  in  the  far  wake,  as  the 
contours  of  w'-jc/Ure;  do  not  show  the  local  ma.sima  "islands"  as  in  plane  11.  Further 
evidence  of  tlic  vortex  decay  and  diffusion  is  shown  by  the  changing  .shape  of  the 
contours  of  w‘Vc/C,tf.  Tine  almost  horizontal  "tongue"  of  high,  turbulei.u  fluid  at  X/C  = 
3.C0  recedes  to  Utc  point  where  the  contours  of  u\q  are  almost  identical  in  shape  to  those 
of  U-rc.  Tlie  adverse  pressure  gradient  seems  to  accclenne  this  change.  Note  that  the 
contour  plots  also  show  the  effects  of  tlie  s’.dewall  and  ceiling  boundary  layci's. 
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Westphal  et  al.  (1985)  compared  vortices  embedded  in  turbulent  boundary  layers 
in  the  presence  of  constant  pressure  and  increasing  pressure.  Their  resulting  contours  of 
mean  streamwisc  velocity  show  a  more  defined  distortion  pattern  for  the  constant  pressure 
case.  They  observed  that  the  adverse  pressure  gradie:  t  increased  the  rate  of  vortex  core 
growth,  but  lowered  tlie  peak  vorticity.  The  decrease  in  peak  vonicity,  which  is  expected 
from  Helmholz’s  invnscid  vorticity  theorems,  explains  the  decrease  in  distortion  for  the 
adverse  pressure  gradient  case.  Comparing  the  contours  of  Uxc/Uref  f^  at  .X/C  =  3.00, 
5.89,  and  9.14  in  Figures  4.5b  and  4.6b,  one  can  observe  the  slight  diffusion  of  the 
horseshoe  vortex  distortion  pattern.  The  contours  of  f^  which  mark  the  plume  of  low 
momentum  fluid  become  more  vertical,  indicating  a  weakening  of  the  vortex's 
recirculatory  action.  Comparing  the  similiarities  of  the  contours  of  from  planes  14  and 
15,  one  can  see  that  the  increase  in  the  apparent  distortion  is  due  primarily  to  the  growing 
boundary  layer.  The  rapid  growth  of  the  boundary  layer  makes  it  difficult  to  determine 
if  the  distortion  pattern  is  decaying  more  quickly  due  to  the  adverse  pressure  gradient. 
The  lack  of  spanwise  movement  from  plane  14  to  15,  does  indicate  the  relative  weakness 
of  the  trailing  horseshoe  vortex  in  region  3.  Tlie  dominant  factor  controlling  the  growth 
of  die  contours  of  and  seems  to  !■)€  the  expanding  boundary  layer,  and  from  X/C 
™  5.89  to  11.56,  the  spanwise  dimension  of  the  contour  distortions  appear  to  grow  very 
little. 

It  is  of  irUci'est  to  note  that  Westphal  et  al.  also  objcrved  that  an  adverse  pressure 
gradient  seemed  to  increase  vortex  core  aspect  ratio,  meaning  tfutt  a  more  elliptic  vortex 
■vas  CiCalcd.  They  repon,  that  this  elliptic  shape  does  not  appc;ir  to  be*,  caused  by  core 
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meandering.  This  result  could  not  be  verified  by  the  present  experimental  data,  since  no 
information  was  obtained. 

The  diffused  effects  of  the  horseshoe  vortex  can  be  seen  in  the  slight  deformation 
of  the  UfcAJref  far-wake  profiles  (Pigures  4.26-4.28).  The  profiles  of  show  more 

kinking,  with  the  local  maxima  being  much  further  away  from  txie  wall  (Figures  4..54- 
4.56).  The  profile  kinking,  which  tends  to  produce  closed  contours  of  decays  rapidly 
in  the  adverse  pressure  gradient.  The  peak  magnitudes  of  tire  vc.locit)'  fluctuations  away 
from  the  wall  remain  approximately  constant  downstream  of  X/C  =  3.00,  with 
(«'To^ref)MAX  .08.  Notc  that  Figures  4..54-4.56  shows  the  original  vs.  the 

corrected  values  of  as  discussed  in  section  2.5.1. 

Figures  4.91b  and  4.92b  show  the  far-v/ake  variation  of  576‘app  and  9/0app  in  the 
spanwise  direction.  Examining  these  figures,  and  also  the  far-wake  contour  plots,  one  can 
discern  the  effects  of  another  region  of  secondary  flow  at  a  larger  value  of  |Z/r|.  One 
possible  source  of  the  smaller  secondary  flow  structure  is  a  flow  disturbance  occurring 
far  upstream  of  the  test  section.  An  example  of  such  a  disturbance  is  a  small  comer 
separation  in  the  adverse  pressure  gradient  of  the  tunnel  settling  chamber.  Another 
possible  cause  is  a  secondary  flow  in  the  corner  /"’ormed  by  the  wind  tunnel  side  wall  and 
floor,  dliis  is  an  example  of  Pra,ndtrs  secondary  flov.'  of  the  second  kind,  a-eated  by 
shearing  stresses,  a.s  described  by  Gessner  (1973).  Figures;  4.9U)-4.93b  reveal  the 
disturbances  due  to  this  outer  secondary  flow',  in  terms  of  incre.ased  5’  and  0,  and 
decreased  Q-  The  position  of  die  local  increase  in  b'  is  appro.xiniately  constant  at  Z/f 
=  -3.-5.  Tlie  position  of  maximum  K  and  0  moves  from  Z/f  =;  -1  ,,3  ;u  X/C  -•  3.00  to  2/V 
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=  -1.8  at  X/C  =  5.89,  where  it  remains  nearly  constant  with  increasing  X/C.  Figure  4.94 
displays  the  rapid  increase  in  A(5*./T)  and  A(0/T)  from  X/C  =  9.14  to  11.56,  which  is  due 
to  the  adverse  pressure  gradient.  Surprisingly,  the  data  do  not  show  a  substantial  increase 
in  slope  from  X/C  =  5.89  to  9.14.  Table  4.2  lists  the  maximum  changes  in  5*/T  and  0/T 
for  planes  13-15. 


lixpi'i'intenfal  Rc.sult.s 


65 


5.0  Data  Set  Comparisons 


This  chapter  focuses  on  the  primary  motivation  for  this  experimental  study:  the 
determination  of  the  effects  of  the  incoming  boundary  layer  and  flow  parameters  on  a 
turbulent  junction  vortex.  Tables  1.1  and  3.3  describe  the  data  sets,  and  Table  2.6  gives 
their  experimental  uncertainties.  For  discussion  purposes,  SCD  will  identify  Dickinson’s 
(19S5a,b)  data,  HMM  will  identify  McMahon  et  ai's  (1987)  data,  JS  will  designate 
Shin’s  (1989)  da'n,  JLF  will  denote  the  author’s  present  data,  WJD  will  identify 
Devenport  ai  d  Sin'pson's  (1990a)  LDV  d.-  a,  as  well  as  Devenpon  ei  al.'s  (1990)  hot¬ 
wire  data.  The  coacour  anu  velocity  vectc  piots  all  use  the  scale  1.0"  =  .5(T),  except 
plane  12  (X/C  =  3.00),  where  1.0"  ■-  1.0('.' 

The  symbols  used  to  designau  t*’  .  '^ous  data  sets  are  consistent  for  all  plots. 

Ihe  symbol  denoting  the  author’s  data  is  a  ^  re  tor  Devenport  and  Simpson’s  LDV 
data,  a  triangle,  for  Dickinson’s  data,  a  solid  i  .  McMahon  e[  al.'s  data,  a  diamond, 
and  for  Shin’s  data,  an  open  circle. 
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5.1  Examination  of  Data  Set  Blockage  Effects 

As  discussed  in  section  4,2,  Figure  4.4  shows  comparisons  between  2-D  potential 
flow  results  and  the  freestream  experimental  results  from  the  different  data  sets.  This  was 
done  to  obtain  a  greater  understanding  of  the  possible  blockage  effects  for  each  data  set 
to  avoid  any  possible  misinterpretions  of  the  differences  between  the  flows.  Shin’s  data 
was  corrected  to  account  for  the  variation  in  the  freestream  test  section  velocity  from  the 
reference  velocity  (UyEsr^ref  =  1.077). 

Dickinson’s  measurments  were  performed  using  a  24  inch  wide  tunnel,  with  a 
model  that  was  2.4  inches  thick,  creating  a  physical  blockage  of  10%.  No  allowances  for 
streamline  curvature  were  made  in  his  case.  Shin’s  data  was  obtained  in  a  36  inch  wide 
tunnel,  with  a  2.8  inch  thick  model,  creating  a  physical  blockage  of  7.8%.  Again,  no 
allowances  for  streamline  curvature  were  made.  Devenport’s  and  the  author’s  data  was 
taken  in  a  36  inch  v/ide  tunnel  with  a  2.4  inch  thick  model,  producing  a  physical  blockage 
of  6.7%.  However,  as  described  in  section  2.2.1,  side  wall  inserts  were  used  to  allow  for 
the  outer  streamline  curvature,  in  an  effort  to  create  negligible  blockage  effects.  Figure 
3.6  shows  that  the  blockage  effects,  if  any,  are  very  small. 

McMahon  et  al.'^  experimental  data  was  obtained  with  the  appendage  located  in 
a  free  jet  12  inches  downstream  of  the  tunnel  exit  plane.  This  configuration  is  not 
actually  free  of  interference  effects,  "^rhe  free  jet/ambient  air  interface  is  constrained  to 
follow  a  flow  streamline,  and  the  static  pre?  rre  along  this  line  must  'oe  equal  to  the 
ambieni  pressure.  However,  the  pressure  along  a  streamline  is  not  constant,  so  these  two 
interface  conditions  conflict,  leading  to  some  form  of  flow  interference.  No  oil-flow 
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visualizations  of  McMahon  et  al.'s  data  were  available,  but  the  mean  flow  contours  at 
X/C  =1,00  show  evidence  of  a  large  separated  flow  region  in  the  junction  comer,  which 
may  be  a  consequence  of  the  free  jet  test  flow  conditions, 

Figuie  4,4  shows  good  general  agreement  for  the  velocity  magnitudes  for  all  data 
sets,  especially  when  one  considers  that  the  jwtential  flow  calculation  did  not  account  for 
the  appendage  boundary  layer.  The  notable  exceptions  arc  in  plane  8,  where  McMahon’s 
frccstrcam  velocities  arc  considerably  lower  (5-10%)  than  the  calculatexl  results.  This 
result  is  surprising  considering  that  in  plane  10,  McMahon’s  freestream  velocities  deviate 
only  a  few  percent  from  the  2-D  potential  flow  results.  This  apparent  anomaly  may  i>e 
a  consequence  of  the  separated  flow  in  the  appendage  trailing  edge  region. 

Tnc  results  from  plane  9  show  another  unexpected  result.  Figure  4.4c  reveals  that 
die  author’s  freestream  velocities  deviate  more  from  die  calculated  results  than  those  of 
Dickinson’s  case,  which  had  a  higher  physical  blockage  and  also  did  not  allow  for 
streamline  curvature.  The  majority  of  the  data  shows  both  Dickinson’s  and  the  author’s 
data  to  agree  surprisingly  well,  with  Dickinson’s  ffecstream  velocities  usually  a  few 
percent  higher. 

The  effects  of  blockage  are  further  revealed  by  examing  the  freestream  skew 
angles.  Plane  5  data  shows  a  general  overestimation  of  otj-c  near  the  wing,  and  values  of 
Otjc  that  are  too  low'  away  from  the  wing.  This  latter  effect  is  easy  to  see,  as  blockage 
will  cause  th<'  flow  angle  to  tend  toward  zero.  Of  all  the  data,  Devenport’s  LDV  data  and 
Shin’s  data  appear  to  agree  the  best,  especially  near  the  maximum  values  of  Plane 
t  shows  good  Ore  agreement  for  both  the  author’s  and  Shin’s  data  sets,  within  1°  at 
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almost  all  locations.  Plane  7  shows  excellent  agreement  for  the  author’s  data. 
Dickinson’s  results  here  are  not  as  good,  showing  some  possible  blockage  effects.  Plane 
8  tells  the  same  story,  with  Devenport’s  and  McMahon’s  data  agreeing  very  well  with  the 
numerical  results,  but  D'  Idnson’s  data  is  constantly  offset  by  2°.  This  can  be  explained 
by  tunnel  blockage,  but  one  cannot  discount  a  bias  error  in  the  measurem.ents.  At  plane 
9,  the  author’s  data  shows  some  scatter,  but  shows  general  agreement  with  the  potential 
flow  results.  Dickinson’s  data  again  does  not  exhibit  freestream  flow  angles  as  large  as 
the  potential  flow  solution  predicts.  Skipping  to  plane  11,  the  Ofc  results  for  all 
concerned  data  sets  agree  well,  with  the  author’s  data  showing  some  scatter.  The  good 
agreement  here  is  most  likely  a  result  of  small  flow  angles  and  reduced  blockage  effects 
since  the  flow  has  passed  by  the  wing  at  this  location  (X/C  =  1.50). 

Considering  plane  10  (Figure  4.4f),  one  observes  two  things  right  away.  First,  the 
freestream  flow  angle  near  the  centerline,  as  predicted  by  potential  flow,  changes  radically 
from  X/C  =  .93  to  1.05.  Second,  Dickinson’s  results  are  not  even  close  to  comparing 
with  tlie  potential  flow  results,  or  with  the  results  of  the  other  data  sets,  'llie  other  data 
sets,  although  displaying  some  scatter,  follow  the  general  trend  set  by  the  numerical 
results  for  freestream  Ct^Q  quite  well.  Dickinson’s  data  shows  angles  that  differ  by  to 
8°.  This  discrepancy  is  difficult  to  explain  simply  by  blockage  effects,  and  section  5.3 
will  further  discuss  this  problem. 
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5.2  Expected  Effects  of  the  Approach  Boundary  Layer  Properties 

Before  comparing  data  sets  to  examine  the  effects  of  varying  the  approach 
boundary  properties,  it  is  helpful  to  use  intuition  to  gain  any  possible  a  priori  knowledge, 
so  that  one  knows  what  to  look  for  and  what  differences  to  expect.  This  section  discusses 
what  effects  different  approach  boundary  layers  may  have  on  the  wing-body  flow 
structure. 

From  the  previous  discussion  of  Baker’s  (1980,  1985)  and  Belik’s  (1973)  work, 
we  sec  that  the  initial  vortex  structure  in  front  of  a  cylinder  (equivalent  to  the  nose  of  an 
appendage)  is  primarily  a  function  of  two  parameters,  Re^  and  0/T.  The  displacement 
thicknes.s  6*  is  interchangeable  with  ti.  Reg,  (or  RCg )  effects  the  positions  of  Xy  and 
Xj,  and  the  basic  vortex  structure  seems  to  be  dependent  on  5*/D  (or  SYf)  for  a  given 
range  of  Rcg*. 

5.2.1  Expected  Re^  Effects 

Baker  (198U,  1985)  found  that  Xy  increases  as  Reg*  .  vh.ich  means  that 

the  vortex’s  horizontal  size  is  increasing.  He  conjectured  that  changing  Reg,  does  not 
change  the  vertical  size  of  the  primary  vortex.  For  Rcg,  <  10,000,  Xg  increases  also,  but 
at  a  slower  rate  tlian  Xy,  which  means  that  tire  distance  between  Xy  and  Xg  is  decreasing 
as  Reg*  increases.  This  distance  decreases  even  faster  when  Reg,  >  10,000,  as  Xj  starts 
to  decrease  at  the  larger  Reynolds  numbers.  One  can  conclude  then,  that  as  Reg, 
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increases,  the  main  vortex  will  increase  in  size  (Xy  increases),  and  will  have  a  more 


energetic  initial  roll-up  (Xg  -  Xy  decTeases). 

Ignoring  the  large  scale  unsteadiness  present  in  the  junction  flow,  a  rather  simple, 
intuitive  argument  would  go  as  follows.  Rce  or  Res*,  for  the  flows  considered  in  this 
study,  is  essentially  a  function  of  (or  U^jS*),  since  v  is  fairly  constant  for  a  given 
fluid  medium  (air  or  water).  The  freestream  fluid  momentum  is  direcdy  proportional  to 
Ufef  and  the  lack  of  high  momentum  fluid  near  the  wall  is  indicated  by  9.  Therefore,  the 
product  U„j6  is  indicative  of  tlie  absolute  relative  difference  in  fluid  momentum  between 
the  boundary  layer  and  the  freestream.  If  the  momentum  difference  is  large  (large  Ree), 
then  the  incoming  high  momentum  fluid  easily  circulates  down  into  the  boundary  layer 
and  fonns  a  stfungly  denned  backfiow  region,  if  the  momentum  difference  is  small 
(small  RCg),  the  backflow  region  is  more  diffuse.  Realizing  that  Urci9  indicates  the 
typical  magnitudes  of  the  inertial  forces  in  the  boundary  layer,  yields  the  same  hypothesis. 

Studies  done  in  a  water  tunnel  within  the  Virginia  Tech  Aerospace  Department  by 
Kim  (1991),  ■'port  these  speculations.  As  previously  mentioned,  Kim  performed 
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cxpcriiuerits  at  2  diffeicni  values  of  Rcq,  330  and  i,iOu,  on  the  same  body  shape  as  the 


author’s.  Mean  velocity  profiles  along  ihe  centerline  for  the  2  test  cases  display  a  thinner 
and  stronger  backflow  region  the  higher  Reynolds  number  case  (Figure  5.1).  The  higher 
Reynolds  number  case  displays  the  enhanced  secondary  flow  "jetting"  effect  due  to  the 
relative  lack  of  momentum  near  the  wall.  Figure  5.2  presents  contours  of  turbulence 
intensity,  showing  more  distinct  local  maxima  concentrated  near  the  wall  for  Rcq  =  1,100. 


The  lower  Reynolds  number  case  shows  several  turbulent  "pockets",  suggesting  a  more 
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complex  mean  flow  pattern.  The  turbulence  is  spread  over  a  much  wider  area  for  the 
Ret,  =  330  case,  displaying  a  more  diffused  secondary  flow  structure. 

Realizing  that  U„f0  is  not  actually  the  relative  momentum  difference  between  the 
outer  and  inner  boundary  layer  regions,  the  author  has  come  up  with  a  new  parameter, 
referred  to  as  the  momentum  deficit  factor  (MDF).  The  momentum  flux  through  a 
boundary  layer  is  (p)(U„f)‘(6-0)  per  unit  width.  The  "momentum  deficit"  is  equal  to 
pOdref)^(Q)  P®r  ttttit  widtli.  If  T  is  used  as  a  characteristic  widtli,  then  a  form  of  MDF  is 
P(U,ef)^(6)(T).  This  is  dimensional,  however,  so  the  final,  non-dimensional  form  for  MDF 
is  (ReT)^(0/T)  which  is  also  equal  to  (Rep)(ReQ).  Flows  with  large  values  of  MDF  will 
tend  toward  smaller  flow  distortions  and  tlie  vortical  and  turbulent  flow  structures  will  be 
more  intense  near  the  wall,  since  the  low  momentum  fluid  region  acts  as  a  type  of  "relief 
valve".  Table  3.3  shows  that  for  the  flows  examined  in  tiiis  study,  Rcq  usually  varies  as 
MDF.  The  exception  to  this  trend  is  at  the  lower  Re^  cases,  where  McMahon  et  al.'s 
data  has  the  lowest  MDF,  despite  having  a  higher  Reg  than  either  the  author’s  or 
Devenport  et  al.'s  thin  boundary  layer  case. 

The  downstream  effects  on  vortex  position  and  re.sulting  flow  distortion  due  to  the 
resistance  of  the  inner  boundary  layer  to  the  intrusion  of  freestieara  fluid  can  be  predicted 
by  considering  Shizawa  and  Eaton's  (1990)  study.  As  previously  discussed  in  .section  4.3, 
two  cases  of  a  voilex  embedded  in  a  skewed  boundary  layer  were  studied.  Case  1,  where 
the  near  wall  .skewing-induced  velocity  was  in  the  same  direction  as  die  vortex-induced 
velocity,  corresponds  to  a  case  with  ver)'  low  resistance  to  the  iatnision  of  outer  fluid  into 
the  boundary  layer  (high  MDF).  Case  2,  where  the  near  wall  skewing-  and  vonex- 
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induced  velocities  pposed  one  another,  corresponds  to  a  case  with  high  resistance  to  the 
intnision  of  outer  fluid  into  the  boundary'  layer  (low  MDFj.  Sliizawa  and  Eaton  found 
that  the  vortex  associated  distortion  patterns  for  case  I  were  not  as  strong  and  decayed 
more  quickly  than  case  2.  This  was  because  the  primary  tlow  distonlnn  mechanism  (the 
upward  transport  of  low  momentum  fluid  due  to  separation)  was  enhanced  in  case  2.  The 
vortex  core  position  moved  further  away  from  the  wall  for  case  2.  More  information 
concerning  the  vorticity  distributions  for  the  two  cases  is  discussed  in  section  4.3. 

S.2J  Expected  5* /T  Effects^ 

As  mentioned  previously,  the  ratio  dVf  (or  6/T)  may  effect  the  initial  voitex 
structure.  Based  on  the  research  of  previous  workers  (Baker  (1980,  1985),  LaRsur  and 
Langston  (1990),  and  Eckerle  and  Langston  (1987)),  it  can  be  postulated  that  the  number 
of  vorticies  along  the  centerline  increases  as  5*/T  increases.  An  argument  supporting  this 
hypothesis  is  that  as  5*  increases,  the  corresponding  decrease  in  near  wall  momentum 
effects  the  primary  separation  mecnanism,  causing  a  larger  numtier  of  discrete  vorticies 
to  be  fonned  in  the  plane  of  symmetry.  Using  flow  topology,  researchers  have  suggested 
several  initial  vortex  coniigurations  (2,  4,  6  vonex  models),  and  for  certain  ra.iges  of  5*yT 
values,  the  number  of  vorticies  seems  to  be  constant.  A  sketch  showing  the  various 
vortex  configurations  is  given  in  Figure  5.3. 

However,  the  relevance  of  the  initial  vortex  structure  on  the  dov.  nstreain  flow  is 
moot  issue,  based  on  the  conclusions  drawn  from  the  surface  oil-flow  visualizations 
(section  3.1.3).  -Any  insiantaueous,  discrete  vortex  structure  present  in  the  nose  region 
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appears  to  fonn  a  tinie- averaged  primary'  vortex  structure  downstream  of  the  body’s 
maximum  thickness  location.  The  primaiy  effect  of  5'‘n  may  be  its  influence  on  the 
bimodal  f  ow  structure  at  the  appetidage  nose,  which  would  in  turn  affect  the  downstream 
turbulence  characteristics. 

The  ratio  5‘/T  is  also  an  indication  of  the  mass  entrainment  potential  of  the  vortex. 
Since  previous  measurments  (Devenport  e:  aL,  1990)  show  that  titc  primaiy  vortex 
dimensions  scale  on  T,  this  parameter  relates  the  relative  size  of  the  boundary  layer  to  the 
horseshoe  vortex.  Smaller  values  of  6’/T  will  enable  the  vonical  flow  to  have  a  greater 
effect  on  the  entire  boundary  layer,  and  on  tlie  availability  of  high  momentum  fluid  which 
energizes  the  secondary  flow. 

5,2.3  Significance  of  V„JIT  and  t/,,y5 

U,yr  and  Ij,^6  form  two  additional  parameters  of  interest.  These  ratios  are 
dimensional  (units  of  sec'*),  and  as  such,  they  are  not  true  flow  parameters  by  themicives. 
Ihese  parameters  are  useful,  however,  for  non-dimensionalization  of  vorticity,  and 
understanding  their  effects  gives  some  insight  into  the  basic  flow'  physics. 

The  average  voitex  stretching  rate  for  a  particular  appendage  geometry  is 
controlled  by  U„f/T.  According  to  potential  flow  theory,  as  the  vortex  stretching  rate 
increases,  the  peak  dimensional  vorticity  should  increase,  due  to  tlie  shrinking  of  the 
vortex  tube.  The  parameter  U„./5  is  related  to  the  average  spanwise  vorticity  in  the 
incoming  boundary  layer.  Therefore,  increasing  U„(/5  should  ulso  increase  the 
dimensionalized  .strength  of  the  horseshoe  vortex,  but  not  the  non-din  ensioniilized 
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strength.  Determining  the  which  of  these  two  parameters  are  most  appropriate  for  non- 
dimensionalizing  the  vorticity  increases  understanding  of  the  junction  vortex  flow. 

5,2.4  Parametric  Space  of  Data  Sets 

Figure  5.4  clearly  displays  the  varying  parameiers  of  each  of  the  data  sets 
examined,  in  this  study.  Shin’s  parameters,  when  measured  under  conditions  compai-able 
to  the  other  data  sets,  may  fail  within  the  region  enclosed  by  die  dashed  line  in  ihe  upper 
right  corner.  This  is  further  oiscussed  at  the  beginning  of  section  5.3. 

Tlie  lines  of  constant  MDF  in  Figure  5.4  show  its  variation  with  changing  Ree  and 
G/T.  Tlie  influences  of  Rcg  and  0/T  on  the  juncture  flow  can  be  easily  compared  to  the 
cifecrs  of  MT)F  using  this  figure.  Figure  5.4  also  can  be  used  as  an  atd  m  the  design  of 
future  parametric  studies  of  wing- body  junction  flows. 


5.3  Data  Set  Comparison  Resnlts 

The  contour  and  secondary  velocity  plots  will  be  the  main  tools  for  comparing  the 
data.  The  profile  plot  conipari.sons  provide  much  more  detaii,  and  in  some  instances  this 
will  help  to  sec  the  effects  of  the  approach  bounda,ry  layer.  One  must  be  cautious, 
however,  tiiat  me  experimental  uncertainties  ao  not  become  a  factor  when  comparing 
subtle  details  between  the  flows. 
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One  fact  is  readily  apparent  when  comparing  the  data  sets.  The  secondary  flow 
structure  scales  on  T,  the  maximum  appendage  thickness,  in  both  the  Y  and  Z  directions, 
as  can  be  seen  by  the  profile  plots  of  (Figures  4.32-4.38),  the  contour  plots  of 

turbulence  intensity  (Figures  5.40-5.46),  and  secvJidary  velocity  vector  plots  (Figures 
5.20-5.26).  The  profile  plots  of  fp  show  that  the  distortion  function  also  scales  on  T.  The 
data  sets  do  not  compare  precisely,  but  one  can  easily  observe  that  T  is  more  appropriate 
as  a  scaling  factor  than  5.  The  importance  of  T  as  a  scaling  factor  reflects  the 
dependance  of  the  secondary  flow  on  the  pressure  field  created  by  the  appendage 
geometry.  All  significant  flowfield  characteristics  originate  from  the  geometry  of  the 
wing-body  junction. 

A  few  things  to  note  before  continuing  i.'e  listed  as  follows: 

■'  The  momentum  deficit  factor  for  Shin’s  data  may  be  effectively  larger  than  the 
listed  value,  even  larger  than  Dickinson’s  value.  This  is  because  the  actual  test  section 
velocity  is  higher  than  the  given  reference  velocity  (by  a  factor  of  1.077),  and  also  the 
given  value  of  6  was  determined  with  no  appendage  in  the  tunnel.  Correcting  MDF  for 
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create  a  value  of  MDF  equal  to  Dickinson’s. 

•  Shin’s  secondary  velocities  near  tne  wall  usually  show  an  extremely  large  V 
component,  which  is  not  seen  in  flows  with  similiar  boundary  layer  parameters  (the  SCD 
data  set).  The  author  believes  this  apparent  anomaly  was  caused  by  the  probe.  This  was 
the  only  data  obtained  using  a  5  hole  pressure  probe,  and  the  relatively  large  size  of  this 
probe  could  have  interfered  with  the  flow  and/or  caused  bias  errors  to  occur  near  the  wall. 


m 
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The  probe  was  3.18  mm  in  diameter,  and  measurements  were  taken  to  within  3.7  mm  of 
the  wall. 

•  Dickinson’s  secondary  velocity  and  vorticiiy  distribution  data  appear  to  be  in  error 
in  plane  10  (Figures  5.25  and  5.35).  The  flow  shown  by  his  data  is  markedly  different 
from  the  other  data  sets,  and  no  reason  can  be  found  to  explain  it.  The  data  in  planes  9 
and  11  compare  well  with  otlicr  data  sets,  which  leads  the  author  to  believe  that  the 
unusual  secondary  flow  pattern  in  plane  10  is  due  to  an  error.  The  streamwise  mean  and 
fluctuating  velocities  in  plane  10  do  not  appear  to  be  significantly  effected  by  this 
problem. 

•  A  relatively  large  separated  flow  region  is  apparent  in  McMahon  et  a/.’s  flow  at 

the  trailing  edge.  This  can  be  seen  in  contours  of  in  plane  10  (Figure  5.14). 

5.3,1  Mean  Velocity  Distribution 

Contour  plots  of  UTgAJref  and  the  associated  contours  of  distortion  function  (f^) 
are  shown  in  Figures  5.5-5.19.  Secondary  velocity  plots  are  given  in  Figures  5.20-5.26. 
These  figures  reveal  the  distribution  of  the  mean  velocity  components  within  the  junction 


flow.  When  comparing  "distortion  levels",  it  is  the  plums  of  low  miomientum  fluid  near 
the  line  of  low  shear  location  that  is  being  referred  to.  Greater  flow  distortion  is  indicated 
by  lower  values  of  f^. 

In  general,  the  SCD  and  JS  data  sets  show  broader,  more  diffused  distortion 
patterns  for  the  contours  of  Ujcr  while  contours  of  fo  reveal  that  the  HMM  data  has  the 
highest  distortion  rates,  with  the  JLFAVJD  data  falling  in  between.  This  is  the  trend 
predicted  by  the  variation  in  MDF.  Using  Rcg  as  a  predictor,  the  JLF/'A^JD  data  would 
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have  the  highest  distortion,  followed  by  the  HMM  data,  and  then  the  SCD  and  JS  data 
sets. 

The  validity  of  the  concept  behind  the  momentum  deficit  factor  is  perhaps  best 
shown  in  plane  12,  the  farthest  downstream  location  where  comparisc 's  vere  made. 
Figures  5.18  and  5.19  show  the  contours  of  U^c  ^ind  f[j.  The  WJD  thin  boundary  layer 
case,  has  the  lowest  Reg  (4,500  compared  to  6,300  and  6,800  for  the  JLF/WJD  and  HMM 
data  sets  respectively).  However,  the  distortion  function  results  clearly  show  that  the 
HMM  data  has  the  largest  distortion  levels  (lowest  values  of  y,  while  the  distortion 
pattern  in  the  "tongue"  of  low  momentum  fluid  is  approximately  the  same  for  the  WJD 
thick  and  thin  boundary  layer  cases.  These  latter  two  cases  differ  in  MDF  by  1 1%,  while 
the  WJD  thin  boundary  layer  case  has  an  MDF  which  is  42%  greater  than  the  HMM  case. 

From  the  examinarion  of  the  WJD  thin  and  thick  boundary  layer  data  in  plane  12, 
and  also  the  HMM  and  WJD  data  in  plane  8  (Figure  5.10),  one  can  see  that  f^  does  scale 
on  T  in  the  Y  direction.  One  can  conclude  that  if  die  HMM  and  WJD  data  had  more 
siniiliar  values  of  MDF,  then  the  fp  contours  would  look  more  alike. 

Plane  by  plane  comparisons  consistendy  yield  the  same  results.  Tlie  distorted 
contours  of  Utc  ‘ire  more  peaked  for  the  cases  with  low  values  of  MDF,  and  are  broader 
for  the  cases  with  high  MDF.  The  contours  of  f^  show  this  trend  in  even  greater  detail. 

The  effects  of  boundary  layer  thicknesses  to  body  thickness  ratios  are  seen 
primarily  in  the  region  between  the  body  and  ‘he  LOLS,  where  tlie  local  boundary  layer 
is  thinned  by  the  sweeping  of  high  momentum  fluid  into  the  junction  by  the  horseshoe 
vortex.  The  increase  in  relative  size  of  the  horseshoe  vortex  compared  to  the  boundary 
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layer  for  flows  with  low  values  of  6/T  results  in  higher  momentum  fluid  being  available 
for  entrainment  into  the  junction.  Figures  5.6-5.19  show  slightly  larger  values  of  (  > 

I. 0)  inboard  of  the  horseshoe  vortex  for  flows  with  low  values  of  b/T.  This  is  a  more 
subtle  effect  than  the  distortions  due  to  the  uplifting  of  low  momentum  fluid  away  from 
the  wall,  and  in  some  instan'.:es,  such  as  plane  1  i,  the  contours  of  fp  in  the  thinned  region 
look  almost  identical  for  diffetent  va'ues  of  5/1. 

As  the  flow  moves  downstream,  another  consequence  of  the  momentum  deficit 
factor  is  seen.  The  distortion  patterns  of  tire  flows  with  higher  MDF  are  centered  further 
outboard,  at  larger  values  of  |Z/T|.  This  is  an  expected  trend,  and  can  be  explained  as 
follows.  For  flows  with  higher  values  of  MDF,  tlie  flows  are  similiar  to  Shirawa  and 
Eaton’s  (1990)  case  1  flows  (see  sections  5.2  and  4.3).  The  vortex-induced  velocities 
between  die  wall  and  the  vortex  core  are  larger  for  case  1  than  .Vo  ■  case  2,  which  is  more 
like  a  flow  with  high  MDF.  This  condition  leads  to  larger  vorti-ciiy  values  for  the  high 
MDF  cases  (Shizawa  and  Eaton’s  case  1  data  show  this),  and  the  increased  voitex- 
induced  velocities  move  the  horseshoe  voitex  further  from  the  centerline.  The  observed 
similarities  to  Shizawa  and  Eaton’s  flows  suggests  that  changing  MDF  results  in  a  change 
in  the  effective  gsometry-inducefl  skewing  magnitudes. 

In  planes  10  and  1 1,  die  secondary'  velocity  vector  plots  show  the  hoi’sesKoe  vonex 
structure  centered  farther  from  the  centerline  for  flows  with  larger  MDF  values.  In  plane 

II,  it  appears  that  the  author’s  earlier  hypodicsis  that  MDF  for  the  JS  data  may  actually 
be  slighUy  higher  than  for  the  SCD  data  is  correct,  since  Shin’s  data  shows  the  most 
outboard  vortex  location.  Surprisingly,  the  spanwise  positions  for  the  primary  distonion 
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do  not  vary  much  between  the  data  sets  in  plane  12.  The  position  for  the  HMM  data 
seems  farthest  from  the  centerline,  contrary  to  the  previous  discussion.  This  contradiction 
may  be  due  to  the  trailing  edge  separation  present  in  HMM’s  case. 

The  HMM  secondary  velocity  results  in  plane  10  reveal  a  quite  circular  secondary 
flow  pattern  near  the  centerline,  unlike  that  observed  for  the  flows  with  higher  MDF 
values.  The  proximity  of  this  flow  structure  to  the  appendage  trailing  edge  may  be  the 
cause  of  the  significant  flow  separation  for  the  HMM  data  set.  The  other  flows  all  have 
a  rather  flat,  elliptically  shaped  pattern,  with  no  substantial  V  velocities  except  perhaps 
near  the  LOLS  location  or  the  body.  These  observations  substantiate  the  predicted  "relief 
valve"  effect  provided  by  the  inner  boundary  layer  for  flows  with  high  MDF.  The  near 
wall  channelling  of  this  fluid  marked  by  positive  vorticity  is  discussed  in  the  next  section. 

The  profile  plots  of  mean  velocities  reveal  additional  information  concerning  the 
con(,;'..t  scaling  of  uie  mean  flow  distribution.  ITie  "kinking"  of  the  UTcAlrcf  profiles 
(Figures  -I-.  10-*^. 24)  tends  to  scale  on  T,  although  the  magnitudes  are  generally  not  the 
some.  The  distordon  function  filters  out  the  effects  of  the  "background"  flow. 

The  apj/ropra  e  scaling  for  the  distribution  of  is  more  difficult  to 

determine.  Some  profile  s  disjiay  an  apparent  5  scaling  (Tigure  4.29),  and  otheis  seem 
to  scale  on  T  (Figure  4.30).  Hows  with  lower  MDF  show  higher  magnitudes,  which 
is  expected  due  to  the  observed  higher  distortion  levels  for  these  flows. 

Most  profiles  show  a  very  strong  dependence  on  T  (Figures  4.32- 

4.38).  The  variation  in  ffecstream  flow  angles  makes  it  difficult  to  compare  the  relative 
skewing  magnitudes  for  the  data  sets.  The  hodograph  plots  (Figures  4.77-4.83)  reveal 
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very  sharply  defined  crossover  profiles  for  data  sets  with  low  5/T  ratios.  The  skewing 
magnitudes  in  general  appear  to  be  slightly  greater  for  these  flows  also.  This  is  indicative 
of  the  availability  of  higher  momentum  fluid  closer  to  the  horseshoe  vortex  secondary 
flow  structure.  When  this  fluid  is  swept  into  the  boundary  layer,  it  causes  larger  skewing 
angles  and  steeper  gradients  in  Wxc- 

5.3.2  Vorticity  and  Helicity  Distribution 

From  observing  the  contours  of  (0(W)/3Y)5AJref  helicity  (defined  by  equation 
(4.3)),  the  concentration  of  voiticity  near  the  wall  is  immediately  evident  for  flows  with 
higher  values  of  N'lDF  (Figures  5,27-5.38).  JS  vorticity  levels  are  shown  in  the  individual 
profile  plots  of  (0Cw)/0Y)6/u\^c  (Tigures  4.84-4.90),  and  are  most  similiar  to  the  SCD 
distribution  patterns,  reaching  higher  values  near  the  wall  in  some  cases.  This  effect  is 
a  result  of  the  relative  lack  of  resistance  to  the  instnjsion  of  the  secondary  flow,  creating 
a  channelling  effect  near  the  wall.  The  profile  plots  of  vorticity  also  indicate  that  the 
peaks  in  voiticity  scale  more  appropriately  on  T  than  on  6.  The  no-slip  condition  at  the 
wall  was  not  used  for  caicuiating  (d(W)/d  Y)  near  the  wall,  since  this  technique  can  lead 
to  significant  errors  for  estimating  derivatives. 

The  trend  in  the  distribution  of  positive  and  negative  vorticity  between  the  flows 
with  high  and  low  MDF  are  again  similiar  to  the  flows  studied  by  Shizawa  and  Eaton 
(1990).  The  lower  MDF  flows  (liMM,  JLFAVTD),  previously  showing  characteristics  of 
Shizawa  and  Eaton’s  case  2  flow,  show  a  relatively  stronger  pocket  of  negative  vorticity 
outboard  of  the  primary  vonex  (Figure  5.31).  A  stronger  negative  vorticity  region  is  also 
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observed  for  Shizawa  and  Eaton’s  case  2  flow,  when  compared  to  case  1.  The  strengths 
of  the  negative  vorticity  legions  are  reflected  by  the  YfT  locations  of  the  vortex  cores. 
The  variation  in  MDF  between  the  data  sets  correctly  predicts  the  relative  Y/T  positions 
of  the  vortex  cores  (as  MDE  increases,  the  core  is  closer  to  the  wall). 


By  comparing  the  JS  and  SCD  data  to  the  JLFAVJD  data  sets,  the  correct  vorticity 
scaling  parameter  appears  to  be  T/U^^f  (remember  that  (T/5)jlf  =  2,  and  (T/5)sc[jj5  ==  1). 
This  may  be  because  the  average  vortex  stretching  rate,  which  is  proportional  to  U„/T, 
is  a  dominant  factor  affecting  vorticity  magnitudes  around  the  appendage  nose. 
Comparisons  between  the  HMM  and  WJB  data  in  plane  8  (Figure  5.31),  and  also  between 
the  JS  and  SCD  data  sets  (Figures  4.84-4.90),  reveal  that  5/U„f  is  the  appropriate  scaling 
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10  (Figure  5.35) 


shows  similiar  results,  although  the  HMM  contours  show  no  well  defined  "core".  The 


regions  of  high  positive  vorticity  near  the  wing-body  corner  seen  for  the  HMM  data  is 
likely  due  to  the  separation  observed  in  tliat  region. 

The  determination  of  which  flows  scale  on  the  same  parameters  may  be  an 
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(3(W)/3Y)5/U,ef  seen  in  the  SCD  and  JS  data  are  app.roximately  twice  as  large  as  for  the 


HMM  or  JLF,YVJD  data  sets.  The  ratios  of  5/T  for  the  JS  and  SCD  cases  are  twice  as 


large  as  for  JLF/WJD,  but  a^c  approximately  the  same  as  for  HMM’s  case,  indicating  that 
neither  T/'Uj.^f  or  8Al,ef  would  scale  the  SCD/JS  and  HMM  flows  to  the  same  levels.  Note 
that  Rcg  for  the  JS  and  SCD  data  is  twice  as  large  as  JLFAVJO’s  and  HMM’s  values. 
Ob.serving  that  the  vorticity  levels  for  the  SCD  and  JS  data  compare  when  scaled  on 
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and  that  the  vorticity  levels  for  die  HMM  and  JLFAV'JD  data  are  also 
approximately  scaled  by  SAir^f,  one  niay  postulate  that  the  vorticity  levels,  non- 
dimensionalized  by  5/Uref  arc  dependent  on  Rcg  or  MDF.  The  contours  of  helicity 
support  these  observations,  indicating  that  flows  with  high  values  of  MDF  have  more 
coherent  large  scale  structures. 

The  decay  of  vorticity  in  the  horseshoe  vortex  core  vs.  X/C  is  shown  in  Figure 
5.39.  Note  that  the  maximum  core  values  of  0(W)/aY)6AJ,cf  at  X/C  =  1.05  and  1.5  for 
the  JLF/WJD  and  HMM  data  sets  were  estimated  from  the  contour  plots.  The  semi-log 
plots  show  that  a  power-law  relationship  may  describe  the  decay  of  voriticity.  The  decay- 
rate  increases  as  MDF  increases,  but  this  trend  may  be  the  result  of  the  scaling  effects  of 
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with  increasing  X/C,  suggesting  a  growing  dependence  on  instead  of  T/U„f. 


5,3.3  Distribution  of  Turbulent  Quantities 

In  general,  the  five  turbulent  quantities  measured  agree  quite  well  between  the  data 
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show  that  the  turbulence  structure  is  largely  governed  by  the  appendage  geometry, 
therefore  scaling  on  T  in  the  Y  and  Z  directions. 

Kim’s  (1991)  contours  of  taken  in  the  plane  of  symmetry  at  the 

appendage  nose,  show  diat  increasing  Reg  confines  and  concentrates  the  turbulence 
stmeture  in  a  near-wall  region  (Figure  5.2).  Tliis  result  seems  to  indicate  that  the  higher 
RCg  flow  (SCD)  would  show  a  higher  concentration  of  turbulent  stnictures  near  the  wall, 
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but  this  is  not  necessarily  true.  There  is  no  obvious  concentration  of  turbulent  n.uid  near 
the  wall,  but  in  planes  8-1 1,  the  contour  plots  show  slightly  larger  values  of  wVc/Urcf  near 
the  wall  for  the  SCD  data.  The  profile  plots  of  u'jc  reveal  how  well  T  scales  the  flow. 

The  JUF/WjTD  contours  of  u\q  shows  more  distinct  "islands"  of  local  turbulence 
peaks.  This  raay  be  an  effect  of  the  ratio  6/T,  since  a  thinner  boundary  layer  has  an 
initially  thinner  region  marked  by  highly  turbulent  fluid,  resulting  in  a  less  diffused 
turbulence  structure  in  the  wing-body  junction.  The  HMM  contour  s  show  neither  local 
"islands"  of  pealc  turbulence  intensity  nor  relatively  high  regions  of  turbulence  near  the 
waU.  This  is  indicative  of  the  low  Rcg,  high  d/T  approach  flow  parameters.  Note  the 
large  separation  region  at  the  trailing  edge  for  the  HMM  flow,  revealed  by  the  high 
turbulence  levels  as  Z/T  goes  to  zero  (Figure  5.44). 

By  plane  12,  at  X/C  =  3.00,  the  flow  does  not  scale  on  T  quite  as  well  as  it  does 
near  the  appendage  (Figure  5.46).  The  vertical  growth  of  the  turbulent  structures  as  X/C 
increases  appears  to  be  limited  by  the  boundary  layer  thickness,  since  they  cannot  grow 
as  quickly  into  the  freestream.  When  one  considers  that  5/T  roughly  doubles  for  each 
contour  plot  in  Figure  5.46,  the  scaling  effects  of  T  arc  qulLC  remarkable. 

Data  for  v  was  compared  in  planes  5,  8,  and  10  (Figures  4.57-4.59).  The 
data  from  plane  5  clearly  shows  tlie  scaling  effects  of  T.  The  cause  of  the  larger  peaks 
in  v'tc  ttt  this  location  for  the  WJD  data  is  not  obvious,  as  plane  8  shows  sirniliar 
turbulence  levels  for  the  SCD  and  WJD  data  sets.  The  initial  separation  for  the  lower 
Rcg  (or  MDF)  flow  may  be  "snonger"  in  a  time-averaged  sense,  causing  higher  rms  V 
velocity  fluctuations.  The  HMM  data  in  phmes  8  and  10  show  significantly  higher  levels 
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of  v't-c-  This  appears  to  be  a  combination  of  low  RCg,  high  5/T  effects,  since  the  SCD 
and  WJD  data  sets  are  comparable.  Judging  from  profiles  of  V-rtAJrcf  in  plane  10  (Figure 
4.31),  high  V  flnctuations  may  be  a  direct  result  of  large  mean  V  velocities,  caused  by 
the  combination  of  low  Rcg,  high  5/T. 

The  spanwise  velocity  fluctuations  are  different  from  the  others,  in  that  the 
distribudon  of  seems  to  scale  on  5  instead  of  T.  This  may  be  because  does 

not  show  much  variation  through  the  boundary'  layer,  and  in  most  instances  the  only 
identifiable  change  "^^curs  when  the  turbulence  levels  drop  to  their  frecstream  values.  The 
location  of  this  drop,  of  course,  scales  on  5.  In  planes  8  and  10  (Figures  4.63  and  4.65), 
the  w'  levels  from  die  HMM  data  set  are  slighdy  lower  than  the  other  2  data  sets  (SCD 
and  JLF/WJD).  This  again  appears  to  be  a  low  Rcg,  high  5/T  effect,  since  the  SCD  and 
WJD  data  sets  agree  very  well  in  the  inner  boundary  layer.  It  is  perhaps  linked  to  the 
generally  lower  spanwise  mean  flow  velocities  observed  in  the  HMM  flow,  particularly 
evident  in  plane  8  (Figure  4.63). 

In  a  fashion  similiar  to  the  distribudon  of  the  mean  Reynolds  stress 

uviq/(IJj.c)^  appears  to  generally  scale  on  5,  not  T  (Figures  4.61  and  4.69).  This  is  .seen 
by  comparing  the  local  maxima  and  minima,  and  nodng  that  the  SCD  and  HMM  data 
compare  better  to  each  ether  than  to  the  WJD  data.  Tlie  most  negative  values  of  tre 
usually  seen  for  the  HMM  data  set.  It  is  difficult  to  draw  any  definite  conclusions 
regarding  iTT'tc-  due  to  the  scatter  in  the  data  and  the  lack  of  data  for  comparison. 

The  anal  turbulence  quantity  compared  is  iIvVj(;^/(U,pf)^.  Unlike  uv-y^,  the 
distribution  of  clearly  scales  on  T  (Figures  4.70-4.76).  Except  for  the  sepiu’ated 
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regions  present  in  the  HMM  flow,  the  magnitudes  of  ususally  agree  quite  closely 
between  the  data  sets.  ITiis  similiarity  shows  tliat  iwjc  is  controlled  by  the  horseshoe 
vortex  flow  structure,  which  is  driven  by  the  junction  geometry.  This  is  contrary  to  mvtc> 
which  appears  to  depend  more  on  the  initial  boundary  layer  thickness. 


5.3.4  Boundary  Layer  Thickness  and  Cj-  Distributions 

Figures  4.91  and  4.92  show  results  of  575*app  ^/^app  for  all  the  data 

planes  (5-15).  These  results  are  in  the  ffeestieam  coordinate  system.  The  substantial 
differences  shown  by  JS’s  data  in  plane  5  (5*)  is  an  indication  of  the  approach  boundary 
layer  properties  being  determined  with  no  appendage  in  the  tunnel.  The  5*  and  9  data, 
in  general,  shows  that  the  data  sets  with  lower  MDF  values  have  a  larger  peaks  near  the 
T  OLS  locations,  while  the  data  sets  with  higher  MDF  values  tend  not  to  have  peaks  in 
boundary  layer  thicknesses.  The  ratio  SAT  appears  to  play  some  importance  here, 
since  the  JLI’AVJD  data,  while  similiar  to  the  HMM  data,  has  steeper  gradients  of  5*  in 
the  spanwise  direction.  The  larger  values  of  5*  and  6  near  the  body  for  HMM’s  data  is 
due  to  the  flow  separation  for  that  case.  The  JS  and  SCD  data  sets  show  no  local 
maxima  in  the  majority  of  the  measurement  planes,  but  5*  and  0  do  decrease  as  |Z/r| 
decreases  toward  the  wing. 

Surprisingly,  the  distortive  effects  when  scaled  on  T  (A(5*/T)  and  A(0,''r))  compare 
very  closely.  The  values  for  the  HMM  data  are  highest,  but  are  lowest  for  the  JLF/WJD 
data,  so  nothing  definite  can  be  said  relating  MDF  or  Ree  to  the  variation  in  A(6‘A’)  or 
A(0/T)  vs.  X/C. 
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Results  of  (Cj)ps  ■vs.  Z/T  are  show-n  in  Figure  4.93  for  all  data  sets.  The  skin- 
friction  coefficient  in  the  freestream  direction  was  estimated  using  equation  (4.10).  As 
expected,  the  data  sets  with  the  largest  boundary  layer  thicknesses  had  the  smallest  Cf 
values.  These  data  sets  also  had  a  much  smaller  spatiwise  variation  in  Cf.  A  local 
minima  in  Cf  was  usually  located  jus:  out.side  the  LOLS  as  deteimined  by  the  oil-flows, 
and  Q  increased  for  ail  data  sets  once  inside  the  LOLS  location.  Plane  5  is  an  exception, 
as  no  local  minima  was  found  in  Cf  at  this  X/C  location,  Cf  decreased  monotomically  as 
[Z/Tj  increased.  The  LOLS,  however,  is  not  well  defined  at  the  maximum  thickness, 
so  this  result  is  not  unexpected.. 

S.3.5  Comparison  of  Hot-Wire  and  LDV  Results 

Overall,  diese  results  agree  well.  The  profile  plots  comparing  the  JLF  data  to  the 
WJD  1X)V  data  in  planes  5  and  10  show  the  closeness  of  the  data  sets.  The  mean 
streamwise  velocities  compare  quite  w'ell,  even  at  X/C  ~  1.05,  Z/T  =  0.0,  which  one 
would  not  (,xpeci.  'fhe  data  sets  slightly  disagree  at  some  locations  near  die  wall,  which 
is  due  to  a  shallower  slope  in  the  log  region  for  the  JLF  hot-wire  data.  This  increases  the 
values  of  L’tc  near  the  wall  for  the  hot-wire  data  by  a  few  percent  at  the  most. 

'Hie  differences  aic  greater  for  the  spanwise  velocities.  Tltere  are  4  factors  which 
influence  these  dilfcrcnces.  These  are  1)  variation  in  initial  probe  aligiment,  affecting 
the  free, stream  flow  angle,  2)  eiTors  in  aligning  the  hot-wire  probe  from  point  to  point, 
3)  differences  in  vertical  position,  and  4)  apparent  hot-wire  probe  interference  near  the 
wall,  which  is  seen  as  odd  variations  in  Wf.(.  as  Y  decreases.  Fur  the  most  par,  the 
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effects  of  probe  interference  appear  to  be  rather  limited.  From  the  hodograph  plots 
(Figure  4.77),  some  differences  in  the  relationship  between  Upg  and  Wpg  can  be  observed, 
which  would  tend  to  rule  out  vertical  positioning  as  a  source  of  the  discrepancies. 

The  streamwise  and  spanwise  rms  normal  stress  values  are  always  less  for  the  JLF 
data  near  the  wall  if  the  hot-wire  and  LDV  data  do  not  agree.  These  differences  are 
considerable  for  the  Wtc/Uref  data,  especially  in  plane  5  (15%  variations).  The  general 
agreement  for  the  w'-xc/^rti  much  better,  usually  within  5%.  Section  2.5.1 

discusses  some  inherent  errors  in  the  hot-wire  data  reduction  equations.  Judging  by  the 
results  of  the  corrections  in  planes  13-15,  where  the  local  turbulence  levels  were  much 
higher  than  near  the  wing,  the  data  reduction  errors  can  be  ruled  out  as  a  cause  of  the 
observed  differences.  A  likely  cause  is  a  damping  effect  of  the  hot-wire  probe  on  the 
flow  as  it  nears  the  wall,  which  would  lower  the  measured  turbulent  fluctuations.  This 
hypothesis  is  supponed  by  the  apparent  linear  growth  in  the  variations  with  ln(Y/T),  a 
functional  form  which  is  present  in  the  solution  describing  a  laminar  flow  near  an 
oscillating  wall,  which  is  analogous  to  a  turbulent  flow  near  a  stationary  wall  (Stokes 
(iS5i),  and  Van  Driest  (1965),  both  referenced  in  Schetz,  iyS4).  Pursuing  this  discussion 
is  beyond  the  scope  of  this  research. 

The  similiarities  in  MvvTc/(Un.()^  are  remarkable  considering  the  uncertainties 
involved.  However,  a  hot-wire  probe  damping  effect  is  evident  in  Figures  4.70  and  4.75 
as  well.  These  plots  show  that  the  turbulent  magnitudes  obtained  using  a  hot-wire  were 
usually  less  than  those  obtained  using  a  LDV. 
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5.3.6  Comparison  to  Ailingeris  Cj-  Results 

The  Cf  results  obtained  using  oil- flow  interferometry  by  Ailinger  (1990)  were 
converted  to  frcestreani  coordinates  and  compared  to  lesults  obtained  using  equation 
(4.10)  (Figure  4.93).  Tlie  wail  skew  angles  reported  by  Ailinger  were  used  to  convert  his 
data  to  freestream  coordinates. 

Knowing  that  this  is  a  highly  3-(iimensional  flow,  witli  separation  regions  and 
large  spanwise  velocities,  one  would  not  expect  equation  (4.10)  to  estimate  Q  to  within 
10-15%  of  the  actual  value.  Also  knowing  that  Ailinger  reported  some  difficulties  in 
obtaining  low  uncertainty  skin-friedon  data,  especially  near  the  line  of  low  shear,  it  is 
surprising  tliat  the  C.  results  compare  usually  to  within  10%.  The  largest  difference  in 
skin  friction  estimates  occurs  in  the  highly  skewed  flow  inboard  of  the  LOLS  in  Plane 
5.  A  cautious  optimist  would  view  the  good  agreement  between  Cf  estimates  as  an 
indication  of  the  robustness  of  skin-ftiction  correlations  based  on  Rcq  and  the  shape 
factor,  H.  A  Q  estimate  accurate  to  within  10%  for  highly  3-D  flows,  based  on  a 
relatively  simple  formulation,  would  certainly  aid  rough  engineering  calculations. 
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6.0  Vortex  Panel  Calculations 


In  order  to  get  a  feel  for  the  qualitative  differences  inherent  in  horseshoe  vortices 
of  various  appendage  shapes,  a  vortex  panel  calculation  was  used  to  examine  the  inviscid 
flow  around  several  different  appendage  geometries.  This  chapter  discusses  the  goals, 
methods,  and  results  of  the  author’s  numerical  research. 


6.1  Numerical  Research  Program  and  Objectives 

As  mentioned  in  section  1.2.2,  the  distribudon  and  circulation  present  in  the 
horseshoe  vortex  is  affected  by  the  appendage  in  two  ways.  The  geometrical  bluntness 
increases  die  average  vortex  stretching  rate,  which  increases  the  peak  vorticity  and  affects 
the  vortex  structure.  The  appendage  bluntness  also  affects  the  initial  vortex  roli-up,  whic  t- 
increases  the  overall  strength  of  the  secondary  flow.  The  author’s  findings  supi a.';1 
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Mehta’s  (1984)  earlier  research  on  the  effects  of  nose  shape  on  the  horseshoe  vortex.  He 
found  that  a  non-dimensionalized  estimate  of  the  streamwise  circulation,  defined  as 

r  =  f®  dr  (6.1) 

^  Jo 

increases  as  the  appendage  nose  becomes  blunter.  In  evaluating  equation  (6.1),  Mehta 
used  the  same  Z  location  for  all  nose  shapes,  chosing  one  which  passed  near  the  vortex 
center.  His  results  were  used  to  obtain  a  correlation  between  apperdage  nose  geometry 
and  the  horseshoe  vortex  strength. 

One  o^*  the  main  quantities  of  interest  is  tlie  partial  derivative,  3(Vs  )/3S,  which  is 
the  inviscid  vortex  stretching  rate  in  the  direction  locally  tangent  to  the  body  surface.  Vg 
is  the  local  tangential  velocity,  and  .S  represents  the  disLonce.  or  arc  length,  along  the 
body’s  surface.  Integrating  this  derivative  around  the  nose  to  the  maximum  thickness: 


^Tmax'^o 


!. 


Dnax 


dS 


(6.2) 


Tmax 


one  can  calculate  the  average  vortex  stretching  rate.  A  major  part  of  the  numerical 
research  work  was  spent  examining  the  dependence  of  this  quantity  on  the  appendage 
nose  geometry.  Calculating  equation  (6.2)  was  performed  numerically  using  a  2-D 
inviscid  vortex  panel  method  from  Kuethe  and  Chow  (1986).  This  method  uses  panels 
with  linearly  varying  vortex  strengths  and  solves  for  no  flow  through  the  surface 
(tangential  flow  only)  at  each  panel’s  control  point.  The  control  points  bisected  their 
respective  panels.  The  computer  code  from  Kuethe  and  Chow  output  the  tangential 
velocity  (Vg )  and  coordinate  infonnation  for  each  panel’s  control  point. 
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The  last  numerical  research  objective  involves  the  temporal  flow  characteristics 
of  tlie  wing-body  junction  flow.  Rood  (1984)  concluded  i.hat  the  frequency  of  large  scale 
stmctures  in  the  junction-appendage  wake  flow  ?je  dependent  on  the  appendage  nose 
shape.  Also,  many  researchers,  notably  Devenport  and  Simpsc.i  (1990b),  have  shown 
that  a  bimodal  flow  structure  exists  in  the  stagnation  region  of  the  appendage-body 
junction.  The  bluntness  and  resulting  velocity  gi'adients  at  the  stagnation  point  may 
greatly  affect  this  bimodal  flow  structure.  The  ratios 


/  / 

I  as  I 

form  two  possible  Strouhal  numbers  characteriiing  the  bimodal  structure,  where  /  is  the 
bimodal  switching  frequency.  The  subscripts  o  and  Tmax  indicate  where  the  expressions 
are  evaluated.  Tmax  indicates  the  maximum  thickness,  o  indicates  the  appendage  2-D 
stagnation  point.  These  parameters  may  also  be  related  to  the  surface  pressure 
fluctuations  along  the  stagnation  line  of  the  wing-body  flow.  While  f  must  be  determined 
by  experimental  data  (surface  pressure  fluctuations  or  velocity  measurements),  the  terms 
in  the  denominators  of  equation  (6.3)  can  be  estimated  numerically,  aiding  in  future 
research  of  the  junction  temporal  flow  characteri.stics. 


(6.3) 
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6.2  Appendages  Shapes  Studied 


A  total  of  30  appendage  shapes  were  studied  at  0°  angle  of  attack,  while  15  shapes 
were  studied  at  12°  angle  of  attack.  Only  the  geometries  with  rounded  leading  edges 
were  studied  at  12°  angle  of  attacK,  since  a  sharp  leading  edge  would  cause  separation  in 
a  viscous  flow.  The  results  from  an  inviscid  calculation  in  such  a  case  would  be 
meaningless.  The  appendage  shapes  studied  were 

•  NACA  OOXX  family  (0006  through  0024) 

•  Mehta’s  body  shapes  (1984)  with  faired  on  tail  sections 

•  688  submarine  sail  cross  section 


•  CircuIsT  nose  with  strs2.ir!lincd  fsiring 

•  3:2  elliptic  nose  with  0020  tail 

•  Sand  1850  airfoil  section 

•  Shin’s  body  shapes  (1989) 

The  modelling  of  the  appendages  was  done  one  of  two  different  ways.  For 
Mehta’s  body  shapes  and  the  circular  nose  shape,  analytical  e.xpressions  were  used  (Z  = 
f(X))  to  achieve  the  final  panel  resolution  of  132  panels.  For  all  the  other  shapes, 
coordinate  data  was  input,  and  then  cubic  splines  were  used  to  obtain  the  final  panel 
resolution.  In  all  cases,  a  modified  circular  arc  point  distribution  was  used  to  establish 
the  final  X/C  values.  This  method  automatically  clusters  control  poi '  ‘.s  at  the  leading  and 
trailing  edges,  where  higher  densities  are  needed.  Figui'es  6.  l(a-e)  show  the  results  of  the 
panelling  schemes.  Errors  inherent  in  this  numerical  modelling  process  resulted  in 
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geometrical  deviations  from  the  actual  appendage  shapes.  Geometric  characteristics  such 
as  maximum  thickness  location  and  leading  edge  radius  (Plb)  are  usually  in  error  by  at 
least  few  percent  These  errors  are  not  significant  when  examining  the  relationship 
between  numerically  calculated  quantities  (as  in  section  6.3.2).  When  applying  the  same 
relationships  to  experimental  data,  however,  one  should  be  careful  to  use  the  correct 
geometric  parameters. 

Mehta  and  Shin  both  used  equations  in  tlie  fonn 


to  define  the  nose  shapes  of  their  "wedge  elliptic"  geometries.  If  n  <  2,  the  nose  comes 
to  a  point,  while  increasing  n  leads  to  a  blunter  nose.  Figure  6.2  shows  how  n  affects  the 
nose  shape  for  these  appendages.  In  Mehta’s  experiment,  the  actual  appendages  had 
"chopped"  tails  (no  trailing  edges).  A  tail  was  numerically  faired  onto  these  appendage 
shapes  for  the  computations,  utilizing  equation  (6.4).  The  tail  length  was  set  at  200  mm 
(a  =  200)  and  n  was  always  1.2.  For  the  numerical  work,  values  of  n  from  1,5  to  3.5, 
and  T/C  ratios  of  5  and  10  percent  were  used  to  generate  a  "family"  of  Mehta-like 
appendage  shapes.  To  calculate  the  nose  geometry,  the  actual  values  that  Mehta  used, 
a  =  150  mm  and  b  =  21  mm,  were  used.  Shin  (1989)  used  values  of  1.5  and  2.0  for  his 
two  wedge  elliptic  shapes,  witli  a  =  142.24  mm  and  b  =  35.56  mm.  The  appendages  used 
had  NACA  0020  tails  faired  on  to  the  nose  sections  at  the  maximum  thickness,  and  the 
coordinates  listed  in  Shin  (1989)  were  used  to  define  the  appendage  geometry. 
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6.3  Vortex  Panel  Calculation  Results 


To  accomplish  the  goals  set  for  this  work,  the  inviscid  velocity  distribution  around 
the  nose  of  varying  airfoil  sections  was  examined  by  calculating  8(Vs  )/3S,  which  controls 
the  horseshoe  vortex  stretching  rate.  This  derivative,  and  the  derivatives  used  to 
approximate  ,  were  estimated  numerically  using  central  differencing.  After  much  trial 
and  error,  a  correlation  between  the  average  vonex  stretching  rate  and  appendage 
geometric  parameters  was  established.  A  relationship  between  a  vortgx  transpon 
parameter  (\'TP)  and  Mehta’s  experimental  circulation  estimates  (T^ )  was  also  found. 

6.3.1  Preliminary  Numerical  Results 

To  verify  the  vortex  panel  code  results,  the  output  velocity  distribution  results 
from  the  NACA  OOXX  sections  were  compared  to  wind  tunnel  data  from  Abbott  and  Von 
Doenhoff  (1959).  As  Figures  6.3(a-f)  indicate,  the  inviscid  approximation  is  quite 
accurate  compared  to  wind  tunnel  data,  especially  for  the  thinner  sections.  The  largest 
error  encountered  was  about  1.5%  of  the  maximum  value  of  ,  This  error,  which 

increases  with  appendage  thickness,  is  likely  due  to  the  presence  of  the  boundary  layer 
in  the  viscous  flow.  The  effect  of  the  boundai^  layer  is  to  increase  the  apparent  thickness 
of  the  wing,  therefore  causing  the  ’ocal  velocity  to  increase.  The  spurious  results  at  the 
trailing  edge  seem  to  be  affected  by  the  relative  length  of  the  trailing  edge  panels.  The 
shorter  these  panels  became,  the  greater  the  jump  in  the  calculated  Vs/U^cf .  This  anomaly 
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was  judged  to  be  inconsequential  with  regard  to  the  final  results,  as  the  nose  region,  in 
the  range  0  <  X/C  <  ,  was  the  region  of  interest. 

Figures  6.4  and  6.5  present  the  velocity  distributions  for  the  other  appendage 
shapes.  Figures  6.6(a-j)  show  the  results  for  the  12°  angle  of  aitack  cases.  The  sharp 
peaks  in  the  velocity  distribution  for  the  thin  appendages  are  unlikely  to  occur  in  the 
presence  of  a  viscous  fluid.  The  astute  observer  will  notice  that  the  velocity  distributions 
around  the  688  sail  and  Sand  1850  sections  are  not  as  smooth  as  the  other  appendage 
results.  This  is  because  the  coordinates  for  these  appendages  were  obtained  by  graphical 
methods,  and  small  inaccuracies  were  inevitable. 


6,3^  Average  Vortex  Stretching  Rate  vs.  Body  Geometry 

The  results  of  attempting  to  correlate  the  average  vonex  stretching  rate  (non- 
dimensionalized  by  Plh/T1»)  and  appendage  nose  geometry  are  shown  in  Figui'es  6.7(a-d). 
This  correlation  involves  only  numerically  calculated  parameters.  The  quantity  BF  is 
defined  as 


BF  = 


1  P,.,  (  T _  /  ^ 


(6.5) 


This  "bluntness  factor"  was  found  through  trial  and  error,  and  represents  a  great  deal  of 
calculations.  More  rudimentary  parameters  were  tried,  without  much  success.  Examining 
its  terms,  it  takes  into  account  the  nose  leading  edge  sharpness  (Pl^  )?  the  "stubbiness"  of 
the  appendage  nose  ((S/X>j-^j),  and  an  indicator  of  the  average  appendage  surface  slope 
(T^„ /S-imax  )•  For  the  appendages  with  pointed  leading  edges,  p^;  was  not  set  at  zero. 
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but  was  estimated  numerically  using  the  control  points  generated  by  the  inviscid  code. 
The  results  of  the  estimate,  however,  are  still  quite  small.  The  correlating  parameters  of 


the  non-dimensional  stretching  to  appendage  geometry  was  detennined  for  a  -  0'^ 
independent  of  a  =  12”,  When  calculated  for  the  a  =  12”  cases,  the  average  vortex 
stretching  rate  was  calculated  from  tlie  stagnation  point  around  the  suction  side,  to  the 
maximum  thickness  location.  At  a  =  12°,  the  correlation  still  held  remarkably  well,  but 
not  quite  as  well  as  for  the  0°  cases.  It  should  be  noted  that  BF  is  not  a  function  of  the 
angle  of  attack. 

Three  types  of  curve  fits  were  applied  to  the  VsPle/SU,,  vs.  BF  data,  linear  (y  = 
mx  +  b),  and  two  types  of  log  curve  fits  Gog(y)  =  m  log(x)  +  b).  The  log  curve  fits  were 
minimized  on  linear  least  squares,  (erroi^  =  (y  -  10^  x""  j  unu  log  least  squares  (error” 


=  (log(y)  -  m  log(x)  -  b)^.  The  linear  least  square  fit  was  done  to  minimize  the  error  of 
the  curve  fit  at  high  values  of  BF.  Table  6. 1  shows  the  values  obuiined  from  these  three 
curve  fits.  The  log  curve  fit  based  on  linear  least  squares  gives  the  best  results  in  terms 
of  the  correlation  coefficient  (Pxy)>  variance  estimate  (a).  Figures  6.7(a-d)  also  show 


the  icSUliS  of  llic  Cui'vc  fils. 


Velocity  profile  data  from  Shin’s  (1989)  three  appendage  nose  geometries  were 
used  to  evaluate  the  {di\W)/dY)d/U^[  distribution  at  the  maximum  thickness  location  for 
each  geometry.  The  span  wise  location  of  the  profiles  was  ZJT  =  -.901,  approximately  the 
vortex  center  location  as  defined  by  the  peak  vorticity  contour.  Figure  6,8  shows  the 
results  of  the  vorticity  calculations.  The  baseline  geometi7  has  by  far  the  largest  vorticity 
levels.  Note  that  the  vertical  location  of  maximum  vorticity  for  each  case  roughly 
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coincide  with  one  another,  revealing  that  this  dimension  may  only  be  weakly  effected  by 
nose  shape.  Also  note  the  steeper  gradient  in  vorticity  for  the  baseline  case,  indicating 
a  greater  relative  concentration  of  vorticity.  Figure  6.9  displays  how  the  maximum 
vorticity  increases  with  BF.  This  is  an  expected  result  since  increasing  BF  increases  the 
average  vortex  stretching  rate.  However,  with  the  limited  number  of  data  points,  a 
definitive  correlation  cannot  be  found. 

6.3.3  Stream  wise  Circulation  vs.  Vorticity  Transport  Parameter 

The  correlation  of  average  vortex  stretching  rate  vs.  BF  involves  only  numerical 
data,  but  the  second  correlation  found  relates  a  numerically  estimated  quantity  to 
experimental  results.  Mehta  (1984)  derives  a  simplified  vorticity  transport  equation  which 
models  the  transport  of  transverse  vorticity  into  the  streamwise  direi,tion.  By  ignoring 
viscous  terms,  the  convection  of  vorticity  by  V  and  W,  and  the  contributions  due  to  0>y, 
Mehta  presents  the  following  equation: 

(6.6) 

dx[  u )  f/2  ax 

This  equation  shows  that  geometries  generating  large  values  of  a(W)/9X  and  small  values 
of  U  will  create  the  strongest  horseshoe  voriicies.  An  appendage  with  a  high  bluntess 
factor  meets  both  these  criteria. 

From  equation  (6.6),  a  vortex  transport  parameter  (VTP)  can  be  derived,  relating 
the  transported  streamwise  vorticity  to  inviscid  quantities.  Integrating  botli  sides  of 
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equation  (6.6)  along  the  X  direction,  and  setting  W  =  V^sin(9),  where  9  is  the  local  panel 
angle,  the  following  equation  is  obtaired: 

~  dCF^sinCe))  (6.7) 

By  assuming  that  co^  «=  Ug/S,  and  realizing  that,  U  =  Ue  =  VsCos(9),  the  following 
relationship  can  be  derived: 


diVs  sin(0)) 


(6.8) 


Evaluating  the  right  hand  side  of  equation  (6.8)  from  the  stagnation  point  to  the  maximum 
thickness  location  forms  the  vorticity  transport  parameter  (VTP),  which  is  related  to 

Figure  o.iO  shows  the  correlation  between  Mehta’s  non-dimensional  circulation 
estimate  and  tlie  VTP.  Note  that  the  results  for  the  5%  and  10%  thick  appendages  are 
shown.  To  correlate  the  two  cases  with  different  thickness  ratios,  was  divided  by 
3(Vs/U,2f)/3(S/C)  evaluated  at  the  stagnation  point.  This  quantity  is  defined  as  F'x. 
Doing  this  docs  not  correlate  tlie  two  set  of  data  completely,  as  the  tv/o  lines  have  slightly 
different  slopes,  r'x  was  found  to  have  some  dependence  on  the  bluntness  factor  as  well. 
The  results  in  Figure  6.10  are  only  preliminary,  since  not  much  data  is  available  for 
different  appendage  shapes.  The  reader  should  note  that  F'^  depends  strongly  on  the 
approach  boundary  layer  properties,  and  Mehta’s  results  correspond  to  one  particular 
boundary  layer.  Mehta’s  experimental  values  for  F^  were  .076,  .102,  and  .124  for  n  = 
1.5,  2.0,  and  3.0,  respectively.  The  experimental  parameters  were  U„f  =  25  m/s,  6/T  = 
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.595,  and  Rcq  =  58,000  in  the  undisturbed  boundary  layer  at  the  wing  leading  edge 
location. 

Fx  and were  calculated  for  Shin’s  geometries,  using  the  (3(W)/8Y)  component 
of  Qx  results  tabulated  in  Table  6.2.  Results  for  F'x  are  plotted  in  Figure 

6.10  as  well.  One  can  see  a  different  relationship  between  F'x  and  VTP  for  Shin’s  data. 
This  variation  is  due  to  the  different  approach  boundary  layer,  and  may  also  be  due  to  the 
different  geometries  involved. 

Appendix  C  lists  all  results  in  tabular  form  for  both  the  0”  and  12°  angle  of  attack 
cases,  including  the  ratios  relevant  to  equation  (6.3). 
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7.0  Conclusions 

This  chapter  attempts  to  collect  and  list  the  most  important  and  significant 
conclusions  of  this  study.  When  discussing  the  results,  the  phrase  in  general  was  used 
quite  often  to  describe  certain  trends  or  tendencies,  'fhis  wording  was  used  because  a 
contradiction  to  any  statement  concerning  the  data  comparisons  could  be  found  if  the  data 
is  examined  long  enough.  The  conclusions  were  based  on  observations  of  all  the 
available  data,  and  reflect  the  prevalent  trends  discovered. 

7.1  Basic  Flow  Physics 

Conclusions  concerning  the  basic  flow  physics  resulted  from  examining  the 
author’s  hot-wire  data  and  Devenport  and  Simpson’s  (1990a)  LDV  data.  These 
conclusions  are: 
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•  The  primary  horseshoe  vortex  flow  structure  is  elliptically  shaped  and  near  the 
wall,  similiar  to  what  Dickinson  ( 1986a, b)  observed.  This  characteristic  shape  may  be 
due  to  vortex  meandering,  caused  by  the  unsteadiness  in  the  nose  bimodal  flow  region. 

•  The  vorticity  in  the  junction  vortex  is  dominated  by  the  change  in  spanwise 
velocity,  and  the  term  (8(W)/3Y)  defines  the  characteristic  distribution  of  vorticity  in  the 
flow. 

•  The  wall  no-slip  condition  creates  a  very  thin  region  of  high  negative  vorticity 
near  the  wall  underneath  the  primary  vortex,  which  thickens  outboard  of  the  LOLS. 

•  The  line  of  low  shear  (LOLS)  position  defines  locations  of  peak  turbulence 
intensity  and  boundary  layer  thicknesses.  The  time-averaged  position  of  the  horseshoe 
vortex  core  is  inboard  of  the  LOLS,  and  coincides  with  the  location  of  large  spanwise 
gradients  in  Urc/U^f. 

•  A  quantity  termed  the  "distortion  function"  (fp)  reveals  the  effects  of  the  secondary 
flow  distortions,  which  enables  more  conclusive  flow  comparison  results  in  chapter  5. 

•  From  examining  and  drawing  conclusions  from  the  results  of  Shizawa  and  Eaton’s 
(1990)  study,  the  importance  of  the  geometry-induced  skewing  was  seen.  I  he  conflicting 
directions  of  tlie  skewing-induced  and  vonex- induced  flows  near  the  wall  controls  th(. 
flow  separation.  Therefore,  the  skewing  affects  the  growth  of  the  streamwise  mean 
velocity  distortions,  and  also  the  secondary  flow  structure  and  associated  vorticity 
distribution. 

•  Flow  distortion  increases  with  X/C,  and  grows  rapidly  when  enhanced  by  a 
favorable  boundary  layer  skewing  direction.  Adverse  pressure  gradients  also  enhance  the 
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giowth  of  flow  distortions.  The  flow  over  the  aft  portion  of  the  appendage  subjects  the 
horseshoe  vortex  to  skewing  and  a  mild  adverse  pressure  gradient,  which  work  together 
to  increase  the  boundary  layer  distortions. 

•  The  characteristic  peaks  in  streamwise  normal  stresses  can  be  explained  by  the 
increased  mixing  of  the  boundary  layer  fluid,  and  a  destabilizing  of  the  boundary  layer 
caused  by  the  lifting  action  of  the  secondary  flow. 

•  The  far  wake  flow  distortion  patterns  reveal  a  more  diffused  vortical  fiow 
structure.  However,  the  growing  boundary  layer,  particularly  in  the  adverse  pressure 
gradient,  is  the  primary  cause  of  flow  distortions  growth  X/C  increases. 

•  "Hre  spanwise  location  of  maximum  boundary  layer  thickness  (and  minimum  skin 
friction)  remains  fairly  constant  at  jZ/Tj  -1.8  for  X/C  values  greater  than  6. 

7,2  Effects  of  the  Approach  Boundary  Layer 


The  conclusions  of  trie  main  focus  of  this  study  are  rather  qualitative  in  na'  '.'re,  but 
this  does  not  diminish  their  significance.  These  conclusions  were  obtained  by  comparing 
the  data  sets  listed  in  Tables  1.1  and  3.3,  and  are  as  follows: 

•  The  secondary  flow  structure  and  distortion  function  (fp)  scale  on  T  in  the  Y  ind 
Z  ducctions.  This  fact  reveals  that  die  horseshoe  vortex  flow  structure  is  primarily  diiven 
by  the  aj  pendage  geometry  and  associated  pressure  gradients. 


( 
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•  The  effects  of  Reg  can  be  easily  seen  in  the  nose  region  using  Kim’s  (1991)  LDV 
measurement  results.  The  larger  Reynolds  number  increases  the  local  mean  flow 
distortions  and  gradients  near  the  wall,  and  also  concentrates  and  amplifies  the  near  wall 
turburlence  levels. 

•  From  the  observing  the  effects  of  Reg  in  the  nose  region,  the  author  postulated 
that  these  effects  may  be  more  effectively  described  by  a  parameter  termed  the  momentum 
deficit  factor  (MDF).  liiis  factor,  equal  to  (ReT)^(0/T),  is  indicative  of  the  relative 
decrease  in  momentum  flux  near  the  wall.  Flows  with  large  values  of  MDF  would 
e-xpc-rience  an  enhanced  near  wall  jetting  of  the  secondary  flow  due  to  the  "relief  valve" 
effect  of  near  wall,  low  momentum  fluid. 

•  Comparison  of  the  data  sets  reveals  that  the  MDF  correctly  predicts  the  variation 
in  mean  flow  distortion  magnitudes  and  associated  horseshoe  vortex  characteristics.  Rcg 
inccmectly  predicts  the  relative  distortion  magnitudes. 

•  The  two  vortex  configurations  studied  by  Shizawa  and  Eaton  (1990)  were  ftmnd 
to  have  chaxacteiistics  which  corresponded  to  flows  with  large  MDF  (case  1)  and  flows 
with  small  MDF  values  (case  2).  One  can  conclude  that  the  effects  of  the  relative 
momentum  near  the  wall  are  to  modify  the  effective  flow  skewing  around  the  appendage. 
Low  r/IOF  increased  the  effective  flow  skewing,  and  high  MDF  decreased  the  effective 
flow  skewing. 

•  As  MDF  increases,  the  following  effects  are  observed: 

-•  the  mean  streamwise  velocity  distortions  ai'e  not  as  large 
-  the  secondary  flow  patterns  are  more  elliptic 
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-  the  positive  vorticity  is  increasingly  concentrated  in  a  near  wall  region 

-  the  regions  of  negative  vorticity  outboard  of  the  LOLS  are  less  significant 

-  the  movement  of  the  vortex  core  away  from  the  wall  is  decreased 

-  the  downstream  spanwise  location  of  the  vortex  core  is  slightly  further  from  the 
plane  of  symmetry 

•  The  vorticity  and  helicity  levels,  non-dimensionalized  by  6/U„f,  are  greater  for  the 
high  MDF  flows.  This  may  be  due  in  part  to  increasing  Ure/T,  which  affects  the  average 
vortex  stretching  rate,  and  TAJ^f  may  be  a  more  appropriate  scaling  factor  for  vorticity. 
Definite  conclusions  regarding  the  vorticity  scaling  could  not  be  reached  due  to  the 
variation  of  MDF,  SAJren  ^d  TAJref  between  the  data  sets. 

•  The  effects  of  decreasing  5*/7.’  arc  seen  primarily  by  increased  boundaiy'  layer 
thinning  between  the  appendage  surface  and  the  line  of  low  shear.  The  increase  in 
relative  size  of  the  vortex  to  the  boundary  layer  height  increases  the  availability  of  high 
momentum  fluid,  leading  to  larger  skin  friction  values.  Also,  die  contours  of  m  xc/^ref 
show  more  distinct  local  maxima  for  lower  values  of  6*/T,  due  to  the  initially  thin  layer 
of  turbulent  fluid, 

•  In  general,  the  turbulence  quantities  show  good  agreement  between  the  data  sets. 
The  higher  Rcq  flows  show  only  slight  increases  in  strcamwisc  turbulence  intensity  near 
the  wall,  which  is  surprising  in  light  of  Kim’s  (1991)  observations  at  the  appendage  nose. 
'*  The  quanitities  u  tc>  v'xc-  ^^^d  mv^  all  appear  to  scale  on  T  in  the  venical 
direction,  while  w  jc  and  tJvxc  seem  to  scale  more  appropriately  on  5. 

•  The  maximum  spanwise  variations  of  S*  and  0  scale  on  T. 
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•  The  LDV  results  of  Devenport  and  Simpson  (1990a)  generally  compare  well  to 
the  author’s  hotwire  results  (to  within  2-3%).  The  near  wall  turbulence  measurements 
disagree,  which  may  be  due  to  damping  of  the  velocity  fluctuations  by  the  hot-wire  probe. 

•  White  and  Ludwieg-Tillmann’s  correlations  predis.  !  (Cf)ps  to  within  approximately 
5-10%  of  Ailinger’s  measured  values.  This  is  remarkabU  considering  the  highly  3- 
dimensional  nature  of  this  flowfield. 

To  gain  further  understanding  of  the  effects  of  the  appioach  boundary  layer  fiow, 
a  more  rigidly  controlled  experiment,  with  less  "interfering  inputs",  could  be  designed. 
This  experiment  would  use  the  same  appendage  geometry,  facii.des  and  experime -tal 
apparatus.  The  only  variations  would  be  in  the  approach  flow  conditions. 


7.3  Numerically  Predicted  Effects  of  Body  Geometry 

The  vortex  panel  calculations  produced  namcrical  results  which  revealed  a  definite 
relationship  between 

•  average  vortex  stretching  rate  and  a  parameter  describing  the  geometric  bluntness 
of  the  body  shape  (the  blunmess  factor,  BF) 

•  streamwise  horseshoe  vortex  circulation  and  a  vortex  transport  parameter  (VTP) 
Calculations  were  also  done  to  aid  in  the  further  study  of  the  junction  temporal  flow 
characteristics. 

All  of  the  numerical  results  are  preliminary,  and  more  work  needs  to  be  done  to 
fully  understand  the  effects  of  nose  shape  on  the  junction,  flow.  In  light  of  the  author’s 
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experimental  conclusions,  the  examination  of  the  2-D  potential  flow  skewing  away  from 
the  body  may  yield  further  insights  into  the  effects  of  geometry.  Any  joint 
experimental/numerical  study  of  various  appendage  shapes  would  certainly  produce 
valuable  information.  An  even  more  in  depth  program  could  further  study  and  account 
for  the  effects  of  changing  the  approach  boundary  layer. 
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Table  1.1. 


Comparison  of  measurement  locations  and  quantities  measured,  3:2  elliptic 
nose,  NACA  tail  junction  flow  studies 


Fleming 

(JLF) 


.18  (5) 

.42  (6) 

Measurement  .64  (7) 
Planes  (X/C) 

.93  (9) 
1.05  (10) 

I. 50(11) 

5.89  (13)* 
9.14(14)* 

II. 56  (15) 


Dickinson 

McMahon 

(SCD) 

et  al. 

.18  (5) 

(HMM) 

.64  (7) 

.75  (8) 

.75  (8) 

.93  (9) 

1.05  (10) 

1.00  (10) 

1.50  (11) 

3.00  (12) 

Shin 

Devenport 

(JS) 

et  al. 
(WJD) 

.18  (5) 

18  (5) 

.43  (6) 

.75  (7) 

1.00  (10) 

1.05  (10) 

1.50  (11) 

3.00  (12)* 

Notes: 

JLF  measurements:  U,  u',  W,  w\  -uw  taken  with  single,  yawable  hot-wire,  except 
locations  designated  ()*  indicate  measurements  of  U  and  u'  only. 

SCD  measurements:  U,  u',  W,  w',  V,  v',  -mv,  -mv  taken  with  hot-film  probe. 
Dickinson  also  found  mean  pressures  on  wall  and  made  oil-flow 
visualizations. 

HMM  measurements:  U,  u',  W,  w',  V,  v',  -tZP,  -uw,  -m  taken  with,  rotating,  slant- 
wire  probes. 

JS  measurements:  U,  V,  W  taken  with  5-hoIe  probe. 

WJD  measurements:  U,  u',  W,  w',  V,  v',  -ilV,  -Uw,  -vvv,  taken  with  LDV 
(Devenport  and  Simpson,  1990a).  Location  designated  ()*  indicates  hot¬ 
wire  measurements  of  U  and  u  only  (Devenpon  et  a/., (1990)).  Also, 
measurements  were  made  for  botii  thick  and  thin  approach  boundary  layers 
at  X/C  =  3.00. 
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Table  2.1. 


Offsets  for  3:2  elliptic  nose,  0020  tail  appendage 


x/c 

Y/C 

0.0 

0.0 

0.015 

0.047 

0.029 

0.065 

0.059 

0.088 

0.088 

0.102 

0.118 

0.111 

0.176 

0.118 

0.294 

0.114 

0.412 

0.104 

0.529 

0.089 

0.647 

0.072 

0.765 

0.051 

0.882 

0.028 

1.0 

0.0 
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Table  2.2.  X,  Y,  and  Z  location  uncertainties  (95%  confidence  levels) 


Uncertainty  in  X: 
Uncertainty  in  X/T; 

Uncertainty  in  Y: 
Uncertainty  in  Y/T; 
Uncertainty  in  AY: 

Uncertainty  in  Z: 
Uncertainty  in  Z/T: 


±(1  +  (Y)sin(2°))  (mm) 
±(.014  +  (Y/a’)sin(2°)) 

±  .08  mm 
±  .001 
±  .03  mm 

(Same  as  X) 

(Same  as  X/T) 
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Table  2,3.  Nominal  experimental  uncertainties  (95%  confidence  levels) 


Temperature  drift  and  uncertainty; 
Pressure  drift  and  uncertainty: 
Manometer  uncertainty: 

Amplifier  offset  drift: 

Amplifier  gain  drift: 

Axial  coefficient  uncertainty; 
Probe  angular  positioning  error: 
Mean  voltage  enor: 


±l^C 
±  5  mb 
±  .012  in 
±  .001  V  (±  .2%) 
±  .1  (±  .5%) 

±  .05 
±  .5° 

±  .5% 


Note:  Velocity  gradient  effects  were  not  considered. 
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Table  2,4.  Reduced  data  uncertainty  (95%  confidence  levels) 


Uncertainty  in: 


Uncertainty: 


u/u„f 

±  1.5% 

w/u,,f 

±  U/U„f  sin(.5°) 

±  1.5% 

±  5% 

-i3vvAJ^„f 

±5% 

«r.s 

±,5° 

a,. 

±  2.5°  (worst  casi 

Note:  Data  in  tunnel  coordinates 
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Table  2.6.  Quoted  experimental  uncertainties  of  data  sets  used  for  comparisons  of 
wing-body  junction  flows 


Quantity 

SCD’ 

WJD^ 

HMM^ 

JS" 

U/U,,f,  V/U„f,  W/U,,f 

±2.5% 

±.035 

±2% 

±.007 

N/A 

±.003 

±8% 

- 

-«v/(U,,f)^  -im-/(U„f)^ 

N/A 

±.0002 

±6% 

- 

fHow  angle 

N/A 

N/A 

±1.0° 

- 

Notes: 


These  figures  are  all  nominal  uncertainties 
’From  Dickinson  (1986a) 

^LDV  uncertainties  estimated  from  Devenport  and  Simpson  (1990a)  (note  also  that 
these  uncertainties  are  not  given  as  percentages) 

"Fstinijited  from  McMahon  et  al.  (1987)  (these  figures  are  the  upper  uncertainty 
limits) 

'‘Nominal  upper  limit  on  mean  velocity  uncertainties,  from  Shin  (1989) 
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Table  3.1.  Adverse  pressure  gradient  boundary  layer  profile  properties 


X/C  location 

-  -■ — 

10.135 

11.08 

12.06 

.853 

.817 

.782 

S.995  /T 

1.427 

1.718 

1.910 

d'rr 

.273 

.325 

.400 

e/T 

.152 

.166 

.184 

H  (=  870) 

1.80 

1.96 

2.18 

Cf 

.00142 

.00124 

.00101 

Notes: 

Q  estimated  using  Spalding’s  wall  law  (equation  3.2) 

Although  T  aird  C  do  not  have  any  physical  meaning  for  these  profiles  (wing  not 
mounted  in  tunnel),  the  results  are  non-dimensionalized  by  T  and  C  to  be 
consistent 
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Table  3.2.  Line  of  low  shear  locations 


Plane  #,  X/C 

(ZVT), 

5,  .18 

-.88* 

6,  .42 

-.99 

7,  .64 

-.94 

8,  .75 

-.90 

9,  .93 

-.87 

10,  1.05 

-.89 

Notes: 

*LOLS  location  in  plane  5  is  not  well  defined,  location  given  is  an  estimate. 
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Table  3.3. 


Approach  boundary  layer  parameters,  3:2  elliptic  nose,  NACA  0020  tail 
wing- body  flow  studies 


Data  Set  ID 

JLF,WJD‘ 

SCD^ 

HMM^ 

JS* 

WJD- 

T  (cm) 

7.17 

6.1 

7.1 

7.1 

7.17 

C  (cm) 

30.5 

25.9 

30.9 

30.2 

30.5 

T/C 

,235 

.234 

.230 

.235 

.235 

U„f  (m/s) 

26.7.5 

30.5 

15.24 

20.9 

32.0 

5/T 

.513 

1.197 

.947 

1.15 

.263 

6*/T 

.0779 

.1345 

.1467 

.148 

.043 

e/T 

.0548 

.1014 

.1003 

.1227 

.031 

H  (  =  570) 

1.42 

1.33 

1.46 

1.206 

1.387 

ri 

^  t  J' r\r\ 

iijDW  ! 

/'  r\f\r\ 

1 1 ,  / 

4,500 

Reg, 

8,900 

15,400 

10,000 

14,100 

6,200 

Re.j^ 

115,000 

1 14,400 

67,800 

95,400 

145,200 

Cf 

.00247 

.00247 

.00232 

.00256 

.0028 

U„f/T  (1/sec) 

373.1 

500.0 

214.6 

294.4 

446.3 

U„f/5  (1/sec) 

T113 

417.7 

226.7 

256.1 

1,697 

MDP  ('vin®'! 

1  - -  wf  y 

7  24 

1'^  'X 

4.6! 

■Bl 

A  ca 

- - - J| 

Notes: 

'JLF  hot-wire  and  WJD  LDV  and  hot-wire  data,  thick  boundary  layer  case,  approach 
boundary'  layer  profile  measured  at  X/T  =  -2.146,  Z/T  =  0.0,  with  wing  in  place. 

^SCD  hot-film  data,  approach  boundary  layer  profile  data  corrected  from  X/T  =  -3.167, 
Z/T  =  0.0  to  V/DJ  location,  wing  in  place. 

^HMM  hot-wire  data,  approach  boundary  layer  profile  data  corrected  from  X/T  -  -2.087, 
2jfY  =  0.0  to  WDJ  location,  wing  in  place. 

‘‘JS  5-hole  probe  data,  approach  boundary  layer  data  at  X/T  =  0.0,  no  wing  installed. 

*WJD  hot-wire  data,  thin  boundary  layer  case,  plane  12  (X/C  =  3.00)  only,  approach 
boundary  layer  profile  measured  at  X/T  =  -2.146,  Z/T  ■-  0.0. 
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p# 

Data  Set 
ID 

Z/T  limits 

0(W)/aY)MAx 

D 

A(5*A') 

A(eA') 

5 

JLF 

.775,1.525 

2.178 

2.137 

.0220 

.013 

WJD 

.55,1.525 

2.243 

2.049 

.0203 

.029 

SCD 

.55,1.525 

4.917 

.022 

6 

JLF 

.692,2.00 

1.112 

1.112 

.0320 

.021 

■ 

JLF 

.555,2.00 

.6583 

.483 

.0499 

.032 

SCD 

.555,2.00 

.8048 

.7003 

.0649 

.045 

8 

WJD 

.308,1.48 

.4637 

.391 

.0772 

.049 

SCD 

.308,1.48 

.5157 

.388 

.0617 

.044 

HMM 

.308,1.48 

.350 

.30 

.0939 

.065 

9 

JLF 

.327,2.00 

.3500 

.273 

.0818 

.048 

SCD 

.327,2.00 

.9767 

.670 

.0835 

.057 

JLF 

.200,2.00 

.300 

.200 

.0921 

.055 

WJD 

.200,1.60 

.300 

.200 

.0944 

.055 

SCD 

.200,2.00 

- 

.0894 

.063 

HMM 

.200,2.00 

.300 

.150 

,0772 

.063 

II 

11 

JLF 

.200,2.00 

.250 

.200 

.0933 

.061 

SCD 

.200,2.00 

.715 

.696 

.0991 

.064 

r 

WJD 

(thin) 

.354,4.603 

- 

.0969 

.0626 

WJD 

(thick) 

.354,4.603 

- 

.1086 

.0689 

HMM 


.354,2.833 


.1559 


.0919 


Table  6.1.  Curve  fitting  results  for  (Vj /S>pj^^PLEAJref  vs.  BF  data  (a  =  0°  and  12“) 


y  =  l(y’x'" 

log(y)=m  log(x)+b 

y  =  mx  +  b 

m 

.88157, .80749 

.94010,  .86617 

1.11.50,  1.0061 

b 

.048231, .014203 

.10727,  .057897 

.016204,  .03953 

PXY 

.99969,  .99806 

.99871,  .997165 

.99747,  .99356 

o 

.005777, .01755 

.031187,  .038464 

.016294,  .03201 

Notes: 


Columns  contain  both  0°  and  12°  cases,  separated  by  a  comma 
28  points  were  used  for  the  a  =  0°  data  (J.  Shin  appendages  not  used) 
14  points  were  used  for  the  ot  =  12°  data  (J.  Shin  appendage  not  used) 


i: 


Tables 


Table  6.2.  Effects  of  nose  shape  at  maximum  thickness  locations  for  Shin’s 
geometries 


Appendage  Geometry 

Tx 

w 
^  X 

0(W)/aY)MAX 

3:2  elliptic  nose,  NACA 
0020  taU 

.140 

.0066 

2.819 

4:1  elliptic  nose,  NACA 
0020  tail 

.0722 

.00134 

.9568 

4:1  wedge  elliptic  nose, 
NACA  0020  tail 

.0605 

.000068 

.7724 

Note:  is  non-dimcnsionalized  by  SAJ.ef- 
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Figures 
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of  Wing-Body  Junction,  Showing  Measurement  Planes  and  Tunnel 
nate  System 


Flow  (from  Devenport  and  Simpson,  1990a) 
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Figure  2.5.  Sketch  of  Diverging  Flow  Section  Used  to  Produce  Adverse  Pressure 
Gradient  Region  (not  to  scale,  all  dimensions  in  cm) 


Scoop 

Leading 

Edge 

Figure  2.6.  Portion  of  Oil-Flow  Visualization  From  Leading  Edge  of  Separation, 
Showing  Lack  of  Separation 
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Figure  2,7.  Variation  of  UAiref  in  Regions  I  and  2  of  Wind  Tunnel  (from  Devenport 
and  Simpson,  1990a) 
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Yjc>  ^FS»  ^TC»  ^FS 


^TC 


Figure  2.8.  Sketch  of  Tunnel  and  Frecstrcain  Coordinate  Systems 
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Yw. 


Ys.  Vs 


b) 


Arbitrary 
Coordinate 
System  (  ) 


Figure  2.9.  Sketches  of  a)  Wire  Coordinate  System  and  Velocity  Components  and  b) 
Local  Streamwise  Coordinate  Systerti  and  Velocity  Components 
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Figure  2.10.  Schematic  of  F,xperimental  Instrumentation 
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Probe  Support 


Figure  2.11.  Sketch  of  Hot-Wire  Probe  Traverse  System 
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d(Cp)/d(X/C)  vs.  X/C 
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V 
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Adverse  pressure  gradient  results 


□  —  Hot-wire  traverse  dato 
solid  line  —  1— D  inviscia  calculotion 
dashed  line  -  l-D  inviscid  calculation  allowing  for 
growth  of  boundary  layer  along  floor 


0 


# 


0 


0 


Figure  3.1.  Characteristics  of  the  Potential  Core  Row  in  the  Adverse  Pressure  0 

Gradient  Section  (Region  3) 
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Adverse  pressure  gradient  boundary  layer  profiles 


□  ~  JLr  data  (X/C  =  10.135,  Z/T  =  0.0) 

•  -  JLF  data  (X/C  =  11.08,  Z/T  =  0.0) 

o  -  JLF  data  (X/C  =  12.06,  Z/T  =  0.0) 

Not*:  data  taktn  in  empty  tunnel 


Figure  3.2.  Boundary  Layer  Piofiles  of  U/U,^  and  in  Adverse  Pressure  Gradient 
Section,  No  Appendage  in  Tunnel 
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Distributions  of  Time-Mean  Static  Pressure  Along  Z/T  =  ±3.17,  Compared 
to  Ponjntial  Flow  Calculation  (from  Devenport  and  Simpson,  1990a) 


Profiles  of  (w'/U„f)^  Measured  on  Either  Side  of  Wing,  1  mm  Downstream 
of  Wing  Trailing  Edge  (from  Devenport  and  Simpson,  1990a) 
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Figure  3.8.  ProHles  of  U/U,,f,  Measured  at  Various  Y/T  Locations  1  mni  Downsteam 
of  Wing  Trailing  Edge  (from  Devenport  and  Simpson,  1990a) 
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Figure  3.9.  Comparison  of  Profiles  of  Measured  at  Various  Y/T  Locations 

1  mm  Downstream  of  Wing  Trailing  Edge  (from  Devenport  and  Simpson, 
1990a) 


Figure  3.10a.  Results  of  Oil-Flow  Visualization,  Showing  Measurement  Locations  in 
Planes  5-10 


Figure  3. 1 L  Simf  lified  Sketch  Showing  Skewing  of  Outer  Streamlines  and  the  Primar}' 
Separation  (from  White,  1974) 
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Approach  boundary  layer  profiles 
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Y/T 


□  -  WJD  data  (X/T  =  -2.146,  Z/T  =  0.0) 

•  -  5CD  daiG  (X/T  =  -3.167,  Z/T  =  0.0) 

o  -  HMM  data  (X/T  =  -2.087,  Z/T  =  0.0) 


Figure  3.13.  i  .  .'lies  of  UAJ^/  for  the  Approach  Boundary  Layers  of  Devenport  and 
Simpson  (1990a),  Dickinson  (1986a, b),  and  McMahon  et  al.  (1987), 
(author’s  approach  flow  is  identical  to  WJD  data) 
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Approach  boundary  layer  profile  results 


Y*" 


a  -  WJD  da1u  (X/T  "  -2.U6.  Z/T  =  O.O) 
solid  line  -  Spalding  equation  for  inner  region 


Figure  3.1‘!a.  Spalding  Wall-Uw  Cm-ve  Fii  Applied  to  the  WJD  Approach  Bow 
Data 
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Approach  boundary  layer  profile  results 


□  -  SCD  dofc  (X/T  =  -3.167.  Z/T  -  0.0) 
solid  fine  —  Spalding  aquation  for  inner  region 


Figure  3.14b.  Spalding  Wall-Uiw  Curve  Fit  Applied  to  the  SCD  Approach  Flow 
Data 
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Approach  boundary  layer  profile  resuMs 
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Figure  3.Fk'.  Sp.ilding  Wall-Law  Curve  Fit  Ap 
Data 


plied  to  the  H.MM  Approach  Flow 
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Figure  4,1a,  Profiles  of  f3(W),/aY)6A4f  and  -0(V)/dZ)6AJ,,f,  Plane  5 
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Figure  4.2. 


Threc-Dimen.sional  Boundary  Layer  Pi-ofiJes,  Showing 
Skewing,  and  b)  Bi-directional  Skevnng  (from  White, 


a)  Uni-directional 
1974) 
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Freestream  date  from  plane  6,  X/C  =  .42 
□  -  JLF  dota 
O  -  JS  aata 

solid  line  -  2— D  inviscid  calculation  results 


Figure  4.4b.  Compaiison  of  Freestream  Mean  Velocity  and  Flow  Angle  Results  to  2-D 
Potertial  Flow  Calculation  Results  (Plane  6) 
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Freestream  data  from  plane  7,  X/C  =  .64 
D  “  JLF  data 
•  -  5CD  data 

solid  line  -  2-*D  inviscid  calculation  results 


Figure  4.4c.  Comparison  of  Frecstream  Mean  Velocity  and  Flow  Angle  Results  to  2-D 
Potential  Flow  Calculation  Results  (Plane  7) 
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Figure  4.4d. 


Comparison  of  Freestteam  Meaji  Velocity  and  Flow  Angle  Results  to  2-D 
Potential  Flow  Calculation  Results  (.Plane  b’) 
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0.4  0.8  1.2  1.6  2.0 

-Z/T 

Treestream  data  from  plane  9,  X/C  .93 
C  -  JLF  data 
*  -  SCO  data 

solid  lin«  ”  2  'D  inviscid  calculation  results 


Figure  4.4c.  Comparison  of  Frcestrcam  Mean  Velocity  and  Flow  Angle  Results  to  2-D 
Potential  Flow  Calculation  Results  (Plane  9) 
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Freest  re  am  data  from  piano  10,  X/C  =  1.05 
O  -  JLF  data 
A  -  WJD  LDV  data 
•  —  SCO  data 
O  '  JS  aata 
o  -  HMM  data 

solid  line  —  2  — D  inviscid  calculotion  results 


Figuie  4.4f,  Comparison  of  Freestream  Mean  Velocity  and  Flow  Angle  Results  to  2-D 
Potential  Flow  Calculation  Results  (Flare  10) 


E 


Figure  4.4g.  Comparison  of  Freestream  Mean  Velocity  and  Flow  Angle  Results  to  2-D 
Potential  Flow  Calculation  Results  (Plane  11) 
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Figure  4.5a. 


-Z/T 


Contours  of  U-ycA-^ru-  Pianes  5-11  (All  Jaia  uoin  present  study  except  XJC 
=  ,18  and  .75  from  Devenpon  and  Simpson  (1990a)) 
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Figure  4.6a.  Contours  of  Distonion  Function,  Planes  5-12  (All  data  from  present 
study  except  ^170  =  .18  .nnd  .75  from  Devenport  and  Simpson  (1990a)  and 
X/C  =  3.00  from  Devenport  ei  al.  (1990)) 
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X/C  -  .93 
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X/C  =  1.05 
Scale:  V  = 


X/C  =  1.50 
Scale:  1'  = 


Figure  4.7.  Secondar>'  Velocity  Vectors,  Planes  5-11  (All  data  from  present  study 
except  X/C  =  .18,  .75,  and  1.05  from  Devenport  and  Simpson  (I99()a). 
JLF  data  presents  only  W  component  of  secondaiy-  velocities) 
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Figure  4.11.  Profiles  of  Distortion  Function,  Plane  5 
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figure  4.13.  Pi'ofiles  of  Distortion  Function,  Plane  6 
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Figure  4.16.  Profiles  of  Utc/U^  Plane  8 
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g’jre  4.16  (cont.).  Profiles  of  LVc/Uc.  Plane  8 
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Figure  4.17.  Profiles  of  Distonion  Function,  Plane  8 
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Figiirc  4.22.  Profiles  of  UTf-AJ„f,  Plane  1 1 
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Figure  4.24.  Profiles  of  Utc/U„„  Plane  12 
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A— WJD  profile  data  (thin  bl) 
profile  data 

Figure  4.25.  Profiles  of  Distortion  Function,  Plane  12 
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Figure  4.30.  Profiles  of  Vtc./U_  Plane  8 
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Figure  4.32.  Profiles  of  Plane  5 
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Figure  4.35.  Pronie.s  of  WTc/U„f,  Plane  8 
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Figure  4.35.  (cont.).  rvofiles  of  WTc/U,ef'  Plane  8 
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Figure  4.36.  Profile.s  of  WV-AJ,,,,  FI;me  9 
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Figure  4.3-7.  (cont.).  Profiles  of  Plane  10 
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Figure  4.41.  Profiles  of  (Xj-r,  Plane  7 
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Figure  4.42.  (cont.).  Profiles  of  aj-c,  Plane  8 
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Figure  4.43.  Profiles  of  Ore,  Plane  9 
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Figure  4.54.  Profiles  of  u\c/\]^.  Plane  13 
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Figure  4.58.  Profiles  of  v'^cAJ^,  Plane  8 
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Figure  4,58.  (cont.).  Profiles  of  vVc/Uref,  Plane  8 
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Figure  4.59.  Profiles  of  v't-cAJ..,,  Plane  10 
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Figure  4.65.  (cont.).  Profiles  of  >v'TCAJref>  Plane  10 
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PiguiC  ...66.  Profiles  of  Plane  11 
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Figure  4,69.  FYofucs  of  10 
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Figure  4.70,  IVufiles  of  Plane  .S 
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Figiire  4.75.  (com.).  Profiles  of  Plane  10 
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Figure  4.78.  Boundaiy  Layer  Hodograph  Plots,  Plane  6 
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Figure  4.79.  Boundary  Layer  Hodograph  Plots,  Plane  7 
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Figure  4.80.  Boundary  Layer  Ho^dogranh  Plots,  Plane  8 


o’ao 


^Fs/'^re#  ^Fs/^ref  '^Fs/ 


X/C  =  0.18  (plane  5) 

C-JLF  profile  data 
•“SCO  profile  data 
^-WJD  LDV  profile  data 
C— JS  profile  data 


Figure  4.84.  Profiles  of  0(W)/aY)5/U^.  Plane  5 
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Figure  4.87.  Profiles  of  0(W)/aY)5/U^f,  Plane  8 
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Figure  4.87.  (cont,).  Profiles  of  0(W)/aY)6/U„f,  Plane  8 
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Figure  4.90.  Profiles  of  (3(W)/3Y)6/XJ,-f,  Plane  1 1 
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Figure  5.6.  Comparison  of  Contours  of  Distortion  Function,  f^,  Plane  5 
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Figure  5.21.  Comparison  of  Secondary  Velocity  Vectors,  Plane  6  (JLF  data  presents 
only  W  component  of  vectors) 
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Figure  5.23.  Conuparison  of  Secondary  Veiocity  Vectors,  Plane  8 
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Figure  5.25,  Comparison  of  Secondary  Velocity  Vectors,  Plane  10 
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Figure  5.27.  Comparison  of  Contours  of  (9(W)/3Y)6AJ„f,  Plane  5 
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a)  SCD  results 


b)  JLF  results 


Figure  5.29.  Comparison  of  Contours  of  O(W)/0Y)6/U„f,  Plane  7 
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a)  SCD  results 


b)  JLF  results 


Figuie  5.41,  Comparison  of  Contours  of  u  jr/U^,  Plane  7 
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Figure  6,1a.  Paneling  Schemes  for  NACA  Airfoils,  132  Panels 
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Figure  6.1e.  Shin  Paneling  Schemes,  1.32  Panels 
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Figure  6.3a,b.  NACA  Appendage  Velocity  Distribution  Results  0®  AoA 
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Figxire  6.6g.  Mehta  Appendage  Velcx-ity  Results,  12°  AoA 
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Figure.  6.6h,i  Ba.sclin.';  imd  Teai^iiop  ApjKndagc  Velocity  Dis’^ributions,  12"  AoA 
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Fig'-ire  6,7a,  Average  Voitex  Sa’ctchin^  vs.  Bluntncss  Facior,  0°  AoA,  Log-Log 
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Figure  6.7b.  Average  Vortex  Strctchuig  Rate  vs.  Bluntess  Factor,  0“  AoA,  Linear 
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Figure  6.7c.  Average  Vortex  Stretching  Rate  vs.  Blumness  Factor,  12°  AoA,  Log-Log 


Figure  6.7(1.  Average  Vortex  Stretching  Rate  vs.  Bluntess  Factor,  12°  AoA,  Linear 
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Figure  6.10.  Tx  vs.  Vonex  Transport  Pttrameter  for  Mehra’s  and  Shin’s  Appetdag 


Appendix  A.  Hot-Wire  Data  Redaction 


As  px'eviously  mentioned  in  section  2.5,  the  calculation  of  Us,  m's>  ^  s>  '^’s* 

Og  is  taken  from  Devenpcit  and  Simpson  (1990a).  Note  that  in  the  following  equations, 
an  overbar  designates  a  time  averaged  quantity.  Using  equations  3.2,  3.3,  and  3.4,  f(E0 
can  be  expanded  in  a  Taylor  se  'os  ./cout  E  (mean  voltage): 

U'  =  f(£')  =  f®  .  ...  (A.1) 

21 

Dropping  terms  higher  than  second  order  and  taking  the  time  average  of  this  equation 
yields: 

=  f(E)  +  (A.2) 

Calculation  of  u'^  follows  along  similiar  lines: 

^  (A.3) 
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Squaring  equation  A.3  and  taking  the  time  mean  average,  the  result  is; 


(A.4) 

Note  that  these  equations  are  in  terms  of  mean  voltages  only.  The  mean  fluctuating 
voltages  can  be  expressed  in  terms  of  mean  powers  of  voltages.  Tlie  equations  needed 
are: 

E^  -  E'^  - 

^  =  F  -  ‘ZFe  +  2^  (A.5) 

=  E^  -  4E^E  -  6E^^  “ 

Utilizing  these  equations  enables  the  reductiu..  >£  ’J, ad  u  by  looping  through  the 
array  of  voltages  only  once. 
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Appendix  B,  Wing-Bodj  Junction  Approach  Flow  Data 

This  appendix  cfiniains  listings  of  the  approach  flow  boundary  layer  data  for  the 
data  sets  of  Devenport  and  Simpson  {i990a),  Dickinson  (i986a,b)  and  McMahon  et  ah, 
(1987).  The  approach  flow  for  Devenport  fm.d  Simpson’s  case  is  the  same  as  for  the 
author’s  hot-wue  data. 

The  author's  daui  data  is  available  on  magnetic  disk,  and  includes  all 
racasurements  taken  hi  pkuies  3-1  i.  !  3~13  in  both  tunnel  and  ffecstream  coordinates.  The 
text  file  "FILET JST.DOC"  contains  infv  "ation  regarding  the  organization  of  the  data 
files  and  their  contents,  lliis  file  is  shown  in  this  appendix  as  well.  The  other  data  sets, 
in  similiar  format,  are  available  on  request.  Tnis  data  is  archived  using  PKZIP  v.  1.02. 
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Devenport.  and  SimpPon  thick  approach  boundary  layer  data 
(Fleming's  approach  boundary  layer  is  the  same) 

X/T  =  -2.116,  Z/T  -  0.0 
wing  thickness  ■“  7.17  cm 
nominal  reference  velocity  =■  26.38  m/s 
temperature  =  23  C  Pre.csTiie  =  940  mb 
nurriber  of  points  =<  30 


Y/r 

U/'Uref 

. 0035411 

.39103 

. 0042493 

.41412 

.0049575 

.42840 

.0063739 

.45372 

. 0074363 

.47239 

.0092066 

.49328 

.011331 

.50678 

.013810 

.52328 

.Gr6997 

.54b9.‘:5 

,020538 

.55712 

.025142 

.57837 

.03080'^ 

.59453 

.03  ’181 

.61043 

. 045326 

.  6  3 '  J  n  0 

.055241 

.  64853 

.067635 

.66889 

.082153 

68786 

.  10021 

.71183 

.  12181 

-s  n  >- 

.14837 

.7  62 

. 18945 

.8021, 

.21990 

.83043 

.26806 

.86820 

.32613 

.90733 

.39766 

.94622 

.48336 

.96978 

. 58853 

.  97587 

.71671 

.97498 

.  87252 

.  97695 

1.0623 

.  97722 

# 
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Dickinson  approach  boundary  layer  data 
X/T  -  -3,1667,  Z/T  -  0.0 
wing  thickness  *  6.1  cm 
nominal  reference  velocity  =  30.7  5  m./ s 


number  of 

points  » 

Y/T 

U/Uref 

. 041667 

.6186 

.08333 

.  6879 

. 12500 

.7253 

,1667 

.7556 

.20833 

.7782 

.3125 

.8221 

.41667 

.  8570 

.  5208 

.8864 

.  6250 

.  9097 

.729167 

.9310 

.8333 

.  9465 

1.04167 

.  9677 

1.2500 

.  9765 

1 .4583 

.9780 

1.6667 

.  97‘'5 

1.875 

.  9778 

2.0833 

.  9781 

McMchon,  et  al.  approach  boundary  layer  data 

X/T  *  -2.087,  Z/T  -  0.0 

wing  thic.kness  =  7,1  cm 

nominal  reference  velocity  15.24  m/ s 

number  of  points  =-  24 


Y/T 

U/Uref 

0.051355 

0,533333 

0.066386 

0 . 568957 

0,083922 

0 . 599491 

0.101458 

0,629008 

0.118994 

0.659542 

0.136530 

0.667684 

0 .152813 

0 . 673800 

0 .1'72855 

0.709415 

A 

O  .  21  U  D  M  ^  ± 

0.715522 

0.237988 

0.731807 

0.274313 

0.765394 

0.309385 

0.782697 

0.344457 

0.784733 

0.380781 

0.809160 

0.450925 

0.858015 

0.519816 

0.879389 

0 . 606243 

0.909924 

0.695176 

0.931298 

0 .779098 

0.950636 

0 .868030 

0.959796 

0 . 956963 

0.968957 

1.042138 

0.972010 

1,112282 

0 . 974046 

1 .180000 

0.973000 
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J.  Fleiing's  data  files  are  organized  in  the  follOHing  tanner: 

fpX](??,dat  -  XI  indicates  the  I  plane  nuaher  (05“i5) 

-  ??  *  "tC  indicates  tunnel  coordinate  syste*  used  to  nt  data 
-  Me"  indicates  freestreai  coordinate  system  used  to  present  data 

data  file  contents; 

file  descriptive  header 
X/T  location  of  YZ  plane 
average  reference  velocity  ii/s) 
nuiber  of  profiles  leasured  in  YZ  plane 

no.  of  pts  for  Ist  profile  te«p(C)  pressurelab)  uref(B/s)  fs  floi*  angle 


(one  line  for  *ac^  ''•ofile) 


no.  of  pts  for  last  profile  te*p(C)  pressureiab)  ureKt/s)  fs  flow  angle 
Y/T  Z/T  U/Uref  u'/Uref  W/Uref  w'/Uref  uw/Uref^2  flow  angle 


etc.... 


Note:  The  far  nne  planes  contain  only  y/'Jref  and  u'-'Ure^  data 
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Appendix  C.  Vortex  Panel  Code  Data 


This  appendix  contains  listings  of  data  discussed  in  Chapter  5.  This  data  is  also 
available  on  disk.  The  listings  containing  the  data  from  the  Mehta  appendage  shapes  are 
"padded"  for  the  and  Ty  /(d(V s)/dS)g  data,  since  Fx  data  was  available  only  for  n  = 
1.5,  2.0,  and  3.0.  The  padding  was  done  so  rhat  the  data  could  be  read  into  a  plotring 
program  correctly. 


3:2  elliptic  nose  *1  0020  tail  (baseline  mnq) 


BF  d(V/Ure1)/d(S/C!o 

ave.  vort.  str. 

VTP 

rhole/C 

.319564  21.2476 

,392097 

1110.91 

.0725128 

at  12  degrees  aoa 
,319564  13.717.0 

.3313B2 

■7  96157 

.0725128 

S8B  sail  section 

BF  d(V/Uret)/ai5/C)o  ave.  vort.  str.  VlP  rhole/C 

,0494162  40.4319  .0787788  388.372  .0258028 

at  12  degrees  aoa, 

,0494182  12.6335  .0813323  -13.2782  .0258028 
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flehts  appendages  (5X  thick,  C  angle  of  attack) 


n 

BF  I 

[l(V/Uref)/d(S/L)o 

ave,  vort. 

str.  VTP 

Ganta 

G/(dV/dS)o 

rhole/C 

1.5 

.00279715 

1683,47 

.00465945 

.856895 

,076 

4.51456-5 

.000932762 

1.7 

.00508421 

695.745 

.008467? 

2.73959 

.0001 

l,0e-B 

.00168978 

1.9 

.00683234 

500.126 

.0110025 

4.03173 

.0001 

l.Oe-8 

.00226667 

1.9 

.00904611 

380.278 

.0144661 

5.39899 

.0001 

l.Oe-0 

.00299552 

2.0 

.0117591 

302.285 

.0186674 

6.96739 

.102 

3.3743e-4 

.0038866 

2.2 

.0167575 

211.107 

.0293085 

10.01652 

.0001 

l.Oe-e 

.00617673 

2.5 

,0332238 

145.409 

.0505372 

15.0588 

.0001 

l.Oe-B 

.0109828 

3.0 

,0668013 

99.564 

.0973200 

23,7873 

.12J 

i.2454e-: 

.02:7128 

3.5 

.109117 

79,5194 

.152306 

32.4570 

.0001 

l,0e-8 

,0352400 

llehta  appendages  (lOX  thick,  0  deg,  angle  of  attack) 


n 

8F  d(V/Uref)/d(S/C!3 

ave.  vort.  str. 

VTP 

Gaiia 

G/dV/dS 

rhole/C 

1.5 

.00177800 

997.312 

.00311492 

1.65820 

.076 

7.6205e-5 

,00122342 

1.7 

,00386479 

37B.020 

.00669791 

4.68615 

.0001 

l.Oe-8 

,00266154 

1,8 

,00557815 

265.324 

.00962287 

6.59602 

.0001 

I.Ob-8 

,00383466 

1.9 

.00784436 

190.335 

.0134593 

8.69654 

.0001 

l,0e-3 

,00538276 

i.O 

.0107324 

155.745 

.0182997 

10.9562 

.102 

6.55e-4 

,00735104 

2.2 

.0185901 

107,118 

.0312528 

15.8482 

.0001 

l.Oe-8 

.012687 

2.5 

.0360552 

72.9092 

.0592239 

23.8528 

.0001 

1.09-8 

,0244796 

3.0 

.0885042 

49.4740 

.  137765 

38.1035 

,124 

2,5064e-3 

.0541469 

3.5 

.155162 

39,3555 

.232173 

52.6940 

,0001 

l.Oe-3 

,0^:4297) 

nehta  appendages  (lOX  thick,  12  deg. 

angle  of  attack) 

n  BF 

d(V/yref)/d(S/C)o 

ave.  vort,  str,  VI? 

6an«a 

G/dV/dS 

r‘iol®/ 

2.0  .0107324 

14,0409 

.0203631 

100) 

.102 

l.Oe-8 

,00735104 

2.'  .0360552 

13.2463 

,0697550 

1000 

.0001 

l.Oe-8 

.0244756 

3,''  .2805042 

17.1085 

,157409 

1000 

.124 

l,0e-8 

,0541469 

3,5  ,1551.62 

16.1884 

.281178 

1000 

.0001 

l.Oe-3 

,0542971 

Sand  1850  section 

BF 

d(V/Uref )/d(S/C)Q  ave.  vort.  str. 

OTP 

rhole/C 

.0132594 

116.567  .0254399 

134.168 

.00904039 

at  12  degrees  AOA 

.0132594 

10.8357  .0265676 

5.47B18 

,00904039 
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NACft  appendages  (0  deg.  angle  of  attack) 


IC  BF  d(V/Uref)/d(S/C)o  ave.  vort.  str 
0006  ,0102643  254.638  ,0107964 

0009  ,0161669  142,179  .0292247 

0012  ,0286692  87.8263  !o502008 

0015  ,0452494  59,1695  .0761781 

0018  .0713588  42,7751  ,1152(0 

0024  .144320  25.8816  ,23r0 


VTP 

86.5650 

140.193 

220.685 

328.404 

464,314 

840.104 


-•hole/C 

.00525416 

.00859113 

.0142638 

.0210323 

.0309974 

.0547926 


MACfl  appendages  (12  deg  angle  of  attack) 


(i(V/llref )/d(S/C)o  ave.  vort,  str 
0006  .0102643  11.9686  .0213268 

0009  ,0161689  13.6685  .0314099 

0012  ,0286892  15.5696  .0510680 

0015  .0452494  13.1528  ,0803768 

0018  ,0713568  14,4878  115761 

0024  .144320  12.8868  ,'213215 


.  VTP 
-8.28003 
5.10029 
-8,52875 
-7.86633 
-2.13583 
-11.3552 


rhoIe/C 

.00525416 

.00859113 

.0142630 

.0210323 

.0309974 

.0547926 


0.  Shin  4:1  elliptic  appendage  (n  =  2.0) 


BF  d(V/Uref)/d(S/C)o 
•0429412  53.0657 

at  12  degrees  aoa 

ave.  vort, 
.0749548 

str.  vrp 
397.908 

rhole/C 

.023866 

•0429412  12.2378 

.0777543 

-9.91513 

.023866 

J.  Shin  4:1  Nedge  elliptic  appendage  (n  =  1.5) 

BF  d(V/lJref)/d(S/C)o  ave.  vort.  str.  VTP  rhnle/'' 
■»«2S8<27  885.7,8 


Circular  nose  appendage  (teardrop) 


BF 

(f(V/Urei)/diS/C)o 

ave.  vort.  str 

.  VTP 

rhoJe/C 

1.0709V  9.18305 

it  12  degrees  aoa 

1.13362 

73.6825 

.210591 

1.070B9 

0.66207 

1.09447 

-27,5322 

.210591 
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